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Abstract 


The  work  carried  out  under  this  contract  was  part  of  a  wide  ranging 
investigation  into  the  effects  of  particulate  contamination  on  the  performance 
of  aerospace  hydraulic  control  sub-systems.  This  phase  of  the  work  included 
the  devising  of  techniques  for  studying  this  problem  and  the  application  of 
these  techniques.  Amo-gst  the  subjects  studied  were  methods  of  particle  size 
analysis,  the  theoretical  prediction  of  the  efficiency  of  sampling  devices, 
the  development  cf  a  mathematical  modelling  technique  to  simulate  particle 
behaviour  in  complete  systems,  a  discussion  of  the  meaning  of  particle  size 
and  how  this  may  be  related  to  analytical  techniques  and  a  new  method  for 
characterising  the  efficiency  of  filters. 

A  brief  analysis  of  the  current  situation  with  regard  to  particulate 
contamination  of  hydraulic  systems  and  a  proposed  fr* ~ework  within  which 
work  should  continue  is  also  included. 
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Chapter  1 


Introduction 


PART  I 


CHAPTER  1 


Introduction 


1 . 1  Layout  of  Report 

This  report  is  an  account  of  an  investigation  of  some  of  the  techniaues 
used  to  study  the  concentration  and  behaviour  of  particulate  contamination  in 
hydraulic  flight  control  sub-systems  and  a  basic  investigation  of  some  of  the 
phenomena  associated  with  the  presence  of  these  contaminants.  The  introductory 
chapter  will  contain  a  brief  resume  of  the  situation  as  it  is  seen  to  exist  at 
present  followed  by  an  account  of  the  approach  the  investigators  have  taken  to  the 
problem. 


The  following  sections  of  the  report  are  then  devoted  to  particular  topics, 
these  sections  often  being  the  work  of  one  particular  investigator  or  group  of 
investigators.  The  final  section  of  the  report  is  an  account  of  some  of  the  tech¬ 
niques  whi ch  have  either  been  developed  in  the  course  of  this  work  or  are  considered 
of  utility  in  a  study  of  the  subject. 

The  personnel  involved  in  this  work  were  Principal  Investigators  - 
Dr.  R.  J.  Akers  and  Mr.  B.  Scarlett,  Associate  Investigators  -  Mr.  P.  J.  Lloyd, 

Mr.  J.  I.  T.  Stenhouse  and  Mr.  A.  S.  Ward,  with  the  invaluable  assistance  of 

Dr.  P.  Allen,  C.  R.  G.  Treasure,  Miss  C.  Forsling,  Mr,  I.  Sinclair, 

Mr.  C.  A.  Troakes,  Mr.  J.  Walker,  Mr.  R,  Buxton,  Miss  C.  Bcstock,  Mrs.  M.  Warren, 

Mr.  L,  Moore  and  Mr.  A.  Saffell,  all  of  the  Chemical  Engineering  Department, 
Loughborough  University  of  Technology,  England. 

The  appropriate  references  are  given  at  the  end  oC  each  section. 


1.2 


Historical  Background 


When  hydraulic  systems  were  first  introduced  into  aircraft  it  was  to  aid 
the  performance  of  utility  functions  such  as  undercarriage  retraction,  bomb  door 
opening  and  flap  deployment.  The  earliest  systems  were  simple  hydrostatic  circuits 
consisting  of  a  force  transmitting  piston  and  a  linear  actuator.  With  the  increased 
use  of  hydraulic  power  in  aircraft  the  hydrostatic  systems  were  superseded  by  multi¬ 
purpose  systems  consisting  of  main  distribution  and  return  manifolds  pressurised  by 
a  motor  driven  pump  and  having  selector  valves  to  control  the  flow  of  fluid  to 
linear  and  rotary  actuators.  As  a  general  rule  the  actuators  were  only  required  to 
be  in  one  of  two  positions,  i.e.,  fully  extended  or  retracted,  hence  the  selector 
valves  needed  only  to  be  of  the  on-off  type.  If  more  than  two  positions  v/ere  neede4 
as  in  flap  operation,  a  double  actuator  was  used.  As  such  systems  were  only 
operated  intermittent ly  and  because  considerable  force,  applied  either  manually  or 
electromagnetically  was  available  to  operate  the  selector  valves,  catastrophic 
failures  were  relatively  few  and  wear  rates,  except  in  the  case  of  continuously 
operated  pumps,  low. 

The  development  of  powered  flight  control  systems  has  meant  the  need  for 
continuously  operating  hydraulic  systems  with  a  very  high  degree  of  reliability  and 
being  capable  of  moving  control  surfaces  with  precision  in  response  to  electrical  or 
mechanical  analogue  signals.  In  order  to  achieve  this,  servo  valves  of  very  great 
mechanical  sophistication  are  employed  and  it  is  these  servo  valves  in  particular, 
and  the  general  reduction  of  available  forces  and  mechanical  clearances  that  have 
made  hydraulic  systems  less  tolerant  to  particulate  contamination. 


1.3 


The  Present  Situation 


Discussions  on  the  effects  of  particulate  contamination  in  hydraulic 
systems  have  been  held  with  hydraulic  component  and  air  frame  manufacturers  and 
with  U.S.A.F.  ir.  their  roles  as  aircraft  users  and  originators  of  specifications 
(see  1-5).  The  picture  that  emerged  showed  much  confusion.  The  airframe  manu¬ 
facturers  considered  hydraulic  contamination  to  be  a  serious  problem  if  only  in 
what  it  cost  both  in  terms  of  precautions  during  manufacture  and  elaboration  of 
the  hydraulic  system.  The  attitude  was  very  much  a  'play  safe'  one  and  it  was 
felt  that  more  precise  information  on  system  tolerances  would  be  of  considerable 
value.  It  is  well  recognised  that  some  vehicles  are  more  prone  tc  contamination 
problems  than  others,  although  the  difficulties  Inherent  in  reliable  sampling  and 
analysis  of  fluids  from  actual  systems  makes  it  difficult  to  draw  firm  conclusions 
as  to  the  best  choice  of  components  for  a  given  system  and  the  specification  of 
its  operating  cleanliness  conditions. 

Among  hydraulic  component  manufacture re.  particularly  servo  valve 
manufacturers,  three  attitudes  were  discemable. 

These  were: 

(i)  No  serious  problem  existed. 

(ii)  That  the  effects  of  contamination  could  be  'designed  out.' 

(iii)  That  a  serious  problem  did  exist. 

It  was  felt  however  that  these  attitudes  were  themselves  open  to 
considerable  interpretation.  At  one  servo  valve  manufacturer  belonging  to  the 
'no  serious  problem'  school  photographs  were  shown  of  servo  valve  first  stage 
filters  and  metering  orifices  blocked  with  fibrous  and  particulate  contaminant. 

This  was  construed  as  an  occasional  but  inevitable  result  of  operation  of  the 
system  rather  than  a  serious  problem  ir.  system  design  or  specification.  Several 
members  of  the  "design  it  out"  school  were  of  the  opinion  that  there  was  an  optimal 
level  of  dirt  in  a  system  and  ..hat  if  'cler.uer  operation'  was  insisted  upon,  valve 
leakage  rates  woutd  increase.  The  reliance  on  a  degree  of  contamination  to  reduce 
leakage  implies  that  the  build  up  of  silt  within  the  valve  with  its  effects  on  both 
valve  wear  and  performance  degradation.  No  studies  have  yet  been  made  on  the  effects 
of  'silt'  on  the  forces  necessary  to  'break  out'  a  spool,  but  whatever  their  magnitude 
they  would  add  something  to  the  hysteresis  of  the  valve.  With  respect  to  valve  second 
stage 


leakage,  no  direct  experimental  information  was  available  as  to  the  separate 
contributions  of  port  overlap  or  underlap  and  leakage  through  the  annular 
spool/sleeve  clearance.  A  recent  study  (AFFBL-TR-69-136)  has  shown  the 
effect  of  particulate  contamination  on  both  second  stage  port  dimensions, 
wear  of  the  flapper  plate  (fig.  1.1),  and  causing  respectively  an  increase 
in  port  overlap  and  jet/flapper  clearance. 

The  increasing  stringency  of  fluid  cleanliness  standards  will  impose 
now  demands  on  filter  manufacturers.  It  is  probable  that  woven  wire  filters 
are  getting  near  the  economic  and  technical  limit  in  the  production  of  fine 
woven  wire  cloth  and  that  the  change  towards  filters  of  the  depth  type 
using  non-woven  porous  media  will  continue.  The  loading/pressure  and 
maximum  dirt  holding  characteristics  of  these  filters  might  also  be  considered 


mors  attractive. 


1,4.  Problem  of  Specifying  Cleanliness  Standards 

In  specifying  rainio'na  fluid  cleanliness  standarc.s  for  a  given  system  by  any 
means  other  than  crbiti-ary  selection  f.t  Till  be  necessary  to  start  from  the  contam¬ 
ination  characteristics  of  those  components  whose  design  is  dictated  yy  the 
operational  requirements  of  the  system,  and  then  select  other  components  such  as 
filters  to  meet  their  cleanliness  specification.  To  Cc  this  much  information  is 
needed  about  the  contamination  characteristics  oi  a  wide  range  of  hydraulic  compo¬ 
nents.  included  in  this  would  need  to  be  a  study  of  the  rates  of  deposition  and 
rsmoval  of  particles  from  the  interior  surfaces  of  pipelines,  hydraulic  accumulators 
aid  other  components  and  the  effects  of  vibration  and  accelerating  forces  on  those 
proces'ves. 

The  operators  of  aircraft  would  like  a  simple  go/no  go  test  of  hydraulic 
cleanliness  but  until  the  order  oi  seriousness  cf  the  facte rs  involved,  e.g., 
particle  size  and  size  distribution,  particle  hardness,  etc.,  any  such  test  would  be 
based  upon  very  unsure  foundations.  In  particular,  with  tfcv  reduction  of  radial 
clearances  in  servo  valves  to  50  micro  inches  li.e.,'/'2  PK)  conside.-u.ble  attention 
must  be  given  to  studying  the  effects  of  par* teles  finer  than  those  counted  by 
existing  techniques.  Unfortunately  there  are  few  particle  size  analysis  techniques 
available  for  particles  <  1  pm  and  those  that  there  arc  vend  to  be  laborious. 


«53fvS^'  *^'1!!‘S**^~&$&tlS?&mt!e&9g&&s®&f&s&Qs&&&G®ti& 


1.5.  General  Effects  of  Particulate  Contamination 

These  may  be  classified  as: 

1.5.1.  Wear 

If  a  component  in  which  sliding  contact  occurs,  e.g.,  an  actuator,  is  d-‘  s- 
mar.tled  after  service,  longitudinal  scratches  will  be  seen  on  the  sliding  surfaces, 
those  being  due  to  particles  trapped  between  the  mobile  anJ  stationary  members. 
Although  wear  ot.  this  type  is  general  it  is  of  particular  significance  in  hydraulic 
piston  pumps,  where  it  mainly  occurs  between  the  valve  plate  the  cylinder  block 
and  between  the  slippers  and  swash  plate.  Less  extensive  wear  will  also  occur  in 
the  pistons  and  cylinder  boras.  The  overall  effect  of  this  wear  is  to  reduce  the 
efficiency  cf  the  pump  by  increasing  leakage  and  also,  more  importantly  from  the 
point  of  view  of  particle  contamination  cf  the  wbo?e  system,  generating  large 
numbers  of  fine  particles,  a*- "erai  studies  of  the  effect  ot  artificially  contaminated 
oils  on  pumps  have  shown  thee  to  be  very  effective  grinding  machines  with  the  pump 
output  containing  sice  reduced  contaminant  particles  and  aatal  from  the  pump  origi¬ 
nating  in  wnai  caused  by  contaminant  passing  tn rough  the  pump.  Wear  of  this  type 
betwee-  sliding  surfaces  is  also  ooserved  in  actuators  where  it  is  not  a  major 
problem  except  in  those  cases  where  an  actuator  is  dithered  for  long  periods  of  time 
about  one  position,  and  in  the  spools  and  sleeves  of  servo  valves  where  it  leads  to 
increase  leakage  and  degradation  of  the  valves  performance  in  terms  c-f  such  para¬ 
meters  as  flow  again  and  hysteresis. 

1.5. ?.  Blockage  of  Orifices 

'pUis  is  a  possible  although  uncommon,  mode  of  failure.  The  fir.ast  orifices 
found  in  conventional  flight  •  xjntrol  systems  arc  the  metering  orifices  in  the  first 
®d  second  stages  of  servo  valves.  However  these  orifices  are  almost  invariably 
protected  by  filters  built  into  the  valve.  It  is  well  known  that  contamination  by 
fibres  is  very  important  in  leading  up  to  the  blockage  of  an  orifice  by  a  bed  of 
particles  aach  of  which  individually  is  smfJLler  than  the  orifice  and  this  type  of 
blockage  cannot  be  discounted  even  in  filter  protected  situations. 
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1.5.3  Sifting 

Tests  of  electro-hydraulic  servo  valves  operating  with  contaminated  oil 
shows  silting  to  be  one  of  the  major  effects  produced.  Silting  manifests  itself 
as  mechanical  binding  between  two  surfaces,  e.g.,  spooJs  and  sleeve;  leading  to 
reduced  flow  gain  and  increased  hysteresis.  The  mechanism  of  silting  is  not 
very  well  unders~ocd.  if  i*:  is  due  to  particles  entering  the  annular  space 
between  the  spc-c  i  act",  sleeve  the  particles  concerned  must  be  smaller  (in  their 
smallest  dimension}  than  the  diametrical  clearance.  When  a  granular  material 
flows  it  has  to  rearrange  its  packing  in  order  to  allow  the  particles  to  move 
freely  over  one  another.  This  bringc  about  a  reduction  in  density,  i.e.,  the 
volumetric  "dilation”  of  Osborne  Reynolds.  In  the  spool/? leeve  situation  this 
dilation  cannot  occur  due  to  the  constrained  geometry.  Very  little  is  known  of 
particle  flow  under  these  conditions  of  fixed  volume. 

Alternate  explanation  cf  silting  is  an  off  axis  tiiting  of  the  spool  and 
subsequent  binding  due  to  particles  jammed  under  one  side,  although  calculations 
of  the  maximum  angular  deflection  psssible  in  a  modern  valve  make  this  unlikely, 
cr  alternatively  the  build  up  of  a  ,:*'illet”  or  particles  between  the  faces  of 
the  spool  and  the  sleeve  wall  which  would  tend  tv  oecome  jammed  in  the  annulus 
when  the  spool  moves.  This  latter  mechanism  is  more  likely  and  would  probably 
occur  along  with  oarticle  build  up  in  the  annular  space. 

1.5.4  Filter  Blinding 

Because  it  is  necessary  to  control  the  number  of  particles  circulating 
in  a  hydraulic  circuit,  filters  are  introduced  into  the  system.  As  a  con¬ 
sequence  of  their  action  filters  gradually  become  ’blinded’  with  particles  and 
have  to  be  replaced  or  deaned.  The  design  of  a  filter  is  an  attempt  to  optimise 
the  performance  in  terms  of  the  fineness  of  particle  retained  with  the  longest 
time  between  services. 
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1.6. 


Aims  &  Methods  of  the  Investigation 


The  object  of  this  programme  is  to  study  the  effects  of  particulate  contam¬ 
ination  on  hydraulic  flight  control  sub-systems  and  ii  possible  define  design 
criteria  and  contamination  levels  to  enhance  che  life  and  reliability  of  such 
systems.  From  previous  studies  it  is  known  that  the  reliability  of  systems  varies 
widely  from  one  design  of  vehicle  to  another  and  also  between  different  examples  of 
the  same  venicle.  These  studies  have  tended  to  be  empirical  in  that  they  have 
attempted  tc  correlate  the  probability  of  the  various  modes  of  failure  with  the 
contamination  level  present  at  one  or  seve.al  points  ..n  tr.e  hydraulic  system. 


This  project  is  an  attempt  tc  apply  analytic:.!  techniques  to  this  problem. 
This  has  been  done  by  constructing  a  mathematical  i;odel  of  the  system  (Chapter  3) 
in  terms  of  the  characteristics  of  the  individual  components  and  the  particle  fluxes 
present.  This  approach  is  illustrated  in  fig.  2  ,  where  it  is  seen  that  the 
ultimate  test  of  the  model  direct  ccmparisor  with  a  real,  if  simplified, 
hydraulic  system. 


In  order  to  use  this  model  it  is  necessary  to  characterise  the  individual 
components  of  the  system.  By  characterise  we  mean  the  number  and  size  of  particles 
leaving  the  component  as  a  function  of  the  number  and  size  of  particles  entering - 
the  flow  rate  and  any  other  relevant  operating  parameters  ,  e.g.,  dirt  loading 
present  on  a  filter.  The  modelling  of  the  individual  r  -nponents  nay  be  either  by 
experimental  study,  theoretical  analysis  or  a  combination  *>f  the  two. 


Whichever  of  these  two  approaches  is  followed  it  is  necessary  to  he  able  to 
describe  particles  in  a  way  *hat  can  be  used  subsequently  in  an  analytical 
description  of  their  behaviour.  Traditionally  the  size  of  a  particle  is  expressed 
in  terms  of  a  sphere  equivalent  in  a  property  such  as  area  (e.g.,  permeability' 
volume  (e.g.,  ''Coulter"  Counter)  projected  area  (e.g.,  Hiac'  Counter}  or 
sedimentation  velocity.  Alternatively  a  statistical  diameter  is  used.  We? 1  knewr 

A 

examples  are  those  of  Martin  and  Feret,  and  of  particular  interest  in  hydraulic 
engineering  is  the  maximum  projected  linear  dimension  used  in  ARP  Z38  and  similar 
methods.  These  various  diantters  are  functions  of  the  particle  chart  and  are  cnly 
equal  for  spheres.  However,  because  so  much  Information  is  lose,  when  an  irregular 
particle  is  defined  in  terms  cf  one  parameter,  alternative  descriptions  are  y.eeded 
if  properties  of  assemblies  of  particles  are  to  be  deduced.  Such  a  method  :.s 
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proposed  in  Chapter  4  and  an  application  of  it  is  shown  in  Chapter  6.  To  exploit 
this  technique  as  a  practical  method  further  dev  iupments  are  required  in  particle 
size  analysis  and  these  are  described  in  Chapter  5.  One  of  the  consequences  of  such 
a  method  is  that  from  the  description  of  the  particle  system  ao  obtained  the  particle 
size  distribution  in  term-:  of  any  of  the  more  familiar  techniques  of  characterisation 
should  be  available.  These  relationships  are  being  investigated. 

Fig.  2  shows  how  the  definition  of  particle  size  is  related  to  the 
characterisation  of  components  through  both  theoretical  and  experimental  techniques 
and  does  m  fact  provide  a  link  between  them. 

In  -rdor  to  carry  out  meaningful  particle  size  analyses,  whether  during  the 
characterisation  of  individual  components  or  in  the  study  of  complete  systems,  it  is 
necessary  to  be  able  to  take  a  representative  sample  from  the  system  or  some  point 
in  the  system.  This  is  discussed  in  Chapter  2  which  describes  mainly  theoretical 
studies  of  the  flow  of  suspensions  through  sampling  devices.  Experimental  work  to 
check  these  predictions  is  in  progress. 
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2.1  Introduction 

In  order  to  study  the  effect  of  particles  on  the  behaviour  of  hydraulic 
systems  or  to  establish  the  degree  of  contamination  in  operating  systems  it  is 
necessary  to  be  able  to  remove  representative  samples  of  fluid  from  the  sysien,  or, 
alternatively,  to  insert  a  contamination  measuring  device  into  the  system.  The 
latter  will  not  be  considered  here. 

Any  sampling  device  has  to  satisfy  several  requirements.  These  are: 

(a)  That  the  sample  should  be  removed  from  the  system  at  such  a  point  that  it 
indicates  the  state  of  fluid  arriving  at  a  contamination  sensitive  element 
in  the  system  or  the  state  of  the  fluid  leaving  an  element  known  to  give 
rise  to  significant  amounts  of  contamination.  As  examples  of  these  are 
servo  valves  and  pumps  respectively.  This  condition  is  complicated  by 
the  fact  that  some  system  components  are  kncwn  to  accumulate  panicles  and 
then  release  them  in  a  short  period  Gf  time.  An  example  of  this  type  of 
device  is  a  hydraulic  accumulator. 

In  any  hydraulic  system  -.lie  rate  of  fluid  flow  at  any  given  point  varies 
widely  and  undergoes  very  rapid  accelerations.  As  a  simplification  the 
flow  rate  may  be  considered  as  being  comprised  of  leakage  with  pulses  super¬ 
imposed  on  it.  An  ideal  sampling  device  would  integrate  the  cont-nunation 
over  a  period  of  time.  This  could  be  done  by  withdrawing  a  repi  esentati ve 
sample  over  a  considerable  period  of  time,  although  this  would  be  unaccept¬ 
able  in  flight  control  systems  on  account  of  tne  volume  of  fluid  withdrawn, 
particularly  in  military  aircraft  or  missile  systems  with  their  small 
reservoir  capacities.  An  alternative  to  this  would  be  the  installation  of 
an  on  stream  monitor.  Only  one  particle  counting  device  capable  of  being 
installed  in  an  operating  system  is  currently  available,  the  Hiac  Counter, 
and  this  has  the  disadvantage  of  being  sensitive  to  flow  rate  (see  5.3). 
Alternatively  a  sample  could  be  withdrawn  from  the  system  at  a  particular 
point  and  from  a  knowledge  of  the  system  as  a  whole  its  significan>  ■  inter¬ 
preted  in  terms  of  the  probability  of  performance  degradation  or  catastrophic 
failure.  To  do  this  requires  a  mathematical  representation  of  the  system  as 
described  in  Chapter  3. 
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t  y)  y.vit  the  sample  piesr.nted  to  the  sampling  device  be  representative  of  the 
fluid  flowing  in  that  particular  line.  It  is  known  tnat  when  a  dilute 
suspension  of  non  deformable  pa  r tides  flo^s  along  a  cylindrical  tube  at 
flow  rates  corresponding  to  Reynolds  Numbers  in  the  laniaa  flow  region,  a 
radial  concentration  gradient  is  set  up  with  the  particles  migrating  towards 
the  tube  wall.  Conversely  de for  .Table  particles  (e.g.  water  droplets)  will 
migrate  towards  the  axis  (Goldsmith,  H.L.  and  Mason,  S.G.  in  Rheology,  Voi. 

4,  Ed.  Eirich,  Academic  Press  1967). 

To  satisfy  this  requirement  either  the  flow  velocity  must  be  great  enough  to 
cause  sufficient  turbulence  to  ensure  efficient  mixing  or  turbulence  must  be 
artificially  inducec  just  upstream  oi  the  sampling  point.  By  normal  design 
criteria  the  Reynolds  Number  in  a  hydraulic  pipe-line  has  a  maximum  value  of 
approximately  7C0G,  with  minimum  values  corresponding  to  leakage  flows  for 
servo  valves  being  in  the  range  200-800.  As  the  high  flow  rates  c-n  only  be 
maintained  (in  systems  containing  'inear  actuators)  for  short  periods  of  time 
a  realistic  Sa?ipling  device  sample  frc-.v,  a  system  flowing  at  the  minimum  rote 

hence  requiring  the  artificial  induction  of  «.ui  oulence.  This  may  be  done 
by  either  changing  the  sh^pe  of  the  flow  channel,  as  in  the  Pi osser  Indust¬ 
ries  .candling  valve  or  vy  placing  structures  within  the  existing  cylindrical 
pipe.  \¥hat‘”,er  the  design  of  turbulence  incucot  it  must  create  the  maximum 
turbulence  with  minimum  pressure  drop.  Also  the  possibilities  of  erosion 
within  t?.a  device  due  to  turbu lence  must  be  considered. 

(c)  That  the  sampling  device  is  capable  of  removing  from  the  syston  the  re  - 
presentative  fluid  as  presented  to  it  (section  2.2.). 

(d)  That  the  sampling  device  does  not  generate  or  retain  contaminant. 

Generation  is  most  likely  to  happen  in  the  valve  portion  where  sliding 
surfaces  nay  be  present.  It  is  desirable  that  valves  bo  designed  to  prevent 
this  kind  of  motion  occurring,  e.g.  the  sampling  valve  manufactured  by 
fairey  Hydraulics,  U.K. 

(c)  That  the  sampling  device  does  not  impair  the  safety  or  reliability  of  the 
system. 
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2.2  ^heury  of  Sampling  Devices 

In  any  sampling  system  the  fluid  flow  is  not  planar.  Under  non  planar 
conditions  the  solid  particles  due  to  their  inertia  will  cross  the  fluid  stream¬ 
lines  such  that  any  fluid  withdrawn  will  not  be  representative  of  the  fluid  as  a 
whole.  The  errors  due  to  this  mechanism  can  be  estimated  theoretically.  The 
procedure  employed  is  to  calculate  the  trajectories  of  particles  of  different 
inertia  levels  in  the  flow  field  pertaining  under  the  particular  sampling  condit¬ 
ions.  A  pre-requisite  cf  such  a  theory  is  an  adequate  description  o*  the  fluid 
flow  field. 

Variations  on  two  basic  geometries  have  been  used  in  sampling  valves. 

These  are: 

(i)  An  open  aperture  in  the  side  wall  of  a  tube  (fig.  3  )  or  in  the  case  of 
the  Prosser  Industries  valve  (fig.  4  )  within  a  chamber  containing  turb¬ 
ulent  flow. 

(ii)  A  probe,  usually  central,  within  the  tube  carrying  the  flow.  Various 
shapes  of  probe  have  been  used,  e.g.  pointing  upstream  (fig,  5  1,  per 
pendrcular  to  the  stream  (fig.  6  )  perpendicular  to  the  stream  and  cut 
obliquely  (fig.  7  ) 

Prom  the  point  of  view  of  their  flow  fielcs  these  latter  two  designs 
p-obably  incorporate  characteristics  of  both  (i)  and  (ii)  above. 

Most  of  the  flow  situations  encountered  in  sampling  valves  are  not  readily 
amenable  tc  analytical  description  so  the  flow  field  in  such  components  must  be 
determiv.ed  using  flow  visualisation  technique*  r.nG  described  empirically.  A  flow 
visualisation  technique  has  beer,  de’^luped  and  is  currently  being  used  to  invest¬ 
igate  the  flow  fields  uruer  a  ru rge  of  conditions.  This  technique  is  described 
in  2.3.1. 

Th?  theoretical  technique  has  been  developed  using  i*n  ideal  model  of  the 
flow  field  ■*-»  a  wimple  sampling  va.ve.  This  technique  and  the  results  which  have 
bran  achieved  -sing  it  are  described  in  2.2.1  below. 

Jince  the  particle  is  moving  in  the  linear  Stokes  regime  equation  (?)  tay  be 
replaced  by  t ue  more  familiar  trajectory  equations: 
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d  y  t  dy  =  v 
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where  u  and  v  are  the  x  and  y  resolutes  of  the  fluid  velocity. 

Equations(2.3)  may  be  reduced  to  a  dimensionless  form  by  using  the 
conversions  in  table  1. 


TABLE  I 


Notation 


i  Parameter 

Derivation 

t 

i 

Velocity 

u 

•■i/bo 

V 

v/Uo 

Co-ordi  nato 

X 

x/1 

y/i 

Time 

T 

Uo  T/i 

I  Stokes  Number 

L 

GO 

L _ 

rp2  p  Uo/9pl 

The  Stokes  number  is  *.ne  ratio  of  the  distance  a  particle  rith  initial 
velocity,  Uo,  will  tra-ol  before  coming  to  rest,  to  some  dimension  3,  of  the 
system.  In  the  case  of  a  tube  withdrawing  a  sample  from  a  streas,:  in  planar  flew, 
as  m  fig.  3  ,1  is  the  radius  of  the  sampling  tube.  Hence  equation  (2.3)  becomes: 


2  Stk  d_X  r  dX  =  U 
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Expressions  for  U  and  V  in  terms  of  X  and  Y  are  ..-ormally  complex  necess¬ 
itating  a  c'wputer  solutior.  for  equation  (2.J).  For  this  purpose  the  equations 
are  expi  ;ssed  in  finite  difference  form  ignoring  terms  higher  Than  those  of  u 
secona  order  in  a  Taylor  expansion.  F.-pt essior.s  5  thus  derived,  are  used  in  a 
stop  by  step  calculation  of  particie  trajectories: 
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2,2.1  Theoretical  Estimation  of  Sampling  Efficiency  of  a  Tube  Probe 

The  sampling  valve  under  consideration  is  shown  diagramatically  in  fig. 

8. 


Central  Probe  Sampling  Valve 

It  consists  of  a  sampling  probe  which  is  inserted  in  the  main  pipe  in  which  the 
fluid  is  flowing  with  a  velocity  UQ.  In  this  system  sampling  will  be  perfect  if 
the  sampling  velocity  is  identical  to  the  main  stream  velocity.  However,  if  the 
sampling  is  anisokinetic  then  the  flow  will  not  be  planar  and  sampling  eriors  will 
arise  due  to  inertial  forces  acting  on  the  particles. 

2 . 2 . 1 . a  Description  of  Flow  Fie  la 

The  flow  field  may  be  approximated  to  a  point  source  in  a  planar  flow. 

The  field  will  be  valid  for  small  values  of  1  and  high  Reynolds  Numbers.  The 
cartesian  velocity  resolutes  are  given  by;- 
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Q  represents  the  p lr'!,:£>.th  of  the  point  sou.  '  2.  Obviously  for  isokinetic 
sampling  conditions  it  is  zero.  The  relationship  between  Q  and  tJ-sS  relative 
sampling  velocity,  V  ,  is  given  byr 


V?.  -  1  + 


2-nd 


VR  is  the  ratio  of  the  mean  velocity  in  the  sampling  piohs  to  the  main  stream 
velocity. 


2 .2 . 1 . b  Calculation  of  Particle  Trajectories 

The  particle  is  considered  as  a  pc.nl  mass,  which  does  not  effect,  the 
flow  field  in  Jie  proximity  of  the  sampling  probe.  The  critical  trajectory, 
for  which  the  particle  just  touches  the  edge  of  th.e  probe  is  computed. 

A  force  balance  or.  a  particle  in  a  fluid  results  an  the  following  vector 
equation: - 


a  du  =  C  T  r  cOt  /  0  -  u\‘ 

—  u  >  ;  ~~  - _  i 


dr- 
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2.2. l.c  Trajectory  Calculation 

It  is  assumed  that  a  particle  approaches  in  the  x  direction  with  the  main 
stream  voli-city  to  a  point  i2  upsti-eam  from  the  probe.  Equations2_5  are  then 
used  to  calculate  the  path  taken  by  the  particle  until  it  is  level  with  the  mouth 
of  the  probe.  New  starting  pointc  on  the  Y  axis  are  taken  for  the  trajectory 
until  the  critica1  trajectory  is  four-d.  The  distance  upstream  of  the  starting 
point  and  the  step  length  'n  the  cal -.elation  have  been  chosen  in  such  a  way  that 
the  calculation  errors  resultin'  iron  the  use  of  a  fmito  different  technique  are 
not  significant.  Tito  calculation  was  performed  using  a  1905  ICT  digital  computer. 

2. 2.1.1  Results  and  Discussion 

The  efficient t  of  sampling  can  conveniently  be  defined  by 


concentrat-o-i  or  c^ntacunant  in  sample 
concentration  of  contaminant  in  main  stream 


The  calculated  efficiency  is  shown  as  a  function  of  V  ,  the  relative 

K 

velocity,  and  St’*.,  the  particle  Stokes  No,  in  fig-  ^  It  is  obvious  that  as  the 

Stokes  No.  tends  to  zero,  —  must  tend  to  1.0.  Again  as  Stk  ,  then 

co 


Agreement  with  these  Uniting  values  is  shown  in  fig.  2.7  which  shows  that 
the  transition  range  of  Sxoke^'  Nos.  for  this  kind  of  sampler  is  0.1  -  10.0. 

On  applying  the  data  :n  fig.  9  to  condition*  in  a  hydraulic  system  with 
fluid  viscosity  2CoS  flowing  through  a  b"  boro  pipe  at  flow  rates  corresponding 
to  fio  ~  200  CO. 6  g.p.m.)  P.nd  Re  =  6,000  {18.2  g.p.v*.),  a  sampling  orifice  1/8”  dia. 
containing  splioricll  aluminium,  copper  end  iungr carbide  particles  of  30V*t 
50  vmand  2C0ji*the  following  sampling  efficiencies  are  found  for  =  2.0  and 

V  =  °*5  ^i‘able  II  >- 
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Efficiency 
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TAJ5LE  II  Calculated  efficiency  of  central  probe  sampler 
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Re  = 

200 

Re  =  P 

2.0 

0.5 

2.0 

0.98<E<1.00 

1.02>E>1.00 

0.9«<E<1.00 

0.98 <£<1.00 

1*02 >£>1.00 

0.95<E«:1.00 

0.95<E<1.00 

1.05 >£>1.00 

0.72 

0.98<E<1.00 

1.02>E>1.00 

0.98 <E<1.00 

0.98<E<1.00 

1.02>E>1.00 

0.85 

0.90 

1.17 

0.5S 

0.98  E  1.00 

1.02>E- 1.00 

0.98<E<1.00 

0.98  E  1.00 

1 

0.80 

i  0.87 

1.21 

0.54 

These  results  show  that  sampling  in  the  range  0.5  <VR<  2.00  from  a  system 
at  leakage  flow  rates  would  be  satisfactory  for  particles  <  100  V  whatever  their 
composition  and  hence  density.  Hence  one  possible  approach  to  sampling  would  be 
to  have  a  turbulent  flow  valve  throttled  for  approximately  isokinetic  operation 
Ci  e.  0.5  <  V  <  2.0)  at  leakage  flow  just  upstream  of  the  sensitive  comoonent. 


2.3  ^-perimental  Jtudy  of  Sampling  Systems  -  Determination  of  Flow  Field 

Flow  visualisation  techuicues  are  used  to  determine  the  flow  fields  in 
models  of  sampling  valves,  the  fields  in  which  are  not  amenable  to  analytical 
description.  The  technique  used  is  that  of  ~treak  photography  of  neutral  density 
particles. 

Figure  10  shows  a  diagram  of  the  apparatus.  Light  from  an  iodide-quartt 
lamp  is  collimated  into  a  two-dimensional  (1/83  thick)  bet t  ;o  illuminate  a  plane 
within  the  sampling  valve.  The  valves  are  transparent  models  of  the  sampling 
valves  encountered  in  practice.  To  avoid  problems  of  parallax  m  observing  the 
flow  inside  the  valves  the  model  is  enclosed  in  a  box  which  contains  liquid  which 
has  the  same  refractive  index  as  the  valve.  In  this  way  a  flat  smooth  surface  is 
presented  fu;  photography. 

The  liquid  (water)  being  u:»«d  iontair.3  a  very  low  concentration  of  small 

neutral  density  particles  (polystyrene)  ohich  Will  follow  the  path  lines  of  the 

o 

fluid.  Streak  photographs  taaen  at  an  angle  of  120  to  the  light  bear;  sl'.vy  the 
path  of  those  particles  clearly  enabling  the  complete  flow  field  to  be  oete'  mined.. 
?he  conditions  which  are  being  studied  are  dynamically  r>o  iar  to  tnnee  which  are 
usee,  iv  sampling  aircraft  hydraulic  systems.  This  technique  is  now  being  used  to 
study  flow  fields  in  a  number  of  sampling  valves  under  a  range  of  conditions  and 
the  full  results  will  be  reported  at  a  I'.ter  date. 

One  of  the  fi-*st  photographs  taken  whilst  making  a  feasibility  study  of 
the  system  is  included  in  tnis  report  as  fig.  II.  This  is  a  photograph  of  the 
flow  in  a  simple  T  valve.  The  path  lines  of  the  fluid  are  clearly  shown  by  the 
idln.  .treaxs,  tn  this  particular  test  IOC  Uicron  classified  particles  of  alum¬ 
inium  uust  wore  adeed  ~c  the  fluid.  feiese  a.O  c'-own  by  the  mo^e  intense  streaks. 
It  can  be  seen  that  these  particluo  which  have  high  inertia  cut  across  the  fluid 
streamlines .  When  this  occurs  in  practice-.,  the  sample  obtained  will  not  be 
representative. 
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-  Glass  model  of  sampling  valve 
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Box  containing  liquid  of  same  refractive 


Index  as  glass 


FIGURE  10  Apparatus  for  flow  visualisation 
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Mathematical  Modelling  of 
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3.1.  Introduction 

This  section  describes  an  attempt  to  model  a  simplified  aircraft  hydraulic 
actuating  .system,  so  that  given  sufficient  information  the  model  might  be  used  to 
calculate  the  probability  of  a  failure  occurring  within  the  system.  This  simplified 
model  system  is  not  normally  found  in  practice  but  however  it  does  illustrate  how 
the  wu'u«cat’c?l  techniques  could  be  employed  to  describe  any  real  system.  No 
actual  calculations  have  been  carried  out  on  the  model,  because  at  present 
sufficient  information  about  the  parameters  involved  is  not  available.  Methods  of 
obtaining  this  information  is  of  course  the  major  function  of  this  contract.  The 
model  used  is  based  upon  the  representation  of  particle  concentrations  by  vectors 
and  the  system  by  matrix  operations  upon  these  vectors.  The  model  will  be  complex 
and  is  therefore  most  useful  if  it  is  in  a  form  suitable  for  computer  solution. 

This  fact  was  considered  throughout  the  report  and  the  final  form  of  the  model  is 
a  ccmpcter  flow  sheet. 


3.2*  Aircraft  Hydraulic  Systems 

3.2.1.  A  hydraulic  system  is  a  means  of  transmitting  newer  which  makes  use  of 
pressurised  fluid.  The  system  must  consist  of  at  least  3  parts,  a  pressure  source, 
a  transmission  line,  and  an  actuator. 

3.2.2.  In  aircraft  systems  the  pressure  source  is  a  mechanical  pump  taking  its 
power  from  either  the  main  aircraft  engines  or  from  auxiliary  power  supplies^  The 
most  usual  type  of  pump  is  c  constant  pressure,  variable  delivery  pump.  However 
the  system  modelled  in  this  report  makes  use  of  a  constant  delivery  pump. 

3.2.3.  The  motor  transforms  the  work  available  in  the  moving  fluid  back  into 
mechanical  peirer.  The  simplest  ar.d  most  usual  form  of  motor  is  a  piston  actuator, 
which  provides  oscillatory  motion  only.  If  rotating  motion  is  required,  then  a 
motor  similar  in  design  to  the  pump  is  employed. 

3.2.4.  A  hydraulic  fluid  must  meet  many  requirements  amongst  which  are  reasonable 
viscosity  over  the  desired  temperature  range,  resistance  to  mechanical  and  chemical 
breakdown,  good  lubricity  characteristics,  compatibility  with  constructional  metals 
and  elastomer  components,  low  specific  gravity  and  reasonable  price. 

3.2.5.  The  basic  system  considered  in  this  report  is  shown  in  Figure  12 

3.2.3.  The  reservoir  is  necessary  to  accommodate  changes  in  the  volume  of  fluid 
within  the  circuit  and  to  provide  a  reservoir  of  fluid  to  replace  that  lost  by 
leakage . 

3.2.7.  A  filter  is  included  in  the  pressure  line  from  the  pump.  In  real  circuits 
there  are  several  filters,  but  almost  invariably  one  in  the  pump  pressure  line. 

C.2.8.  The  system  illustrated  has  a  pressure  relief  valve  in  conjunction  with  a 
constant  delivery  pump.  It  would  not  be  difficult  to  alter  the  model  to  accommodate 
a  constant  pressure  variable  delivery  pump  with  the  relief  valve  only  being  used  to 
relieve  excess  pressures  in  the  pumping  circuit. 

3.2.9.  The  system  is  equipped  with  a  4-way  valve  enabling  an  actuator  to  be 
operated  in  both  directions. 

3.2.10..  In  Figure  12  arrows  show  the  direction  of  ilow  of  the  fluid.  In  the 
system  illustrated  the  flow  is  unidirectional  in  all  compenents  before  the  4-way 
valve,  but  beyond  the  valve  is  in  either  direction  depending  upon  the  valve  position. 
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■"he  pert  of  the  system  before  tl  in  which  the  fluid  flows  in  one  direction 

only  is  -cal  ltd  power  oiruiii. 

The  oppiMtion  of  the  system  shown  in  Figure  li  may  be  understood  by 
studying  the  diagram.  It  will  be  seen  that  the  control  valve  has  3  positions.  In 
the  position  illustyato-3  pressvie  is  applied  to  the  underside  of  the  power  pi.scon 
and  the  fluid  from  above  the  power  pis:  on  flows  through  the  valve,  and  back  to  the 
reservoir.  When  -he  valve  position  xS  changed,  as  shown  by  the  dotted  lines,  fluid 
is  forced  by  the  pump  into  the  actuator  a'-ove  the  piston  -end  fluid  below  the  piston 
passes  through  the  valve  and  back  to  the  reservoir.  In  tfcxt  way  it  is  possible  to 
obtain  a  power  stroke  in  either  direction  using  the  combination  of  4-way  valve  and 
double  acting  actuator. 


3 . 3  The  U.°e  of  Matrices  for  Kc Celling  Part i  ouJ ate  Systems 

3.3.1  In  order  to  construct  a  mathematical  model  for  a  system  in  which  individual 
components  may  alt^r  the  concentration  of  particles  it  is  necessary  to  be  able  to 
describe  efficiently  the  particle  si2e  distribution  and  toe  way  in  which  this  is 
changed  by  the  system.  One  convenient  technique  is  <■  vrp^esent  particle  3iz6 
distributions  as  column  vectors,  and  operatois  »7hic>  change  ti.em  by  square  matrices. 
This  has  the  advantage  tha any  calculation  requires  only  the  u?«  of  matrix  algebra. 

3.3.2  The  simplest  method  of  representing  a  particle  size  distribution  as  a  vector 
is  to  equate  the  elements  of  the  »*yo tor  with  the  olumns  of  a  frequency  'size 
histogram  for  the  distribution  Thus  the  histogram: 


%  by  weight 
of  total 
sample 


X2 


Size  x 


would  be  represented  by  the  vector: 


S1 

= 

10 

S2 

20 

S3 

= 

4C 

S4 

= 

30 

where  is  the  percentage  by  weight  in  the  size  range  x^  to  is  the  percentage 

by  weight  in  tne  size  range  x2  to  x, ,  and  so  on.  Using  this  technique  it  is 


possible  to  represent  the  particle  size  distribution  In  terms  of  any  convenient 
parameter,  e.g.,  weight,  number,  volume  or  any  of  the  other  parameters  used  in 
particle  size  analysis. 

3.3.3  An  alternative  method  is  to  let  s  ,  s2,  etc.,  be  actual  weights  in  each 
siz6  fraction.  The  sum  of  the  elements  of  the  vector  is  then  equal  to  the  total 
weight  of  particles  in  the  system.  Thus  for  example,  a  particle  sample  with 
the  following  sieve  analysis: - 


TABLE  III 
Sieve  Analysis 


Size  range 

Weight  retained  gw. 

300  pm 

0.03 

300-200  pm 

0.28 

200-100  pm 

0.47 

100  pro 

0.41 

Tefal 

1.19  j 

1 
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could  be  represented  in  vector  form  by:~ 


3.3.4.  The  use  of  actual  weights  rather  than  percentages  by  weight  has  one 
important  advantage  wh^n  matrix  '‘■echniques  are  used  to  '^present  particulate  systems. 
For  example,  the  perfect  mixing  of  two  particle  size  distributions  r  and  s  is 
represented  by  addition  01  the  veci -jrs ,  provided  that  both  size  distribution  vectors 
are  in  terms  of  actual  we-'-hts, 


Size  distribution  vector 
for  mixed  material. 


r3 

r4 

»»  m 


If,  however,  the  sise  distribution  vectors  are  expressed  in  terms  of  percentage  by 
weight,  then  they  must  be  multiplied  by  the  sample  weight  p  and  p  before  addition 
can  be  psrformed. 


i,e., 

P  [r. 


<r\  S. 


(  C  +  6*)  1  P  r, 


+ 


p  + 


+  6"s„ 


r3  S3  jl  - p^r - 

r-  I  5+  0r3  +  gsc 

.  1  P  +  r 

8" ‘4  *  g»4 

[  0  +  f 

^  « 

For  this  reason  it  is  itor>  useful  to  express  particle  size  distribution  vectors  in 
terms  of  weights  ra+her  than  o-;vcentagos  by  weight,  and  the  former  method  of 
representation  is  <>ned  in  the  following  uhap^ers  of  this  report. 
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Where  it  is  more  appropriate,  weight  concentrations  in  the  fluid,  or  weight  feed- 
rates,  i.e.,  weight  passed  per  unit  time,  are  used  instead  of  actual  weights. 

3.3.5.  Operations,  other  than  perfect  mixing,  which  affect  particle  size  distri¬ 
butions  may  be  represented  by  square  matrices.  Matrix  models  for  classifying  and 
steady  stale  grinding  operations  have  already  been  developed  to  a  considerable 
extent 3  ref. 3. 1  j and  matrix  modelling  work  in  other  types  of  particle  processing 
should  prove  equally  fruitful. 

In  general,  if  a  size  distribution  r_  is  changed  to  a  size  distribution  j,  by  a 
process,  then  the  process  may  be  represented  by  a  matrix  A  such  that: 


s  =  A  r 


multiplying  s  +  A  r  gives: 


<r  — 

'1 

aiiri  + 

3l2r2 

a  r 

13  3 

-*■  a  r 
14  4 

S2  = 

a2irl  + 

a22r2  + 

a23r3 

+  a24r4 

S3  = 

a3iri 

a32r2 

a33r3 

a34r4 

3  4  = 

a4iri 

a42r2 

a4Cr3 

a44r4 

The  coefficients  ajj  of  A  represent  the  traction  of  particles  of  size  j  which  are 
changed  to  size  i. 

3.3.6-  In  many  cases,  several  of  the  coefficients  of  A  are  Zero.  For  example  a 
classifier  is  usually  considered  not  to  cause  changes  in  particle  sizes.  The  only 
non  zero  elements  a^j  are  thus  those  where  i  =  j.  A  typical  classifier  matrix  C, 
has  the  form: 
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For  a  perfect  classifier  the  coefficients  c  are  either  unity  or  zero,  depending 
upon  whether  a  size  range  is  greater  or  .smaller  than  the  classifier  size.  For  a 
real  classifer,  however,  the  coefficients  represent  t*:e  classifier  efficiencies  for 
the  size  ranges  concerned. 

2.3.7.  In  a  grinding  machine  it  is  reasonablo  to  assume  that  no  size  enlarge¬ 
ment  takes  place,  therefore  a  matrix  representing  grinding  would  contain  zero 
coefficients  for  all  those  terms  representing  enlargement.  The  final  matrix  would 
be  lower  triangular  in  form: 

f 

\i 


"'ll 

0 

0 

0 

m21 

ffl22 

0 

0 

m31 

“32 

m33 

0 

*41 

"42 

m43 

m 

44 

Since  there  is  no  generation  or  loss  of  solids  in  the  mill  then  the  sum  of  all  m 

for  any  constant  .i  is  unity,  i.e.,  the  sum  of  all  fractions  of  particles  from  a 
given  size  range  is  equal  to  the  weight  of  particles  in  that  size  range. 


Thus: 

rl 

52 

m  r 

U  1 

+ 

m2iri  +  ra3iri 

+  m4irX 

r2 

= 

m22r2  -  ™32 

m42r2 

r3 

= 

m33r3 

+  m  r„ 
43  3 

r . 

j. 

nr  ,r 
41  4 

This  property,  of  the  sums  of  the  elements  in  each  column  being  ecual  to  unity,  is 
comaon  to  all  operator  matrices  which  change  only  the  particle  size  distribution, 
and  not  the  total  weight  of  the  particles. 

3.3.8.  A  size  enlargement  device  may  also  be  represented  by  a  matrix.  Ir  this 


case  no 

particle 

is 

decreased 

i  *'  si  ze , 

sc  the 

matrix 

matrix: 

K  ~ 

F*n 

& 

12 

ei3 

ei4 

0 

&22 

e 

23 

e24 

0 

0 

e33 

*34 

°44 


0 


0 


0 


3.3.9.  All  operations  involving  particles  may  bo  considered  to  be  a  combination 
of  the  three  basic  operations  above,  classification,  size  reduction  and  size 
enlargement.  A  matrix  representing  any  operation  can  thus  be  built  up  from  the 
three  basic  processes  by  analysis  of  the  operation  in  terms  of  these  three  basic 
processes.  For  example,  filtration  is  a  form  of  classification;  erosion  is  a  form 
of  grinding  although  the  generation  of  particles  results  in  the  sums  of  the  matrix 
columns  being  greater  than  unity, 

3.3.10.  The  values  of  the  matrix  coefficients  are  bes.  determined  experimentally, 
e.g. ,  the  coefficients  of  a  classifier  matrix  are  the  classifying  efficiencies  for 
each  size  range.  Much  work  is  being  done,  however,  towards  theoretical  determin¬ 
ation  of  matrix  coefficients,  especially  in  the  matrices  associated  with  grinding 
operations. 
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Characterisation  of  Individual  Cof-po cents  in  c  ,ivdu 'ul.  e  Syst ten 

J.4,1  in  order  to  set  up  a  mathematical  modei  foi  the  complete  system  it  is 
necessary  to  find  a  matrix  description  of  each  component  in  the  circuit.  The  main 
items  considered  are  the  pump,  filtei,  pressure  reliif  valve,  reservoir,  the  4-way 
valve,  and  the  actuator.  The  volumes  of  fluid  retained  in  the  connecting  pipework 
must  also  he  considered.  Details  of  models  used  for  each  item  in  the  circuit  to¬ 
gether  witn  tlie  ssumptions  made  in  order  to  devise  each  model  are  given  in  para¬ 
graphs  3.4.3  to  0.4,8  below.  Paragraph  3.4.2  summarises  the  nomenclature  used  in 
this  and  sufcysqcent  chapters,  Fig.  16  shows  the  particle  concentrations  and 
fig..  IT  the  volume  elements  m  the  model  system. 

3.4.2  Consistant  units  must  be  employed  for  the  particle  size  distribution  vectors 
at  all  points  ;n  the  system.  Unless  otherwise  stated,  all  particle  size  distribu¬ 
tion  vectors,  in.  this  and  s  tcceedmg  chapters,  have  the  units  (weight  of  particles) 

(unit  volume  of  fluid).  The  coefficients  of  all  me  trices  are  dimensionless.  The 
notation  used  throughout  is: 

(1'*  Upper  case  letters  -  all  matrices 

(V)  Lower  case  letters,  underlined  -  all  particle  sice  distribution  v^ctois. 

(3)  Lower  case  letters,  not  underlined  and  groek  smhois  -  scalars,  usually 

t cioci ties  cr  volume*:. 

r-  g.  (a  -  -. )  \.H.  £ _ 

Scalars  Matrices  Vectors 

A  list  of  all  the  symbols  used,  'n  alphabetical  order,  is  given  on  pace4  xiu- 
x\  i 

3.4.3  The  Fump 

The  pump  considered  in  t  le  circuit  is  a  constant  displacement  swr.sh  plate 
piston  pump  of  the  type  used  ii  aircraft  systems.  The  major  points  oi'  wear  and 
hence  sources  of  particles  will  be  considered  as  the  valve  plate  and  the  piston 

shoes,  "here  is  comparative)-  little  erosion  of  the  pistons  and  cylinder  barrel 

and  this  may  be  ignored  as  a  first  approximation.  Also  the  simplifying  assumption 
is  made  that  the  rate  of  wear  of  the  valve  plate  and  piston  shoes  is  independent 
of  the  concentration  of  particles  in  the  fluid,  but  depends  upon  the  speed  of 
rotation  of  the  pump  and  to  some  extent  upon  the  age  Of  the  pump.  It  is  also 

assumed  that  there  is  no  leakage  of  particles  in  the  case  drain  fluid 
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(1)  The  eroded  particles  may  because  of  the  constructional  materials  of  the 

pump  ,  be  very  hard  and  nught  be  expected  to  cause  further  erosion  of  the  pipework 
and  other  items  m  the  circuit.  Assuming  again  that  the  pump  is  operated  at  con¬ 
stant  speed,  the  rate  of  erosion  will  he  constant,  the  fluid  flow  rate  will  he 
constant  and  hence  the  concentration  of  particles  m  the  fluid,  as  the  result  of 
pump  erosion,  will  he  constant.  New  m  is  the  particle  concentration  m  the  fluid 
leaving  the  reservoir  and  b^  is  the  increase  m  concentration  as  a  result  of 
breakdounof  the  line  between  the  reservoir  and  the  pump. 

*  the  concentration  at  the  pump  inlet  =  m  +  Ojj,. 

Let  the  concentration  at  the  pump  cutlet  be  p. 

Let  the  increase  m  concentration  as  a  result  of  erosion  of  the  valve  plate  be  a. 

If  it  is  assumed  that; 

(1)  there  is  no  accumulation  01  part’eies  inside  the  pump, 

(2)  there  is  no  loss  of  particles  as  a  result  cf  leakage  into  the  case, 
then  a  particle  mass  balance  over  the  pump  gives: 

p  “HI 

Let  the  particle  concentration  of  eroded  material  from  the  piston  shoes  be  d.  If 
it  i*  assumed  that  there  is  no  accumulation  of  particles  inside  the  case,  men  the 
particle  concentration  of  the  fluid  leaving  the  case  drain  valve  is  d. 

3,4.4.  The  Kilter 

The  filter  acts  as  a  classifier  on  the  particles.  The  particle  concen¬ 
tration  at  the  filter  inlet,  including  particles  in  the  iluid  as  a  result  of  break¬ 
down  of  the  pipework  is  <p  +  bp) . 

Let  the  output  concentration-  be  f. 

Then  f 

where  F  is  a  diagonal  Matrix  of  the  form: 


+  %  +  £• 


F'  \P  f  £p> 


=  ft 


11 


0 


f 


■22 


0 

0 


0 

0 


fS3  0 
0 


l44 


The  matrix  coefficients  may  be  represented  by  the  filter  efficiencies  for  each 
sixe  range.  These,  however,  are  not  constant  but  depend  upon  the  fluid  velocity 
and  acceleration,  particle  loading, and  cake  thickness. 


■?0, 


3.4.5  The  Pressure  Relief  Valve 

In  the  pressure  relief  valve,  erosion  of  the  valve  seat  by  the  hard  part- 
icIgs  occurs  fxrst,  followed  bv  collection  of  particles  inside  the  valve.  Since 
many  of  the  particles  are  nard,  and  there  are  considerable  difficulties  in  the 
practical  differentiation  between  hard  and  soft  particles,  we  have  acsiimed,  for  the 
purposes  of  the  model,  that  all  the  particles  are  hard.  It  nay  even  be  possible  to 
incorporate  this  difference  m  the  coefficients  of  the  erosion  matrix,  although  this 
is  beyond  the  scope  of  this  report. 

Tne  particxe  concentration  at  the  pressure  relief  valve  inlet  including  pipewont 

breakdown  is  (_f  +  +  b^)  .  After  erosion  this  is  changed  to  Epv  (_f  +  ^  +  ]0 

where  C  is  a  lower  triangular  matrix  similar  to  that  used  to  represent  a  grrnd- 
PV 

mg  machine.  E.g.:  F 

"PV 


r*u 

0 

0 

C  ' 

e 

21 

S22 

0 

0 

°31 

e32 

e33 

0 

e 

41 

°42 

e43 

°44 

This  concentration  is  changed,  by  the  collection  of  particles  inside  the  valve,  to 

C  E  (f  f  b  +  b  )  where  C  is  a  diagonal  matrix. 

pv  pv  —  — f  — g  pv 

This  concentration  is  equal  to  the  outlet  concentration  z 


3.4.6 


z  =  z  (£ 

•  * 

The  Reservoir 


-f 


+  b  ). 

s 


The  reservoir  is  assumed  to  be  a  perfect  mixer.  Since  it  is  assumed  that 
there  >s  no  leakage,  and  there  is  no  generation  or  settling  of  particles  in  the 
reservoir,  then  the  reservoir  composition  is  easily  calculated  by  a  mass  balance 
across  the  reservoir.  The  method  used  to  determine  this  mass  balance  depends  entire¬ 
ly  upon  the  model  used  for  the  circuit.  The  methods  used  for  our  final  model  are 
described  in  detail  in  3.7. 

fhe  rese-vcir  outlet  contains  a  filter.  If  the  particle  concentration  inside  the 
reservoi-  is  r,  and  the  outlet  concentration  is  m,  'hen:  m  =  R£ 
where  R  is  a  diagonal  matrix  rep-esenting  the  filter,  e.g.: 


R  = 


1  j 


22 


0 

0 

r33 

0 


44 


41. 


3.4.7 


The  4 -way  Valve 


The  directional  control  valve  is  a  4-way  servo  valve,  but  in  order  to 
simplify  the  model,  is  only  considered  as  having  3  operating  conditions,  the  neutral 
position  and  the  hard  over  position  m  either  direction.  As  an  example  of  r  valve 
for  this  circuit  or.e  could  consider  a  manually  operated  valve  similar  to  the  second 
stage  of  an  electro-hydraulic  spool  and  sleeve  valve.  This  valve  is  illustrated  in 
figs.  13-15. 

It  will  be  assumed  that  particle  build  up,  i.e. silting  occurs  in  those  positions 
within  the  valve  where  high  pressure  fluid  is  leaking  to  the  return,  and  these-  pos¬ 
itions  will  be  called  Xj,  Xg  and  Xg.  The  close  fit  between  the  valve  spool  and 
spool  cylinder  will  cause  the  leakage  points  to  act  as  filters,  collecting  particles 
at  points  x^,  and  x ^  and  X3.  The  particles  remaining  m  the  leaking  fluid  then 
cause  erosion  of  the  spool  and  cylinder  wall. 

The  inlet  particle  concentration  for  the  valve,  including  breakdown  of  pipework  is 

<£  +  bf  +  bjj) 

Let  the  outlet  concentration  of  the  valve  be  £ 

*■  =  E.C  (f  +  b,  +  bh  ) 

.  .  —  vc  vc  —  — f  — 

where  Evc  is  a  triangular  matrix  representing  erosion,  and  Cvc  is  a  diagonal 
matrix  representing  particle  collection.  Eyc  and  Cvc  can  be  multiplied  to  give  a 
triangular  matrix  Lc>  representing  the  overall  leakage  effect  when  the  valve  is 
closed.  This  gives: 

L  =  Lc  <1  +  «>f  +  bj,  > 

Since  the  valve  is  of  symmetrical  design  it  is  reasonable  to  assume  that  the 
action  of  the  valve  on  the  particles  is  exactly  the  same  when  the  actuator  is 
being  raised  as  it  is  when  the  actuator  is  being  lowered.  When  the  actuator  is 
being  raised,  (fig.  15  ),  the  only  action  of  the  valve  itself  is  erosion  by  the 

particles  in  the  fluid.  This  is  quite  small  and  has  been  ignored.  However,  the 
concentration  of  particles  leaving  the  value  is  affected  by  the  particles  collected 
ai  Xg,  x^j  and  x-  during  leakage. 

Total  weight  of  particles  raptured  (T.-Cvc)  (£  +  bj  +  b^) 

where  C  is  the  total  volume  of  fluid  which  has  leaked. 

However,  only  a  fraction,  <J>  ,  of  this  weight  collects  at  point  X3. 

Further,  accumulation  of  particles  at  point  Xg  will  probably  result  in  compacting 
of  the  particles.  Thus  when  the  valve  is  open  some  particles  will  probably  adhere 
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FIG  15  VALVE  ACTUATOR  MOV  ING  TO 


to  ch'\  orifice  walls.  If  tee  fraction  of  partioca  which  does  not  adhere  is  9, 
than  the  weight  of  particles  which  enters  the  f  lwid  being  pumped  into  the  actuator 
■fs  given  by? 

Weight  of  particles  released  torn  X3  *  0d  l£.  (I-C^)  (f  +  +  b^). 

Whilst  the  actuator  is  being  raised,  furthr.r  leakage  occurs  at  points  X4  and 

Xy  The  iSecitanisE  of  this  leakage  is  exactly  the  same  as  for  leakage  when  the 

valve  is  closed,  tne  only  difference  being  that  the  leakage  points,  and  heDce  the  matrix 

coefficients,  are  different.  The  equation  for  the  leakage  effect  whei.  the  -alve  is 
open  is  thus: 

A  *  L0  +  kf  * 

where  LQ  **  Evr  ♦  Cyj,  (the  matrices  being  similar  in  structure  to  ^ c *  ^vc  and  ^ve 
above) . 

The  total  weight  of  particles  captured  as  a  result  of  this  leakage  is  thus: 


£  (I~C  )  (f  +  b+  b^) 

where  o  is  the  total  volume  of  fluid  winch  her  leaked. 

Since  the  high  pressure  fluid  is  in  motion  it  is  reasonable  to  assume 
that  ail  particles  captured  as  points  and  X5  are  carried  into  the  actuator  by 
the  motion  of  the  fluid.  Thus  if  the  fraction  of  the  total  particles  captured 
which  arc  captured  at  points  and  X5  is  ?  then  the  weight  of  particles  bsing  pump¬ 
ed  to  the  actuator  is: 

o 

*  ^  (I-Cvo)  (-f-  %  +  -V 

If  the  concentration  of  particles  in  the  fluid  as  it  leaves  the  valve  before 
entering  the  actuator  is  v,  then:  c 

1  -  (i  +  b£  +  bh>  +  6r  21  a-Cvc)  a  +  hf  +  V 

*-  1 

**  |  (I-C  )  (f.  +  K  +  b  ) 
c.  vo  —  — f  — h 

uhere  is  the  total  volume  of  fluid  pumped  into  the  actuator. 

The  effect  of  the  valve  on  low  pressure  fluid  returning  from  the  actuator 
is  erosion.  Since  this  is  only  slight  it  may  be  ignored.  Thus,  if  the 
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concentration  at  the  valve  inlet,  including  pipework  breakdown,  is  (j_  +  A> , 

then  the  outlet  concentration  is  given  by: 

t  -  (I  +  b,  ,). 

—  "kA 

The  valve  action  when  the  actuator  is  lowered  is  similar  to  that  above,  the 


equations  being: 


=  (f  r  b  +  b,  )  + 
-  —  f  ~h 


H  (i-r.  )  (f  +  b  +  b  ) 

~  vc  f  ~h 


+¥  -s:  (i-c  )  (f  +b  +  b  ) 

1  vo  f  > h 

-  '  (1J  *  V 

The  values  of  the  total  leakage  volumes,  c  and  o,  are  not  necessarily 
the  same  as  for  these  whan  the  jack  is  being  raised. 

3.4.8  The  Actuator 

The  action  of  the  actuator  is  basically  particle  capture,  at  the 
inlet  nozzle,  on  the  piston  rings  and  on  the  piston  gland  (side  B  ^nlvl .  On  side 
15  there  is  an  additional  effect  due  to  ingress  of  fresh  particles  from  the  surface 
of  tnl  piston  rod. 

Consider  side  3,  efter  it  has  been  filled  n  times. 

Total  weight  of  solids  in  actuator  =  (contents  ci  Ime  5)  +  '"maLe-up  from  line  3) 

+  (ingress  f rc j  piston  rod). 

The  volume  of  line  5  is  V,  „  and  t'  u particle  concentratiou  is  i_  -1. 

kd  Bn 

If  Lne  collection  effect  as  fluid  enters  the  actuator  is  Ujgj,  then  the  total 
weight  of  particles  which  enters  the  actuator  from  1 ine  5  including  pipework 

breakdown  is* 

V.:3  CjBi  (^n-l  +  \.-,} 

The  voiume  of  1  iuid  from  line  3  which  is  ^equi^ed  to  make  up  the  actuator  volume 

is  (V  -V  ),  assuming  that  the  actuator  volume  is  larger  than  the  volume  of  line 
iB  kB 

5.  fl.c  *  oncent ration  of  fluid  entering  the  iccuator  from  line  3  is  Cv  +  bj~)  •  Thus 
the  total  weight  c  particles  from  line  3-.  including  pipework  breakdown,  which 
enters  the  actuator  is: 


C. 


Bi 


2 
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The  reason  fcr  a  sum  's  that  v  is  not  constant  but  varies  with  eac-h  stop  change. 


The  weight  of  particles  which  enters  the  actuator  as  external  contamination 
on  the  piston  rod  is  assumed  to  be  constant  and  equal  to  k' 

Let  the  total  weight  of  solids  inside  the  actuator  be  s 

BN 


.  .  -BN 


VkB  CjBl  +  +  CjBl  jE  £  kB  (v  +  b.  _)  +  k*. 


.  .  The  particle  concentration  inside  che  actuator  is  1 


As  side  B  empties,  particles  are  collected  on  the  piston  rings,  outlet  nozzle,  and 

piston  rod  gland.  If  these  effects  arv-  represented  by  the  matrices  C  c 

J  I  J 

and  C  respectively,  then  the  actuator  outlet  composition  j  is  given  by: 
jG  -t  ~Bn 

j  -  j  (c  +  C  „  +  C  )  S  j 

%n  v  3  jB2  jB3  jG  -Bn 


-ft 

C  .  ) 

s  j 

JG 

-Bnj 

(v  -v 

+  b  \ 
— k»  j 

h  '  %' 

/ 

OB  1  "  "  J 

.  .  4n  =  v'~  j*  jB  TVkC  ^is.,-1  +  --k!^  "  ^  £  (v  +  b^Jl+K.k 

«h6«  CJP  -  (C.B2  *  ctB,  *  c.c)  c,B1 

K  ■  Km  "  CJB3  *  V 

If  the  mc-ease  in  particle  concentration  in  the  actuator  as  r  -esult  of  ovtemal 
/ 

contamination  is  !<,  then: 

t 

k*  v  .  k. 

jB  - 


Rearranging  the  equation  for  j.  g’ves: 

_  — Bn 


i  fc  r*  • 

-  jB  \  KB 

vjo  *■ 


J-Bn  -1 


(v  -v 
-ir-£kB 


t  c  ,  b 
IB  -kB 


by  a  similrr  process  the  co.icentrut ion  of  particles  leaving  side  A  of  the  actuator 
after  it  has  bojn  filleo  n  times  \s  given  by: 


i  =  -—  fc  f  j  +  (VjA^VkA>  1]  + 

i\.n  VjA  ^  jA  kA  -An-1  vjj 


C  b 
jA  -kA 


where 


CJA  =  <>2  *  V’  Vl 
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Note:  Since  there  is  no  piston  rod  on  side  A,  k  and  C  are  both  zero. 

—  JG 


3.4.9  Erosion  of  Pipework 

The  erosion  of  pipework  has  been  represented  in  the  sections  above  as 
addition  of  particles  rather  than  t-j  ,r.a trices .  Although  pipework  erosion  is  to  some 
extent  dependent  upon  particle  concentration  it  is  both  a  very  complex  phenomenon 
and  also  occurs  on  a  very  small  scale,  hence  it  will  be  ignored  from  no’.t  on. 
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3.5. 


McJes  of  failur o  Caused  by  Particle  Conta^ira t ion 
1  Sou  rocs  ot  Fluid  Contaoination 

Cl)  Built-in  contamination,  i.e.  dir*  that  ir  riready  in  the  components 

when  the  system  is  assembled.  When  hydraulic  components  are  manufact¬ 
ured  they  undergo  a  rigcrouo  ult^aso^ic  cleaning  process  followed  by 
flushing  with  clean  fluids.  However  it  has  been  tourvi  *hat  when 
assembling  the  most  critical  systems  that  the  act  cf  tightening  a 
screw  thread  will  release  a  considerable  number  of  particles  into  the 
system. 

(Reference  3.5.1)  These  systems  are  now  being  designed  sc  tha;  no 
screw  threads  are  m  contact  with  the  fluid  circuit. 

(2)  Contaminant  m  ne.f  hydraulic  tluid  added  to  the  system. 

(3)  Airborne  dust  infiltrating  the  system,  e.g.  dust  brought  in  on  the 
piston  shaft  and  iLiough  reservoir  breather  valves  if  these  exrst. 

(4)  Contaminant  generated  within  the  system  by  pu.un  wear;  packing  wear, etc. 

(5)  Contamination  from  hand  tools  etc,,  introduced  when  a  system  is  opened 
for  maintenance  reasons. 

(6)  Contamination  introduced  by  intercha-.ge  of  fluid  with  that  of  ground 
maintenance  equipment. 

2,  4,  5  and  6  may  he  largely  controlled  by  effective  maintenance  of  ground 
equipment  and  3  by  efficient  seal  design. 

■ 2  Types  o'  failure 

Contarei nation  induced  failures  of  hydraulic  systems  and  their  components 
bo  divided  into  thro-  cetngo.iss .  {l}efJ3.5.2> 

Malfunction  dus  to  single  large  particles,  or  a  few  large  particles 
arriving  s i.vrl  tansously ,  mterfer * 05  with  the  motion  of  moving  ports,  e.g. 
reseat  tag  of  valves. 

Loss  of  pcrforffiAnoe  and  eventual  fnlure  due  to  woar. 

Loss  of  ponoroonoe  when  frictional  forces  become  significant  with  respect 
tc  drivir--  fui-css,  or  v^en  heavy  contamination  affects  the  chiving  forces, 
e-s  s.Itir<£  in  sc*  vc  control  valves. 

3.5.3  A  particle  large  enough  to  plug  vn  orifice,  prevent  :  valve  from  seating, 

>-£e.,  is  laige  enough  to  be  easi'y  teen  and  you  Id  be  trapped  by  the  coarsest 
filter.  I  to  presence  is  indictive  of  carelessness  on  the  pa.-t  Of 


3.5 

may 

(1) 


(2, 

(3) 
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manufacturing,  overhaul,  or  maintenance  personnel.  Failure  to  debur*-  a  drilled 
passage  at  the  time  of  manufacture,  for  example,  may  eventually  result  in  the  burr 
being  dislcsged  into  the  system.  Good  maintenance  and  the  installation  of  coarse 
filters  at  critical  points  will  prevent  this  type  of  failure. 

Another  failuie  of  this  typo  is  possible.  This  is  the  proba.lility  of  the 
simultaneous  arrival  of  a  number  of  particles-  at  an  orifice,  etc.,  whose  combined 
effect  is  the  same  as-  that  described  above  for  a  large  particle.  CRef.  3.5.3) 

3.5.^  Failure  due  to  a  large  partic  .e  bloching  an  orifice  will  be  dependent  upon 
the  sizes  of  the  orifice  and  the  particle.  If  tno  particle  is  above  a  certain  size 
it  will  cause  block  ~e,  hue  this  is  very  unliHeiy. 

Failure  due  to  a  number  of  particles  causing  blockage  can  be  calculated  on 
a  probability  basis,  the  important  criteria  being  the  orifice  size  and  the  particle 
size  distribution,  i.e.  two  particles  of  one  size  range  could  cause  blockage,  whilst 
four  of  a  small  size  range  might  be  needed,  or  a  combination,  of  particles  of 
different  size  ranges. 

3.5.5  All  moving  parts  are  subject  to  wear  at  points  of  contact,  ft.e  priiiz  s  pci  *. 
exceptions  are  hydraulic  pumps,  and  continuous  duty  hydraulic  Rotoi  s .  Infiltration 
of  the  precision  fits  on  e^uipmer.  ■...  fluid  borne  particles  results  in  veer,  scoring 
and  increased  clearances.  This  sometimes  causes  s  aider:  ftiluro  but  more  often 
results  m  a  gradual  loss  of  efficiency:  >•!  jr!-'*d  oia’  venar.ee  cap.  pr-svent  this  be¬ 
coming  critical. 

3.5.6  Failure  of  lt^ms  due  io  wear  car.  be  calculated  on  a  total  weight  of  solids 
lost  basis,  o . go  *f  tf:  n.ur.p  loses  .tore  ihar  a  certain  wei  bt  of  particles,  it  will 
fail . 

S.n,'  Control  valves,  'n  particular  elect rohynaul ic  servo  valves,  may  nave 
limited  fcrces  available  to  drive  the  valve  speyi.  Frictional  fc~cee  may  tr-comv 
great  enough  to  ozu-'-e  degradation  of  vei"f 0 ■ arje  r.J  thcugrs  little  nr  nr  wear  has 
occur. td . 

jVrvo  valves,  «lso  called  t runnier  mV  ■:  or  hyurAui  ;c  r:nril>  f ie^s  ;  are 
usually  designed  to  proobc..  hydraulic  t !  ow  output  proportio  to  an  electric 
current.  It  is  the  performance  of  these  units,  more  than  any  o-.hc.-'  lector,  -vhich 
has  focussed  attention  on  hydraulic  fluid  contair.irai.on  in  che  last  few  years. 


The  valve-  design  U3©d  m  this  project  is;  nimiiar  to  tiie  second  stage  of  a 
typical  oleetrohydraulic  servo  valve. 


In  a  spool  v«l; 


i'-g:.  13  ,  particles  tend  to  colloct  behind  the  spool 


at  pc-nto  x  ,  x  art*  x  due  to  leakage  effects..  These  particles  increase  resistance 

•*>  *•  J 

tc  motion  of  the  spool,  and  wav  delay  or  even  5top  the  notion,  cuus.  ng  eithei  complete 
failure,  or  a  loss  in  efficiency  Chich  may  be  critical  is'  the  valve  is  part  of  a 
flight  control  loop. 
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3.3.  Models  of  Complete  Systems 

3.6.3.  The  requirements  of  a  model  of  an  hydraulic  system  designed  to  predict 
failure  of  components  due  to  particle  action  arC-:- 

(1)  it  must  be  able  to  calculate  continuously  changing  particle 
concentrations , 

(2)  it  must  be  able  to  handle  continuously  operating  and  intermittently 
operating  circuits, 

(3)  it  must  be  able  to  predict  failure, 

(4)  it  must  be  able  to  handle  changes  in  matrix  coeflxcients.  where 
they  aye  dependent  upon  particle  concentrations,  or  wite~e  they  are 
changed  by  routine  replacement  of  equipment  represented  by  the 
matrices , 

£6}  it  c.ust  be  capable  of  modiiicatlon  to  simulate  other  hydraulic 
systems , 

(6)  it  must  be  in  form  suitable  for  solution  bv  comouter. 

3.6.2.  There  aie  two  distinct  types  of  model  which  fulfil  these  requirements: 

(1)  Continuous, 

r2)  Step  wide  based  on  either: 

(1)  volume  increments  of  fluid,  or 

(2)  time  increments. 

3.6.3.  /,  continuous  model,  which  calculates  particle  concentrations  as  continuous 
functions  of  time,  or  possibly  voxurne  of  fluid  pumped  theory.  However,  matrix 
models  are  ur.suited  to  continuous  solution,  parti  •ulariy  when  the  mc.trix  coeffic¬ 
ients  themselves  --ary  with  time. 

3.6.4.  The  basis  of  the  step-wise  approach  is  to  assume  that  a  variable  function 
is  ~onstant  for  a  small  increment  of  a  fixed  parameter,  and  then  changes  ior  the 
next  increment. 

The  simplest  parameters  to  consider  a*-e  increments  of  volume  and  time.  The  values 
of  the  particle  concentrations  in  any  increment  depend  upon  the  values  in  the 
previous  increment  and  upon  the  sise  of  the  increment.  Each  time  these  values  are 
calculated,  since  they  depend  upon  particle  concentrations. 
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Il'  volume  increments  are  used,  and  the  particle  concentrations  are 
multiplied  by  the  fluid  velocity  to  give  particle  fluxes,  the  model  will  handle 
acceleration.  This  is  discussed  below  as  the  step-wise  flex  model. 

3.6.5.  With  the  step-wise  volume  model-  hize  distributions  can  be  very  easily 
represented,  all  the  requirements  of  3.6.  ?  are  met.  however  ac.c2_err.fion  of  the 
fluid  cannot  easilj*  be  dealt  with  by  this  nodel. 

A  step-wise  time  model  again  satisfies  all  the  initial  requirements,  cut 
once  again  acceleration  cannot  be  handled.  Particle  sire  distributions  must  also  be 
expressed  in  terms  of  volume  and  velocity  vhicn  mokes  this  model  .less  attractive  than 
a  volume  model. 

A  step-wise  velocity  or  flux  -r.odel  is  the  cnly  model  which  can  h<  sidle 
accelerating  fluid.  It  is  also  the  best  model  for  represen; ing  the  continuous 
circuit.  Di ificuj ties ,  however,  arise  when  the  intermittent  circuit  (the  actuator) 
is  brought  into  action. 

3.6.6.  For  the  first  model  it  wr.s;  decided  to  represent  the  continuous  circuit  by  a 
flux  model,  and  the  intermittent  circuit  by  c  volume  model.  The  development  of  this 
model  into  the  final  model  adopted  is  described  below  in  3.7. 
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3.7 


The  Filial  Mathemati cal  Model  of  a  Simple  Hydraulic  System 


3.7.1.  The  first  model  sras  based  on  a  division  of  the  basic  hydraulic  system  shown 
in  Fig.  lc  ,  into  two  distinct  pa:ts,  a  continuously  operating  circuit  and  an 
intermittent  circuit.  The  continuously  operating  circuit  consists  of  flow  from  the 
reservoir,  through  the  pump,  the  filter  and  the  pressure  relief  valve,  and  back  to 
the  reservoir,  except  when  leakage  occurs.  During  leakage,  flow  is  from  the 
reservoir,  through  the  pump,  the  filter,  leakage  in  the  Control  Valve,  and  back  to 
the  reservoir.  The  intermittent  circuit  comes  into  operation  when  the  actuator  is 
moved.  It  has  been  previously  described  in  paragraph  3.4.8. 

The  continuous  circuit  was  modelled  by  means  of  a  mass  balance  on  the 
reservoii ,  cased  upon  fluid  velocities  and  particle  concentrations.  This  enabled 
the  particle  flux  leaving  the  reservoir  to  be  calculated.  The  main  advantage  of 
this  moael  war  that  it  could  handle  simultaneous  flow  through  the  pressure  relief 
valve  and  the  Control  Valve,  (Models  based  upon  volume  or  time  can  only  handle  this 
ov  alternating  the  flow  from  one  circuit  to  the  other  for  successive  volume  or  time 
increments  I . 

In  order  to  calculate  particle  concentrations,  it  was  assumed  that  the 
concentration  leaving  the  reservoir  remained  constant  until  particles  travelling  at 
the  mean  fluid  velocity  had  made  one  complete  circuit.  Matrices  were  then  adjusted 
and  the  process  was  repeated. 

The  intermittent  circuit  was  based  on  an  actuator  volume,  assuming  that 
there  was  only  one  step  change  in  particle  concentration  for  each  actuator  stroke. 

The  model  for  the  continuous  circuit  is  the  best  available  type  since  it 
can  easily  handle  fluid  acceleration.  However,  it  was  found  to  be  extremely  diff¬ 
icult  to  combine  both  circuits  into  a  single  model  without  drastically  modifying  one 
or  ether  of  the  models. 

3.7.2.  A  second  model  was  then  developed  from  the  first  by  substituting  a  step¬ 
wise  volume  for  the  flux  model.  Again  it  was  assumed  that  only  one  change  in  con¬ 
centration  was  made  per  circuit. 

Acceleration  was  handled  by  considering  it  to  be  infinite.  This  resulted 
in  step  changes  in  velocity  which  were  accounted  for  hy  step  changes  in  the  ope.-ating 


FIGURE  18  Model  hydraulic  system 


This  model  inherited  most  01  the*  disadvantages  of  the  first  model,  the 
:wo  most  important  being: 

*jl)  Difficulties  arose  in  combining  the  modexS  fo.  the  continuous  and 
intermittent  circuits. 

(2)  .Although  the  vc.’.ume  mcrowent  was  small,  tl'e  st.e>.  change  in 

narticlo  concentration  cas  quite  large-  since  one  volume  increment 
s-as  assumed  to  represent  the  whole  s:  stem  volume. 

3.7.3,  A  chird  model  was  then  developed  irom  the  second  mode  •  whereby  t’  o  fluid 
was  completely  ujvideu  into  unit  volumes,  and  a  step  change  in  particle  con¬ 
centration  was  assumed  to  occtr  every  time  a  unit  volunc  left  the  reservoir.  iTcch 
line  volume  and  tne  actuator  and  reservoir  volumes  were  assumed  to  be  a  whole  number 
of  unit  volumes.  Plug  flow  was  assumed  in  all  pipework  and  equipment  except  the 
actuator  and  the  reservoir.  Changes  in  particle  cor centration  due  to  the  action  of 
equipment  were  assumed  to  be  instantaneous,  i.e.,  they  took  place  between  adjacent 
unit  volumes.  Since  plug  flow  is  assumed,  there  is  no  interchange  of  particles 
between  adjacent  volumes.  In  order  to  calculate  the  particle  concentrations  ir  the 
actuator  and  reservoir,  both  items  are  assumed  to  contain  perfectly  mixed  fluid. 

Step  changes  in  matrices  were  brought  about  by  considering  the  particle 
concentration  ol  each  unit  volume  entering  the  item  of  equipment  concerned. 

This  model  eliminated  the  two  major  difficulties  of  the  previous  models, 
and  also  handled  changes  in  tne  actuator  and  reservoir  volumes,  but  difficulties  in 
programming  led  to  a  slight  modification,  and  a  'ourth  model  was  finally  developed. 

3.7.f ,  In  this  model,  unit  volumes  of  the  pipework  and  equipment  were  considered 
instead  of  unit  volumes  of  fluid.  The  assumptions  of  plug  flow  in.  pipework,  etc., 
and  the  perfect  mixing  in  the  actuator  and  reservoir  still  held.  A  step  change  in 
every  volume  was  made  each  time  the  fluid  nad  completely  pa>  seu  from  one  unit 
volume  to  the  next.  The  unit  volumes  referred  to  ure  shcrenin  pig.  20  and  are 
referred  to  in  the  text  as  line  1,  etc. 

A  description  of  this  model  is  given  below:- 
3.7.5.  Assumptions: 

(1)  Changes  an  particle  concentration,  between  two  successive  volume 
increments,  are  s_,*.all.  This  justifies  the  use  of  a  step-wise 
model. 
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(2)  The  voluir.e  e  Darticles  is  small  compared  with  the  system 

volume,  and  :.rij  be  neglected. 

; 3 }  The  high  fluid  viscosity  and  small  particle  size  eliminates  particle 

settling  from  all  points  in  the  syscero. 

(4)  Tl.e  flow  regime  m  all  pipework  and  equipment  is  plug  flow,  except 

for  the  actuator  body  and  reservoir,  in  which  perfect  mixing  occurs. 

<5)  Ko  fluid  leaves  the  system  as  a  result  of  leakage, 

(6)  Volume  flow  rates  art  constant  for  each  of  the  three  valve 
positions. 

(7)  Acceleration  of  the  fluid  changing  from  ore  operating  state  to 
another  is  infinite,  i.e.. ,  the  change  of  state  involves  a  step 
change  in  velocity . 

3.1.6.  The  action  of  equipment  is  exactly  as  described  in  3.4  using  the  same 
nomenclature  as  in  Fig.  13  A  summary  of  the  equations  used  is  given  below, 
where  i  is  the  number  of  unit  volumes  which  have  left  the  reservoir. 


Outflow  from  reservoir: 


~  ill  _Li 


Pump  action: 


Filter  action: 


Fi  =  mi  +  bn+a 


£i  =  £i  hpi  bp.i 


Pressure  relief  valve: 


z.  =  Z.  Cfj  +  b , 


Control  valve  uction: 


(11  leakage  when  closed: 

C  .  =  L  .  ( f .  +  b  *f  b  ) 

ci  —i  -f  — h 

(2)  leakage  when  open: 

h  ■=  ici  <£,  *  i,  *  S»> 

(31  flow  of  high  pressure  fluid; 


Eq*  3.7.1 


Eq.  3,7.2 


Eq «  3.7.3 


Eq,  3.7.4 


Eq.  3,7.? 


Eq .  3.7.6 


=  <U 


+  bf  1  b^l  +  v,  +  bh) 


J  j  ’*■ 
,  1**  1 


VC  '  hi  ■  .if  T  tar1 


*  *  £  «-ho>  « *  s,  * 


Eq.  3.7.7 


(4)  flow  of  low  pressure  fluid: 


5 .  =  j  +  b. 
— 1  —  — k 


Eq.  3 ,,  7 . 8 


The  actuator  compositions  after  the  nth  time  it  is  filled  are: 


Side  A: 


...  .f±- 

~HTl  {  V  .  , 


rV  V 

_L_  C  (K  ‘iA_  kA 

Jjfl  JA  kA  --An"'*  ' 


ifV^A*  £q‘3‘7’9 


Side  B: 


+  C  b  +  Kk 
Jfl-kS  - 


Eq,  3.7.10 


Derivations  of  all  the  above  equations  are  given  in  3.4. 


The  basic  mass  balance  for  the  reservoir  is: 


r.  SL+  -  •  (Input) 

pi  P 


Eq.  3. 'Ml 


where  p  is  the  reservoir  volume. 


The  value  of  the  input  depends  upon  the  state  of  the  system. 


There  are  two  possible  values: 


fl)  When  flow  is  through  the  pressure  relief  valve: 
Input  =  (z  i  b  •*•  d) 


Eq.  3.7,12 


(2)  when  flow  is  through  the  control  valve: 


Input  =  (1  x  -r  d) 


Eq.  3.7.13 


The  value  of  1  depends  upon  the  valve  position,  i.e.,  whether  or  not  the 


actuator  is  keihg  Moved. 


However,  rh"‘;  the  actuator  is  cpeiated.  there  is  a  change  in  leservoir 


volume.  This  is  becous**  the  two  actuator  sides  have  been  different  volumes,  due  to 


the  piston  rod  in  side  v.  Thus  wher  side  d  is  filling,  the  reservoir  volume 


increases,  when  siae  A  is  fiilrng  the  reservoir  volume  decreases.  This  effect  may 


be  allowed  for  by  the  method  of  paragraph  3.7.7  below. 
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3.7.7.  When  the  side  B  is  filling,  the  increase  in  reservoir  volume  is  accomm¬ 
odated  by  increasing  the  input  at  regular  intervals,  proportional  to  the  ratio  of 
the  actuator  volumes.  Thus  if  the  ratio  of  the  actuator  volumes  is  5/4  then  every 
fourth  input  is  incxeased  to  (2^  +  2b^  +  d). 

However  d  is  not  doubled  since  it  depends  upon  the  -flow  rate  of  the 
reservoir  outlet. 

lhe  reservoir  volume  P is  then  increased  by  one  unit  before  calculation  of 
the  next  step. 

When  side  A  of  the  actuator  is  filling,  the  decrease  in  reservoir  volume  is 
accommodated  by  reducing  every  fifth  input  to  d  only,  and  decreasing  the  reservoir 
volume,  p  ,  by  one  unit  before  calculation  of  the  next  input. 

3.7.8.  It  is  assumed  that  failure  of  che  system  will  only  occur  in  the  valve,  and 
the  methods  of  3.5  are  used  to  calculate  this.  The  weights  of  parxicles  collected  in 
the  valve  are- 

(1)  behind  the  spool  (points  x^  and  x^): 

I  bu  =.  (1-$  )1(I-C>(f  b„  +  b)  Eq.  3.7.14 

— a  /c  —  —1  — h 

(2)  at  the  inlet  (point  x^): 

Z  bu  =  2(1"  e)$  (I-C  )  (f  +  b.  T  b  )  Eq.  3.7.15 

v- )  inside  the  spool  chamber  (points  x  and  x.,  or  x  and  x  ): 

‘x  o  b  / 

bu  =  4-  (I-C  )  (f  +  b,  +  b)  Eq.  3.7.16 

— f  vo  —  — f  — h 

Failure  in  other  equipment  can  be  handled  using  similar  techniques  to 
those  described  above. 

3.7.9.  The  line  filter  and  the  reservoir  filter  pick  up  particles  as  a  result  of 
their  action.  This  build  up  of  particle  can  effect  the  filter  efficiency  and  must  be 
calculated.  The  build  up  of  particle  after  n  unit  volumes  of  fluid  have  passed  through 
is  equal  to: 

(output  -  input)  Eq,  3.7.17 

o 
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3.60  The  computer  Plow  i  heet  for  In?  Final  iiodel 

3.8.1.  The-  hydr'.ulic  system  is  divided  into  si.<  lines  es  shown  in  Fig.  IS 

The  computer  flow  sheet  lor  predicting  failure  of  a  control  valve  due  t< 
particle  contamination  is  outlined  in  Fig.  19  to  Fig,  27 

3.8.2.  The  input  statement  consists  of  the  values  of  the  following: 


kA 


k3 


jA 


‘JB 


r 

R 

7. 

L 


jA 


JB 


These  are  system  /c  .tines. 


These  ere  all  matrices  which  descviDe  the 
action  cf  rhe  elements  in  +he  system  on  the 
particles. 
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k 


b 

— p 


5 


S  ’T 


"2 


b  i 


— kA 
— kB 


sbua' 

SBUC 

3UF 


Thepe  are  vacto-,  describing  particle 
concentrations;  a  refers  to  the  piston  rod, 
while  the  b's  refer  to  line  breakage. 


These  are  vectors  referring  to  build  up  in  the 
value,  their  initial  value  is  zero. 


This  is  the  number  of  increments  which  have  left 
the  reservoir,  initial  value  is  zero. 

These  are  parameters  used  in  calculating  reservoir 
volume,  initial  value  is  zero. 

Another  input  statement  is  the  value  of  the  particle  concentration  oi  each 

unit  volume  ir.  the  system.  I.*  the  system  is  to  be  modelled  from  start-up,  we  can 

assume  that  each  unit  volume  has  the  same  peptide  concentration  r  whe-c  r  the 

—o'  — o 

particle  concentration  of  the  clean  fluid. 


2.3.2.  The  first  step  in  the  calculation  is  in  calcu^at.*  the  particle  concentra¬ 
tion  of  nt  unit  volume  in  line.  1  a.te-  the  iluid  '•"■as  . -loved  one  ur._t  forv/r.rc, 

(Fig.  20).  The  new  particle  concent rn t -> on  oi  the  first  volume  is  then  calculated 

vEq •  3.7.1).  Other  unit  volumes  are  calculated  by  making  the  particle  concentration 

of  that  unit  volume  equal  to  the  particle  concentration  of  the  ^receding  unit 

volume,  except  when  the  fluid  has  to  ;ass  through  a  niece  of  equipment.  The 

particle  concentration  of  the  v  +  g  unit  volume,  i.e.,  the  first  unit  volume  after 

in 

the  pump,  is  made  equal  to  the  action  of  the  pump  on  the  narticlo  concentration  of 

the  v  unit  volume  (Fq,  3.7.2.) 
m 

The  process  is  then  continuous  until  unit  volume  v  *  v  +  1  is  reached 

p  m 

the  particle  concentration  of  this  unit  volume  is  made  equal  to  the  filter  notion  on 
the  particle  concentration  of  the  v  4-  v  unit  volume.  The  normal  sten  process  is 

p  !.! 

then  con  i.vueri  until  all  the  particle  concentrations  of  the  unit  volumes  in  Hue 
are  known. 

The  particle  build-up  in  the  lino  filter  and  reservoir  filter  is  then 
calculated  (Eq.  3.7.17).  This  is  done  as  the  efficiency  o*'  the  filter,  and  thus  the 
value  ov  .he  matrix  coefficients  is  dependent  on  the  build-up  in  the  filter 

2.3.3.  Depend ; r-g  an  the  value  of  I  the  programme  can  now  go  five  wavs .  (The  five 
ff  statements  <Fir.  20)  r.re  really  oi.e  but  are  shown  as  five  for  clarity). 

These  values  "HI  car  b.„  adjusted  to  simulate  any  operating  conditions  The  five 
paths  a re- 

fl)  Flow  through  the  pressure  leHcf  c-ive  in  the  continuous  circuit. 

(2;  Flow  due  to  leakage  through  the  -  mtrol  valve  in  tne  continuous 

circuit. 

(3j  Flow  due  to  tne  actuator  being  raise".,  i.e  ,  tne  i»«j ermittent 
c i rev 1 t, 

(4;  Flow  due  to  leaKige  . n rough  ths  control  -aiv;  nunnp  -ctuatc  •- 
operatior.. 

(5)  Flow  due  to  the  actuator  b’-ing  ,  owe  ret1 . 

Statements  20  and  20  deal  with  the  continuous  tircr.it  while  statements  40, 
50  and  60  refer  to  the  intermittent  circuit. 


Gd. 


%,-ssi  */*  jesKwisp 


3.G.4.  Statement  20  (Fig,  21)  accounts  for  flow  through  the  p’-ess.u.e  relief 

valve.  The  first,  step  is  to  calculate  the  particle  ccpceni ration  of  the  unit 

volumes  in  line  2.  This  is  done  by  continuing  the  normal  s.'ep  process  until  unit 

volume  v  -  v  +  v-r  v  ^  1  is  reached.  This  is  made  eaual  to  the  action  of  the 
m  p  f  g 

pressure  relief  valve  on  the  particle  con,  ei  tration  of  unit  -oliune  v  -*•  v  t  v  +  v  -i 

m  p  1  ? 

(Eq.  3.7.4).  The  new  value  of  £  is  then  calculated  \Eq  3.7  ii)  ,  I  is  increased  by 
£  and  the  process  repeated. 


3.8.5,,  Statement  30  (fig.  22)  deals  with  leakage  through  the  control  valve 
when  the  actuator  is  not  moving. 

The  first  operation  is  to  calculate  the  particle  concentration  of  the  unit 
volumes  in  line  3,  ana  also  tne  particle  concentration  of  the  first  unit  volume  in 
line  6.  This  must  be  -tone  as  the  fluid  in  line  6  can  be  obtained  from  either  line 
3,  1,  or  5.  The  values  in  ^ine  3  are  calculated  bv  tne  normal  step  process. 

'Ire  particle  concent rat ; on  of  the  first  unit  \ oiume  in  line  8  is  equal,  in 

this  case,  to  the  action  of  the  valve  in  the  particle  concentration  of  the 

v  +  v  +  v  +  v  unit  volume  in  line  3,  (Eq.  3.7.5).  Now  the  oartiole  concentra- 
m  p  f  h 

tion  of  the  v  -  v  -f  v  +  v  unit  volume  can  cause  blockage  of  the  valve  as 
n  p  f  h 

described  In  (3.5)  so  an  IF  statement  is  introduced  and  if  blockage  occurs  a  print 
out  is  cb coined. 


I;  no  blockage  occurs,  build  up  in  the  varve  in  positions 


arid  x 

3 


as  shown  in  Fig.  1!  is  then  calculated.  SBUA  refers  to  the  total  build  up  in 

position  Xj  and  which  will  prevent  the  spool  moving.  S3UC  refers  to  the  build 


up  in  position  and  may  cause  blockage  of  the  orifice.  Expression  for  calculating 
wne*l  ar  these  are  critical  have  been  discussed  in  3.5.2.  The  test  for  SBUA 
is  placed  before  the  test  lor  58UC  as  the  spool  has  to  move  before  fluid  can  flow. 


,f  neither  is  critical  the  nest  stage,  i.e.,  the  90  statement  (Fig. 

„  ,  ls  to  calculate  the  particle  concentration  of  the  unit  volume  m  line  , 
;MS  can  he  done  completely  on  the  step  basis  as  there  are  no  components  m  t  is 

line . 

1  lo  +  n  -  fFn  3  7  13)  an<i  move  the  process  on  one 
The  final  stage  is  to  calculate  (Eq.  '• 


stage. 
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3.8.6.  Statement  40  (Fig.  24)  describes  raising  of  the  actuator,  i.e., 

filling  of  Side  B.  The  iirst  step  is  the  calculation  of  the  unit  volume  of  line  3 

and  the  first  unit  volume  of  line  4.  Line  3  is  calculated  uc  above.  The  particle 

concentration  of  the  first  unit  volume  is  equal  to  the  valve  action  on  the  particle 

concentration  of  the  v  +  v  +  v„  +  v^  unit  volume  (£q.  3.7.7.)  Again  the  next  stage 

m  p  f  h 

is  to  ensure  that  this  particle  concentration  can  flow  through  the  valve  without 
causing  blockage.  If  it  can  pass  through  we  move  on  to  the  next  stage. 

The  first  unit  volume  of  fluid  through  the  valve  will  take  most  of  the 
build  up  in  position  Xg,  i.e.,  SBUC,  with  it,  and  no  other  unit  volume  will  take 
any  part  of  this  build  up.  Thus  the  first  unit  volume  will  take  the  lelt  hand  part 
of  this  build  up.  Thus  the  first  unit  volume  will  take  the  left  hand  paih  after 
the  second  IF  statement.  Here  the  first  step  is  to  calculate  the  particle  con¬ 
centrations  of  the  unit  volumes  in  line  4.  The  firsc  unit  volume  in  line  4  will  be 

given  by  the  valve  action  on  the  particle  concentration  given  by  the  addition  of 

SBUC 

the  narticle  concentration  of  unit  volume  v  +  v  +v„+v  to  (SBUC  - — —  ).  Where 

m  p  f  h  W 

is  'lie  particle  concentration  of  the  particles  left  in  position  x  .  The 
w  3 

program  now  moves  on  to  calculate  the  new  value  of  JE  (Eq.  3.7.10).  After  the  first 

volume  of  fluid  has  passed  through  the  valve,  the  other  volumes  take  the  righ*,  hand 

path.  Wnilst  the  actuator  is  being  is.  .”v2  leakage  still  occurs  causing  build  up  of 

particles  in  the  spool  chamber  at  points  x^.  and  x  .  (Fig.  15).  This  is  known  as 

BUF  and  is  the  particle  concentration  oi  the  build  up  due  tc  leakage  of  one  unit 

volume 0  The  first  unit  volume  through  the  valve  after  a  unit  volume  has  leaked  is 

assumed  to  carry  the  whole  of  DUF  with  it  into  the  actuator  line.  This  is  then  taken 

into  account  when  calculating  the  particle  concentration  of  the  unit  volumes  in  line  4 

(Eq.  3.7.6).  BUF  is  then  sot  equal  to  zero  since  there  is  nc-  collected  solid  reft 

at  points  x  and  x  until  after  another  unit  volume  has  leaked.  BUF  is  calculated 
4  5 

in  otatenent  50. 

The  next  stage  is  to  calculate  the  particle  concentration  of  the  unit 
volumes  in  line  5,  the  first  new  unit  volume  is  calculated  as  the  actuator  action 
on  the  particle  concentration  of  actuatcr  side  A,  i.e.,  (Eq,  3.7.9).  The  particle 
concentration  of  the  first  unit  volume  of  line  C  is  then  calculated  as  previously 
described  (Eq.  3.7.7). 

The  particle  concentration  is  then  tested  to  see  if  it  is  critical  for 
valve  blockage,  if  the  anser  is  negative  the  program  moves  on. 


Another  IP  statement  now  follows  to  account  for  changes  in  the  reservoir 
volume  as  described  in  Eq.  3.7.11.  In  this  program  the  ratio  of  the  actuator 

volumes  V3A  is  taken  as  5/4.  If  parameter  S  is  a  multiple  of  four,  two  unix 

vjB 

volumes  are  added  to  the  reservoir  while  only  one  unit  volume  leaves  and  Statement  70 
deals  with  this.  If  S  is  not  a  multiple  of  four  there  is  no  difference  in  the 
reservoir  volume  input  and  output. 

In  Statement  70  (Fig.  25)  the  first  step  is  to  calculate  the  particle 
concentration  of  the  unit  volumes  in  line  6e  This  can  be  done  on  a  step  basis  as 
the  particle  concentration  of  the  first  unit  volume  of  line  6  is  already  known. 

The  particle  concentration  of  the  final  unit  volume  is  then  set  equal  to  the  particle 
concentration  of  a  dummy  volume  Y, 

Since  at  S  =  4  two  unit  volumes  leave  the  actuator  for  only  one  leaving 
the  reservoir ,  the  fluid  in  lines  5  and  6  is  moved  on  two  unit  volumes  so  the 
particle  concentration  of  the  unit  volume  in  line  5  and  the  first  unit  volume  in 
line  6  is  then  recalculated.  The  test  for  critical  valve  concentration  then  takes 
place  and  if  the  answer  is  negative  the  particle  concentration  of  the  unit  volume 
in  line  6  is  then  calculated  as  previously  described. 

The  reservoir  composition  is  then  calculated  remembering  that  the  input  to 
the  reservoir  is  the  unit  volume  Y  and  the  final  unit  volume  of  line  6,  the  particle 
concentration  of  these  volumes  has  already  been  calculated  (Eq.  3.7.11).  The  new 
reservoir  volume  is  then  calculated  by  adding  one  unit  volume  to  the  old  reservoir 
volume.  The  program  is  then  moved  on  one  stage  and  returned  to  10. 


go  to 


<~oapuier  procr.ua 
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0.8.7.  Statement  50,  (Fig.  26)  deals  with  leakage  through  the  valve  while 
there  is  fluid  flow  through  the  valve,  i.e»,  actuator  operation.  Leakage  is 
assumed  to  occur  at  a  certain  value  of  I.  The  first  calculation  is  of  the  particle 
concentration  of  the  unit  Vorutjes  in  line  3  as  described  before,  and  the  particle 
concentration  of  the  1st  unit  volume  in  line  £  (Eq.  3.7.6).  Again  a  test  far 
critical  vaive  concentration  takes  place  on  the  inflow  to  the  valve.  The  next 
calculation  is  of  the  valve  buiid  up  BUA  and  BUF  as  described  in  3.4.4.  SBOA  is 
then  calculated  and  tested  to  see  if  it  will  stop  the  spool  moving,  if  a  negative 
answer  is  obtained  the  program  moves  on  to  statement  90. 


3.8.8.  Statement  60  (Fig.  2?)  deals  with  fluid  flow  due  to  lowering  of  the 
actuator.  When  the  jack  is  lowered  side  A  is  filled  and  thus  there  is  a  decrease 
in  the  volume  of  the  reservoir.  The  first  calculation  is  to  calculate  the  particle 
concentration  of  the  unit  volume  in  line  3,  as  previously  described,  and  the  first 
volume  of  line  5,  in  a  similar  manner  to  the  way  the  particle  concentration  of  the 
first  unit  volume  in  line  4  was  calculated  in  statement  *10  (Eq.  3.7.7). 

The  test  for  critical  value  concentration  then  takes  place  if  the  answer 
is  negative  the  program  moves  on  to  another  IF  statement.  This  IF  statement  is 
included  to  handle  particle  build  up  in  the  valve,  and  is  of  a  similar  fora  to  that 
in  statement  40.  The  ne:;t  calculation  of  the  particle  concentration  in  actuator 
side  A,  i.e.,  (Eq.  3.7.9). 

The  next  set  of  statement  ir.  designed  to  handle  the  change  in  reservoir 
volumes.  As  we  assumed  that  the  ratio  cf  the  actuator  volume  is  °/4 ,  we  look  at 
the  value  of  parameter  T.  If  T  is  not  a  multiple  of  4  we  calculate  the  particle 
concent  ration  of  the  unit  volume  in  line  4,  and  the  first  unit  volume  in  line  6  in 
a  similar  manner  to  the  calculation  used  in  statement  40. 


Calculate  Cone .  of 
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If  T  is  a  multiple  of  four,  we  move  to  statement  80  (Fig.  28).  This 

calculates  the  new  reservoir  composition  (Eq.  3.7.jlI),  decreases  the  reservoir 
volume  by  one  unit  volume,  advances  T  and  I  by  one  and  goes  back  to  10. 

3.8.9.  The  above  program  will  run  until  one  of  the  critical  tests  is  positive. 

It  will  then  print  out  why  failure  was  caused,  the  value  of  I,  the  number  of  times 
the  actuator  circuit  is  used,  the  reservoir  concentration  and  both  filter  build-ups. 
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Computer  program  flowsheet  -  St clement  SO 


3.9 


Conclusions  and  Suggestions  for  Further  Work 


3,9.1.  The  mathematical  model  for  the  hydraulic  system  is  presented  in  the  report 
in  computer  flowsheet  form.  A  computer  program  based  on  the  flowsheet  would  be 
able  to  predict  failure  in  the  control  valve.  It  could  easily  be  modified  to 
account  for  additional  equipment  in  the  circuit,  and  routine  replacement  of  worn 
equipment  by  new  equipment.  It  could  also  be  used  to  predict  failure  in  other  items 
of  equipment. 

There  are,  however,  several  difficulties  in  writing  the  program  and  these 
are  tabulated  below.  Most  of  then  are  essentially  lack  of  available  information  for 
the  input  data. 

The  types  of  information  i-eauired  are: 

Details  of  the  volumes  of  components  within  a  system,  the  initial  particle 
concentrations ,  the  frequencies  cf  actuator  operation  and  the  number  of  operations, 
information  on  modes  of  equipment  failure  and  above  all  experimental  studies  to 
determine  the  matrix  coefficients  for  individual  components.  To  do  this  would 
require  detailed  studies  of  the  input  and  output  of  particles  through  given 
components  under  controlled  conditions.  Probably  the  most  useful  means  of  obtaining 
this  information  is  to  construct  a  test  circuit  and  make  experimental  observations. 
This  procedure  must  bo  used  to  obtain  the  initial  values  of  the  matrix  coefficients 
since  no  theoretical  work  (to  our  knowledge)  has  been  done  towards  calculation  of 
these  coefficients  from  the  system  parameters.  This  in  itself  would  provide  a 
heavy  program  of  further  work. 

The  complexity  of  the  computer  flow  sheet  would  suggest  that  some  work 
towards  determining  the  comparative  magnitude  of  the  effects  in  each  part  of  the 
syscem.  If  this  information  was  available,  the  flow  sheet  could  probably  be 
simplified  to  a  large  extent,  reducing  considerably  the  amount  of  computer  storage 
required  to  run  the  program. 
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Chapter  4 


The  Characterisation 


of  Particles 
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4.1 


Introduction 


The  behaviour  of  a  particulate  system  depends  on  the  properties  of  the 
individual  particles  which  go  to  make  up  that  system  and  also  on  their  spatial 
arrangements  with  respect  to  each  other.  It  is  often  said  that  particulate  systems 
are  particularly  complicated  because  of  the  large  number  of  variables  involved. 
However,  the  variables  are  usually  not  independent  and  the  basic  properties  of  a 
particulate  system  are  very  few  in  number. 

The  easiest  way  to  characterise  a  powder  is  to  measure  some  macroscopic 
property  of  the  powder  and  to  predict  its  behaviour  from  that  measurement.  If  the 
pcwder  cannot  be  changed  then,  for  engineering  purposes,  this  kind  of  measurement 
may  be  sufficient.  If,  however,  it  is  desired  to  change  the  macroscopic  properties 
j  in  a  controlled  manner  or,  alternatively,  to  find  the  best  possible  system,  then 

this  must  be  done  by  relating  the  macroscopic  properties  of  the  powder  back  to  its 
primary  properties.  However,  rarely  does  the  behaviour  of  a  powder  have  a  unique 
relationship  with  its  primary  properties  since  the  same  powder  can  exist  in  different 
states.  Unless  the  particles  are  so  well  dispersed  that  they  behave  as  individual 
particles,  the  spatial  arrangement  of  the  particles  must  also  be  considered. 

Thus  the  characterisation  of  a  particulate  system  is  summarised  in  Table  IV 
the  macroscopic  properties  listed  being  only  examples  of  the  many  which  can  be 
measured  or  possibly  predicted  from  the  primary  properties. 


TABLE  IV 

Relationship  of  properties  of  particulate  systems 


Primary  Properties 

State 

Macroscopic 

Properties 

Behaviour 

Particle  Size 

Flow  of  Fluid 

Distribution 

Porosity 

j 

Permeability 

Anisotropy 
&  Homogeniety 

Shear  Strength 

Flow  of  Powder 

Properties  of 

Solid  and 

Interstitial 

Moisture 

Flow  of  Heat 

Fluid 

Content 

Surface 

Properties 

Conductivity 

and  Electricity 

The  basic  difference  in  the  behaviour  of  a  subdivided  solid  and  a  homoge¬ 
neous  solid  is  the  interaction  of  the  constituent  particles  and  this  is  dependent  on 
the  distribution  of  size  and  shape  of  the  particles.  The  particle  size  distribution 
of  the  particles  is  thus  the  primary  property  of  a  powder  which  differentiates  it 
from  homogeneous  solid.  If  the  size  and  shape  of  each  particle  in  a  sample  of 
powder  were  known,  then  these  measurements  would  be  sufficient  to  differentiate  it 
from  any  other  powder.  However,  they  would  not  be  sufficient  to  predict  the  macro¬ 
scopic  properties  of  the  system,  since  the  same  set  of  particles  can  be  arranged 
geometrically  in  different  ways  with  respect  to  each  other  to  produce  completely 
different  behaviour.  Thus,  the  porosity  of  the  bed  must  be  known  and  also  the 
moisture  content  of  the  powder  if  a  three  phase  system  is  involved.  A  further 
complication  is  added  in  that,  although  the  porosity  of  a  bed  does  not  uniquely 
depend  on  the  particle  size  distribution  of  the  particles,  they  are  nevertheless  not 
completely  independent.  A  particular  set  of  particles  in  the  packed  state  has  upper 
and  lower  limits  of  porosity  which  are  a  unique  function  of  the  particle  size 
distribution.  The  actual  porosity  may  lie  between  these  two  values  depending  on  the 
manner  in  which  the  particles  have  been  packed.  Under  some  circumstances,  the 
porosity  of  the  bed  may  be  a  unique  function  of  the  particle  size  distribution.  For 
example,  if  powder  is  flowing,  then  it  must  be  at  its  critical  porosity. 

The  complexity,  then,  in  describing  the  behaviour  of  a  particulate  system  is 
usually  because  several  mechanisms  are  involved.  An  attempt  to  correlate  the 
behaviour  of  a  powder  directly  to  its  particle  size  distribution  is  usually  not 
possible  and  the  mechanisms  must  be  separated  and  related  individually  to  the  size 
distribution,  taking  into  account  the  state  of  the  system.  Illustrating  this  point 
with  Table  IV  the  behaviour  of  a  powder  can  only  be  related  to  the  particle  size 
distribution  by  passing  through  the  two  other  columns  of  macroscopic  property  and 
state. 
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4.2 


Particle  Size  Distribution 


The  traditional  methods  of  defining  particle  size  is  as  the  size  of  a 
sphere  equivalent  in  some  geometrical  property  such  as  area  or  volume  or  equivalent 
in  behaviour  such  as  settling  speed.  The  particle  shape  is  then  inherent  in  the 
definition  of  size.  Shape  factors  are  then  the  ratios  of  the  various  defined  sizes, 
either  for  one  particle  or  an  average  shape  factor  for  a  system. 

The  alternative  definition  of  particle  size  is  a  statistical  diameter. 

Well  known  examples  are  those  due  to  Martin  and  Feret. 

In  addition  to  the  problem  of  defining  tne  size  of  an  individual  particle, 
there  is  a  further  problem  of  expressing  the  distribution  of  particle  size,  for 
example  by  number  or  by  weight.  If  particle  size  is  defined  as  an  equivalent  sphere, 
then  the  moments  of  the  various  distributions  which  can  be  plotted  have  significance 
and  may  sometimes  be  related.  On  the  other  hand,  if  the  Martins  or  Feret  diameter 
is  used,  the  physical  significance  of  the  distributions  is  more  obscure. 

It  is  obvious  chat  for  most  aspects  of  particle  behaviour,  there  is  no 
single  equivalent  size  of  spherical  particle  to  which  the  behaviour  of  an  irregular 
particle  can  be  uniquely  related.  Nor  even  can  a  single  statistical  diameter 
contain  all  the  information  about  the  geometry  of  a  particle.  The  limitation  of  the 
approaches  is  that  of  imagining  a  single  particle  to  be  the  primary  unit  of  a 
powder  system.  An  alternative  set  of  measurements  will  now  be  proposed  which 
completely  characterise  a  particle  system. 

Consider  an  irregular  particle  as  shown  in  Fig-,  29.  That  particle  can  be 
imagined  to  be  split  into  laminae  which  are  of  finite  but  small  thickness.  The 
shape  of  each  laminae  is  irregular  but  each  has  a  definite  area.  Therefore,  a 
number  distribution  curve  of  the  area  of  the  laminae  can  be  plotted.  If  the 
orientation  of  the  particle  is  changed  and  the  process  repeated,  then  a  different 
set  of  laminae  will  be  produced.  If  the  laminae  in  every  possible  orientation  or 
direction  are  measured,  then  they  may  all  be  added  together  to  produce  a  distribution 
curve  which  is  characteristic  of  the  particle. 

This  distribution  is  a  unique  function  of  the  particle  but  the  reverse  is 
not  necessarily  so.  Thus,  two  different  particles  could  conceivably  produce  the 
same  distribution  curve.  The  laminae  may  now  be  broken  down  again  into  smaller  units. 


FIGURE  29  Characterisation  of  a  single  particle 


Each  lamina  will  bo  split  into  filaments  which  have  a  square  cross-section  of  small 
but  finite  value.  These  filaments  vary  now  only  in  their  length  and  a  distribution 
curve,  again  by  number,  of  their  lengths  is  a  characterisation  of  the  particle.  The 
filaments  may  be  taken  in  one  direction  or  alternatively  they  may  be  taken  in  all 
possible  directions  and  these  various  distributions  added  together  to  produce  a 
compound  distribution.  This  distribution  must  then  be  normalised  so  that  the  first 
moment  of  the  distribution  is  equal  to  the  volume  of  the  particle  divided  by  the 
cross-sectional  area  of  a  filament. 

If,  instead  of  considering  a  single  particle  a  powder  is  considered,  the 
numoer  distribution  of  laminae  and  filaments  can  be  plotted  as  a  characterisation  of 
the  whole  system.  A  further  distribution  must  be  added  in  order  to  make  this 
characterisation  unique  and  that  is  the  number  to  volume  distribution  of  the 
particles. 


It  is  proposed,  therefore,  that,  a  system  of  particles  is  completely 
characterised  by  three  distributions :- 


number 

number 

number 


volume  of  particles 
area  of  laminae 
length  of  filaments 


In  the  limit,  r.s  the  thickness  of  the  laminae  and  filaments  tends  to  zero, 
the  distributions  transpose  from  a  histogram  to  a  continuous  curve  and  the  laminae 
and  filament  distributions  become  area  and  chord  distributions  respectively. 

This  characterisation  of  a  particulate  system  is  illustrated  in  Fig.  30 
and  is  now  a  unique  characterisation.  Although  two  different  sets  of  particles 
could  produce  the  same  distribution  of  one  or  even  two  of  these  variables,  if  all 
three  distributions  are  the  same  for  two  samples,  then  they  are  tne  same  powder. 


The  curves  must  now  be  normalised  such  that  the  first  moment  of  the  number 
to  volume  curve  is  equal  to  some  convenient  value.  This  will  probably  be  either 
unity  or  the  volume  of  particles  per  unit  volume  of  bed.  The  area  and  chord 
distributions  must  then  be  normalised  such  that  the  first  moment  of  the  distribution 
is  equal  to  the  same  normalising  voluse.  This  point  is  better  illustrated  ip 
Fig.  31  where  each  lamina  has  been  taken  and  they  have  been  shuffled  until  they 
stand  in  order  of  ascending  size.  Then  this  curve  in  fact  represents  the  number 
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distribution  curve  of  laminae  plotted  as  a  cumulative  distribution.  The  area  under 
the  curve  is  the  volume  of  the  particles  divided  by  the  thickness  of  one  lamina 
and  is  also  equal  to  the  first  moment  of  the  primary  distribution  curve.  The  same 
argument  applies  to  the  distribution  of  filaments  as  illustrated  in  Fig.  31 

Thus:- 

n(v)  v  dv  =  Y 


f 


J 

all  v 


n(a) 


da 


all 


nf£)  d 


where 


t-t 


V 

At 

y_ 

(At)2 

thickness  of  o  cherd  or  filament 
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4.3 


Pore  Size  Distribution 


The  pore  space  between  the  particles  is  a  continuum  and  is  not  sub¬ 
divide!  like  the  particles.  It  can,  however,  be  characterised  in  the  sa-ao  manner  as 
the  particles  by  splitting  it  into  filaments  of  small  but  finite  cross-section. 

Each  filament  will  terminate  cn  a  particle  surface.  This  is  illustrated  in  Fig.  32 
Tne  number  distribution  of  filament  lengths  is  then  a  characterisation  of  the  void 
Space  and  it  must  be  normalised  to  a  total  length.  Usually  this  will  be  the  total 
length  of  filament  in  unit  volume  of  ted.  Thus:- 


\ 

J  n(k)  k  dk 
J  all  k 


1  -  - 
(A  t2) 


The  filaments  may  be  taken  in  one  direction  or  may  be  taken  in  ail  possible 
directions.  If  the  particles  are  touching,  the  distribution  curve  must  oass  through 
the  origin,  lair  in  a  disperse  sj :*ea  this  is  not  nece ssarily  so. 

The  mean  filaueni.  size  of  the  pores  is  easily  related  to  that  of  the 
particles  if  the  bed  is  both  random  and  isotropic.  Thus,  if  „i  long  straight  line  is 
passed  through  the  bed  of  particles  then  it  samples  both  the  particle  and  void 
filaments.  The  fractional  length  of  the  line  occupied  by  the  void  filaments  is  equal 
to  the  voidage  cf  the  bed  and  the  total  number  of  solid  filaments  equals  the  total 
number  of  void  filaments.  This,  of  course,  includes  the  v'id  filaments  of  zero 
length  "hen  two  particles  touch.  This  fact  has  been  used  for  some  time  by  met¬ 
allurgists,  biologists  and  geologists  to  measure  the  fractional  volume  o* 

constituent  parts  by  microscopic  traversing  techniques. 

The  relationship  is  then  as  follows:- 


bet  1 
n(l) 
k 

n(fc) 


be  a  particle  filament 

be  the  distribution  function  by  a  number  of  filaments. 

be  the  length  of  a  void  filament 

to  be  the  distribution  function  of  void  filaments. 


Since  in  any  total  length  cf  line 
th»>  number  cf  void  filaments. 


*i 


r.{l)dl 


n(k)dl 


the  number  of  solid  filaments  is  equal  to 
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mount;  33  Partlclo  and  Voidago  filaments 


Since  the  fractional  length  of  void  £i laments  is 


nCDldl 


l  -  S' 


J  n(k)kdk 

Jn(l)ldl 

J  n(l)dl 

=  1 
tr 

J  n(k)kdk 

J  n(k)dk 

k  =  T 

s: 

1 

-  51 

where  the  bar  indicates 

the  Bean  values  by  number. 

Description  of  Anisotropy 


4,4, 

If  the  bed  is  isotropic,  then  the  choid  size  distribution  of  both  the 
particles  and  the  void  space  will  be  the  same  ill  all  directions.  However,  in  an  an¬ 
isotropic  bed,  the  distributions  may  be  measured  in  two  or  more  directions  depending 
on  the  problem  to  be  solved.  Thus,  for  example,  in  calculating  the  flow  of  elec¬ 
tricity  through  the  solid  particles  or  the  flow  of  liquid  through  the  pore  space,  it 
would  be  necessary  to  know  the  distributions  in  the  direction  of  flow  and  also  those 
perpendicular  to  the  flew. 

It  should  be  noted  that  a  bed  of  particles  can  be  anisotropic  due  to 
differences  in  the  shape  of  the  filament  distributions.  Physically  this  means  that 
the  particles  are  preferentially  orientated.  However,  the  fractional  space  occupied 
by  void  and  by  particle  is  the  same  regardless  of  direction  and  the  ratio  of  the 
mean  values  is  always  given  by  equation  (1).  The  distributions  of  an  anisotropic  bed 
are  illustrated  in  Fig.  34 


93. 


w«»W*®sm 


2.  PARTICLE  SPACING 


4.5 


Characterisation  of  Agglomerates 


If  the  particles  are  agglomerated  together  into  clumps,  then  the  behaviour 
of  the  system  is  affected  and  further  information  must  be  added  to  the  description 
of  the  system.  This  is  done  by  expressing  the  number  distributions  of  volume,  area 
and  chord  which  describe  the  agglomerates  as  well  as  those  for  the  individual 
particles.  There  are  now  also  two  porosities  which  are  needed  to  describe  the  state 
of  the  particles.  The  porosity  of  an  individual  agglomerate  describes  the  packing 
of  the  particles.  If  the  agglomerates  are  now  considered  to  be  solid,  then  the 
fractional  volume  of  the  space  between  them  describes  their  spatial  arrangement. 

This  is  illustrated  in  Fig.  35. 
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Conclusion 


(a)  A  system  of  particles  can  be  completely  characterised  by  the 
following  three  distributions 

number  :  length  of  chords 

number  :  area  of  areas 

number  :  volume  of  particles 

(b)  The  arrangement  of  the  particles  can  be  expressed  as  the  porosity 
.  of  a  bed  or  concentration  of  a  suspension  if  the  arrangement  is 

isotropic  and  random.  If  directional  or  spatial  variation  exists, 
the  distributions  also  have  to  be  measured  in  different  directions 
or  different  localities. 

(c)  Although  without  rapid  automatic  techniques  this  method  of  analysis 
is  not  suitable  for  routine  work,  it  does  enable  relationships  to 
be  established  between  the  conventional  particle  size  analysis 
parameters  and  properties  of  the  particle  system.  Work  is  in 
progress  on  the  mathematical  nature  of  these  relationships. 
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Chapter  5 


Particle  Size  Analysis 
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5.1  Introduction 


Many  methods  of  measuring  the  size  of  fine  particles  have  been  advocated 
and  the  principal  ones  are  listed  in  Table  V. 

TABLE  V  Techniques  for  measuring  particle  size 


Glass  of  ! 

Technique 

Technique 

Approx,  size  range 
commonly  measured 
(|f  m) 

Parameter 

Measured 

Field  Scanning 

Optical  Microscopy 

>  0.5 

Longest 

dimension 

equivalent 

area  dia. 

statistical 

dias. 

Electron  Microscopy 

10"3-  5 

as  optical 

Scanning  Electron 

Microscopy 

0*05  -  50 

-3 

mic^pscopy 

Light  Scattering 

10  -1 

surface  area 

Stream  Scanning  ! 

Coulter  Counter 

1-100 

volume 

Hiac  Counter 

2  -  2000 

projected  area 

Light  Scattering 

0.5  —  oO 

surface  area 

Classification 

Sieving 

>  5 

sieve  diameter 

Air  classifiers 

2-50 

Stokes  diameter 

Sedimentation 

Pipettes 

1-60 

Stokes  diameter 

Balances 

1-60 

Stokes  diameter 

Photosediment- 

ometers 

0.2  -  60 

Stokes  diameter 

Centrifugal  methods 

0.05  -  10 

Stokes  diameter 

Surface  Area 

Gas  Adsorption 

Permeability 

>  lOm^g  1 

2  -1 

>  50m  g 

The  parameters  listed  in  the  final  column  are  only  equal  for  spherical 
particles.  For  irregular  particles  their  ratios  are  characteristic  of  the  shapes 
of  the  particles  although  these  relationships  are  very  complex.  Several  workers, 
notably  Heywood  (Chem.  and  Ind.  1937,  56,  149  and  other  papers),  have  defined 
shape  factors  in  terms  of  the  ratios  of  principal  dimensions  of  particles.  This 

method,  whilst  giving  a  very  good  description  of  the  components  of  a  system  doe3 
not  permit  of  an  analytical  synthesis  of  the  system  properties  from  the  croper ties 
of  these  individual  components.  One  such  possible  method  is  described  in  Chapter 
4  of  this  report  where  the  particles  are  considered  as  being  made  up  of  more 
elementary  one  and  two  dimensional  units  which  can  be  expressed  in  terms  of 
distribution  functions. 

Of  the  methods  of  size  analysis  listed  in  Table  V  many  can  be  eliminated 
by  the  nature  of  the  problem  in  analysing  particles  in  mineral  based  hydraulic  oil. 
The  characteristics  of  this  particular  system  are: 

(i)  Low  concentration  of  particles 

(ii)  Size  range  >100  -  5p  a  (and  possibly  less) 

(iii)  A  mineral  based  fluid  containing  several  additives 

(iv)  A  range  of  material  densities 

(v)  A  suspending  fluid  of  relatively  high  viscosity. 

Considering  these  properties  in  relation  to  the  methods  listed  in  Table  V  , 
several  techniques  may  immediately  be  eliminated,  e.g. 

Methods  Undesirable  Characteristic  above 

Sieving  (ii) 

Sedimentation  (i),  (ii),  (iv),  (v) 

Surface  area  generally  impracticable 

Air  classifiers  "  " 

This  leaves  the  field  and  stream  scanning  techniques.  The  standard 
methods  in  use  today  are  manual  field  scanning  by  optical  microscopy  (e.g. 

ARP  598  and  I.P.  (draft)  method  (U.K.)  and  surveys  have  shown  these  methods  to 
be  of  less  than  desirable  accuracy  and  extremely  time  consuming  and  fatiguing. 
A.R.P.  598  is  discussed  in  5.2  below.  Automatic  scanning  systems  have  been 
developed  to  remove  the  manual  effort  from  this  type  of  analysis  and  all  current 
instruments  are  of  the  spot  scanning  type.  These  include  the  Millipore/Bausch  & 
Loab  TIHC  (U.S.A, )  counter  which  measures  a  parameter  approximating  to  the 


maximum  length  specified  in  ARP  598  and  I.P.  (draft)  methods,  the  Quantimet 
Image  Analysing  Computer  (U.K.)  developed  originally  for  inetal3ographic  analysis 
and  German  instruments  manufactured  by  Laitz/Bosch  and  Zeiss.  Ail  of  these 
instruments  rely  on  a  television  camera  to  scan  the  specimen  and  have  digital  or 
analogue  computing  facilities  of  varying  complexity .  The  majoi  problem  in 
designing  instruments  of  this  type  are  grey  level  discrimination,  re-entrant 
particles  and  touching  particles.  The  latter  is  a  very  difficult  problem  beyond 
the  current  range  of  commercially  available  machines. 

For  a  more  detailed  analysis  a  more  complete  collection  of  data  is  required 
and  any  machine  making  use  of  a  television  camera  producing  up  to  say  25,000,000 
items  of  information  per  second  is  operating  too  fact  for  most  practical  data 
logging  systems.  This  type  of  analysis  requires  a  mechanical  or  flying  spot  scanning 
system  and  as  such  is,  together  with  its  computing  facilities,  both  slower  and  more 
expensive  than  the  current  type  of  equipment.  An  outline  of  the  operation  of  one 
system  is  described  below  in  5.5. 

The  best  known  example  of  a  stream  scanning  particle  size  analyser  is  the 
Coulter  Counter.  In  a  stream  scanning  apparatus  particles  are  passed  one  at  a 
time  through  a  sensing  zone  which  produces  a  signal  which  is  then  recorded  and 
interpreted  in  terms  of  particle  size.  In  the  Coulter  Counter  this  transducer  is 
an  orifice  through  which  an  electrically  conducting  fluid  containing  the  particles 
is  drawn.  A  potential  is  applied  across  the  orifice  and  the  change  in  electrical 
resistance  due  to  the  presence  of  a  particle  measured.  Section  5.4  describes 
a  detailed  investigation  of  this  device  and  illustrates  some  of  the  general 
problems  to  be  investigated  in  any  stream  scanning  device.  Because  of  the 
requirement  of  an  electrically  conducting  fluid,  the  Coulter  Counter  is  not  suit¬ 
able  for  Mil-H  5606  systems.  Although  mineral  oils  can  be  made  electrically 
conducting  by  addition  of  certain  suitable  solvents  and  electrolytes  this  is 
extremely  difficult  with  Mil-H .5603  probably  due  to  the  additives  present.  During 
a  short  but  intensive  survey  no  completely  satisfactory  electrolyte  system  was  found 
for  Shell  Aero  Fluid  No.  4,  DTD  585  (approx.  U.K.  equivalent  of  Mil-H. 5506)  and  this 
impression  was  confirmed  by  discussion  with  a  representative  of  the  manufacturer. 

An  instrument  which  is  becoming  increasingly  used  fcr  hydraulic  fluid 
analysis  is  the  HTAC  Counter  (High  Accuracy  Products  Corpn. ,  Clafemont,  Calif.) 

in  which  the  amount  of  light  obscured  by  a  particle  crossing  a  light  beam  is 
measured.  It  is  assumed  that  the  amount  of  light  obscured  is  proportional  to  the 


cross  sectional  area  of  the  particle.  A  brief  report  on  this  instrument  appears 
in  5.3  below,  Because  it  can  bo  used  with  the  hydraulic  oil  as  the  sole  suspend¬ 
ing  fluid  it  Las  possibilities  for  use  in  hydraulic  systems  as  an  on-line  real 
time  analyser. 

Many  other  forms  of  stream  scanning  counter  have  bean  describee  making 
use  of  electrical  capacitance,  light  scattering,  ultrasonic  reflection  or  atten¬ 
uation  amongst  the  properties  studied.  The  production  of  one  of  these  devices 
in  a  form  suitable  for  use  in  an  operating  hydraulic  system  would  involve  an 
intensive  program  of  development  work. 

One  problem  to  be  faced  in  defining  the  size  of  particles  is  that  presented 
by  aggregates.  For  example  it  is  known  that  Mil-H.5606B  fluid  that  was  within 
specification  when  packed  may  be  found  to  contain  particles  >100um  after  storage 
and  that  these  particles  are  floes  of  smaller .primary  particles.  Although  these 
particles  will  be  counted  by  ARP  598,  their  resistance  to  shear  would  be  very  low 
and  they  would  not  be  expected  to  manifest  themselves  under  operating  conditions. 
If  this  problem  were  of  any  magnitude  in  approving  fluids,  it  might  be  advisable 
to  introduce  the  applicati  on  of  some  shear  work  to  the  fluid  prior  to  ar  > lysis 
although  how  this  could  be  done,  apart  from  vigorous  agitation  of  the  sealed 
container,  without  introducing  more  contaminant  is  difficult  to  see.  Thin  problem 
might  be  less  significant  with  stream  scanning  devices  where  shear  occurs  during 
the  course  of  analysis. 
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5.2 


ARP  598 


re;  Procedure  for  the  determination  of  Particulate  Contamination  of 
Hydraulic  Fuels  by  the  Particle  Count  Method 

Objective 

(1)  To  investigate  the  validity  of  the  assumption  implicit  in  the  procedure 
that  the  particles  counted  over  small  number  of  randomly  chosen  grid 
squares  could  be-  used  to  express  the  total  particle  count. 

(2)  To  evaluate  critically  the  experimental  pronsdure  of  AT,?  598. 

Results 

A  closely  sized  fraction  (Alpine  M2  classifier)  of  A.(J.  Fine  test  dust  in  the 
range  35  -  cOym  was  prepared  and  a  suspension  of  the  particl-ss  in  Aero-3heii  4 
hydraulic  fluid  was  prepared.  Samples  of  this  suspension  were  filtered  following 
the  procedures  of  ARP  596.  Two  uncertainties  were  found  in  the  method. 

(1)  Instruction  8. 2. 2. 5  is  a  little  misleading  as  it  did  not  indicate  v/heiher 
the  petroleum  ether  should  be  layered  on  the  oil  or  whether  some  mixing 
should  be  obtained  by  pouring  in  the  petroleum  ether  more  vigorously. 

(2)  Particles  were  observed  to  deposit  on  the  shoulders  of  the  filter  funnel. 

No  quantitative  tests  were  made  of  the  loss  incurred  but  it  was  observed 
that  with  particles  of  the  size  used  only  vigorous  washing  would  dislodge 
them. 

The  results  of  4  separate  experiments  are  presented  in  Tables  VI 
IX.  Care  was  taken  with  each  experiment  to  locate  the  membrane  filter  grid  in 

the  position  with  respect  tc  the  glass  sinter  substrate.  Thus  the  counts  in  each 
square  can  be  added  and  the  result  of  this  addition  is  given  in  Table  X  where 
only  the  complete  grid  squares  are  represented.  The  mean  particles  for  a  square 
is  108  with  a  standard  deviation  of  13.  It  should  be  noted  that  only  three 
results  are  outside  the  range  of  mean  size  +  2  x  standard  deviation  and  all  are 

within  +  3  x  standard  deviations.  This  is  consistent  with  a  random  distribution 
of  particles. 

Table  XI  shows  in  each  corner  of  each  square  the  number  of  counts  x  100 
divided  by  the  mean  for  that  experiment.  No  systematic  high  counts  or  low  counts 
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are  apparent.  If  present,  systematic  high  or  low  counts  might  indicate  heder- 
ogenerties  in  the  porous  membrane  support. 

Conclusions 

Since  no  control  was  maintained  on  the  concentration  of  the  suspension  used 
in  these  experiments,  it  has  not  proved  possible  to  apply  variance  analysis  between 
the  experiments.  The  indications  are  that  on  the  apparatus  used  the  particles  are 
distributed  at  random. 

Two  further  points  should  be  noted; 

(1)  The  approximate  standard  deviation  on  a  number  of  counts  N  is  ^N~and  for 

small  N,^| is  large  (for  N  =  100, ^^T"=  10).  Since  at  95%  confidence 

limits  there  can  be  a  acceptable  range  of  N  -  in  the  number  of  counts 

or  +  i.e.  a  total  range  of  4<^t(e.g.  40%  in  the  case  above).  Thus 

N  N 

the  scale-up  from  counting  the  number  of  particles  on  a  few  squares  may 
lead  to  large  over  or  under  estimates  of  the  true  count. 

(2)  It  is  possible  that  the  washing  procedure  used  during  the  test  might  lead 
to  an  excess  of  particles  at  the  periphery.  As  however  the  number  of 
complete  squares  very  close  to  the  edee  is  very  small,  not  enough  data  was 
collected  to  test  for  this  particular  form  of  maldistribution. 
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5.3  A  Preliminary  Report  on  the  Hiac  Particle  Size  Analyser 
5.3.1  Introduction 


The  laboratory  had  the  use  of  a  Hiac  counter  for  10  days  and  during  that 
time  made  a  brief  examination  of  the  performance  of  the  instrument.  The  tests 
were  carried  out  using  spherical  pollen  particles  42 Pm  dia.  and  the  pulse  shapes 
recorded  on  a  Tektronix  Storage  Oscilliscope  type  564B.  Because  of  the  veiy  short 
time  available  this  report  is  necessarily  extremely  tentative. 

5.3.2  Pulse  Shape 

42  V  m  pollen  was  introduced  to  the  counter  and  the  pulse  shape  presented 
at  the  output  socket  of  the  counter  was  recorded.  Preliminary  investigations 
suggest  that  this  output  is  the  transducer  signal  amplified  and  that  the  pulses  have 
not  been  shaped.  The  pulse  shape  was  not  of  the  expected  flat  topped  type,  i.e. 


Observation  confirmed  the  statement  in  the  instrument  manual  that  the  pulse 
height  was  a  variable  of  particle  residence  time  in  the  counting  head  (figs. .36, 

37  ) .  If  as  has  been  suggested,  particles  spin  whilst  in  the  sensing  zone, 

elongated  particles  would  be  expected  to  produce  approximately  flat-topped  pulses 
with  a  series  of  maxima  on  them,  whilst  spherical  particles  should  produce  a  flat- 
topped  pulse.  The  claim  that  flat-topped  pulses  are  produced  implies  that  the 
rise  time  of  the  dector  is  small  compared  with  the  residence  time  of  the  particle 
within  the  secsor.  Both  the  shape  of  the  pulse  and  the  dependence  of  pulse  height 
upon  particle  flow-rate  would  appear  to  contradict  this. 

5.3.3  Background  Count 

Clean  water  triply  filtered  through  0.2  ym  aperture  membrane  filters  was 
used  to  attempt  to  reduce  the  particle  count  at  the  small  size  end.  Some  reduct¬ 
ion  did  occur  but  the  counter  continued  to  register  counts  including  the  occasional 
particle  in  excess  of  20  ym.  There  are  two  possible  explanations  for  this  (a) 
that  some  particles  remained  from  previous  counts,  (b)  that  the  registered  counts 
were  due  to  noise  in  the  counter.  From  the  work  carried  out  it  was  not  possible 
to  differentiate  between  these  effects  although  the  fact  that  The  counter,  in  spite 
of  making  use  of  integrated  circuitry,  took  2  hours  to  stabilize,  suggested  that  there 
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Pulse  height  vs  flow  rate  of  HIAC  Counter 


might  be  some  electronic  noise  within  the  system. 

5.3.4  :  Pulse  Shapes  Produced  by  Irregular  Particles 

Runs  were  carried  out  using  40b m  sieved  aluminium  flake  particles. 

The  pulses  observed  were  identical  to  those  produced  by  spherical  particles. 

5.3.5  Conclusions 

The  extremely  good  reproducability  obtained  suggests  that  this  instrument 
is  probably  an  extremely  good  comparative  method  of  counting  particles.,  but  its 
value  as  an  absolute  counting  device  has  yet  to  be  confirmed.  More  work  is  needed 
on  analysis  of  pulse  shapes  and  the  effect  of  operating  parameters  of  the  results 
obtained.  The  present  conclusions  assure  that  the  pulses  observed  were  not  shaped 
before  the  output  socket.  This  assumption  needs  to  be  checked  by  experiment  and 
reference  to  circuit  diagrams. 
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5.4 


The  Coulter  Counter 


5.4.1  Introduction 

One  of  the  features  of  new  developments  in  methods  of  particle  size 
analysis  is  the  increasing  importance  of  stream  scanning  methods  whereby  a  particle 
is  sized  and  counted  as  it  passes  through  a  sensing  zone.  This  field  is  domin¬ 
ated  by  the  extremely  successful  Coulter  Counter  (1)  in  which  the  change  of 
resistance  between  two  electrodes  immersed  in  an  electrolyte,  is  monitored  as  a 
particle  passes  through  a  small  orifice  separating  the  electrodes.  Other  develop¬ 
ments  are  in  the  use  of  the  measurement  of  the  light  scattered  from  individual 
particles  as  the  means  to  count  and  identify  particles.  3oth  of  these  types  of 
measurement  have  the  ability  to  count  and  analyse  large  numbers  of  particl.es  and 
so  give  a  statistically  meaningful  analysis.  However,  ho  matter  how  dilute  the 
suspension  used  there  is  always  a  finite  clance  that  two  or  more  particles  will 
arrive  in  the  sbne.'.ng  zone  simultaneously.  The  response  of  the  instrument  to  such 
an  occurrence  is  important  as  it  can  give  rise  to  an  inaccurate  analysis.  It  is 
the  purpose  of  this  paper  to  examine  the  response  of  the  Coulter  Counter  type  of 
instrument  to  two  or  more  particles  within  the  sensing  zone. 

The  earliest  correction  fo?.  coincidence  is  that  developed  from  the  paper 
by  Mattern  et  al  (1)  and  is  the  correction  recommended  by  the  manufacturers. 

The  undercount,  n11  the  number  of  particles  not  counted  is  related  to  the  actual 
count  n1,  the  orifice  diameter  D  and  the  volume  sampled  V 


This  expression  contains  the  realisation  that  a  correction  must  depend  on  the 
orifice  dimensions  and  the  volume  sampled  but  its  origins  are  not  explicit. 

The  second  approach  was  that  by  Wales  and  Wilson  (2)  who  realised  that  the 
response  of  the  instrument  could  give  rise  to  two  types  of  coincidence. 

(a)  All  pulses  from  two  or  more  coincident  particles  add  so  as  to  produce  a 
single  pulse  whose  maximum  is  the  sum  of  the  individual  pulses. 

or  (b)  the  pulses  from  two  or  more  simultaneous  particles  produce  a  single  pulse 
whose  maximum  corresponds  to  the  largest  particle. 


They  cohsidered  the  sensing  of  particles  to  be  equivalent  to  the  successive 
filling  and  emptying  of  the  sensing  volume  v.  The  chance  of  finding  n  particles 
in  such  a  volume  is  obtained  by  the  application  of  F-olsson’s  Law. 


where£is  the  expectation  of  finding  one  particle  ii.  the  volume  v 

^  =  v,  .  No  where  No  is  the  number  of  particles  in 

V 

the  volume  V. 

-€• 

The  probability  of  no  count  is  p(o)  =  £  and  the  probabilitv  of  some 
count  is  therefore  p(n>o)  =  1  -  p(o)  =  1 

They  finally  conclude  that  the  number  of  particles  counted  will  be 

n  =  No  .  (i 

£ 

This  equation  is  in  error  since  if  the  true  number  per  volume  V  is  No  and  n  is  the 

numbsr  counted,  the  probability  of  a  count  is  and  not  — ^  .n  , 

No  Ho  Wales  ana  Wxxson 

continue  to  examine  coincidence  as  a  function  of  response  level  (threshold)  for  each 
of  their  assumptions . 


where 


n(t)  =  Nofc 


+  H  f2  +  4"  f3  + 
6 


where  f^  are  dependent  on  their  assumptions  (a)  or  (b).  A  solution  can  only  be 
obtained  for  assumption  (b)  when  .....  R  (t) 


N  (t)  = 

o 


1  -  J 


n  (t) 


This  paper  was  criticised  by  Kubitschek  (3)  who  stated  that  only  assumption 
(a)  was  valid  and  this  lead  to  a  further  paper  by  Princen  end  Kwlek  (4)  who 
assumed  that  the  'time  distance'  Ti  between  particles  was  distributed  randomly. 

From  the  mean  'time  distance’  %  between  particles,  they  calculated  that  a  fraction 
*  will  be  separated  by  a  distance  Jess  than  t  the  time  soent  in  the  sensing  volume 

2*e 

and  so  they  calculate  the  correction 


n  =  N  -  CN  where  n  is  the  number  ccunted  and  N  is  the  true  number. 

C  is  a  constant  =  v  .  This  relation  has  been  used  experimentally  by  Edmunson  (5). 

2  V 
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It  is  the  purpose  cf  this  paper  to  re-examine  the  assumptions  of  Wales 
and  Wilson  and  to  derive  correct  coincidence  corrections  for  the  Coulter  Counter. 

5, 4.2  Response  of  the  Coulter  Counter 

The  predictions  of  the  resistance  change  in  the  orifice  due  to  the  presence 
of  a  particle  have  oecn  presented  in  papers  by  Batch  (6)  and  by  Gregg  and  Steidley 
(7) .  The  analysis  depends  on  the  particle  being  in  a  region  of  uniform  field  and 
the  approaches  differ  only  in  the  methods  solving  the  integral.  The  solution 
obtained  by  Gregg  and  Steidley  for  the  sphere  is 


This  equation  is  approximated  by: 

-n  =  r  ay  1  +  0.3  K2  -t-  0.13  K4  *  .  .  .  ) 

A2  1  -  K2 

and  illustrates  that  provided  K  is  small  the  resistance  change  is  proportional  to 
the  volume  of  the  particle.  The  difference  between  the  assumption  that  the 

resistance  is  proportional  to  the  volume  and  the  Gregg  and  Steidley  result  has 
been  examined  by  Allen  (8)  and  has  been  examined  experimentally  by  Eckoff  (3) . 

Grover  et  el  (10)  with  an  alternative  approach  apply  Maxwell's  expression 
for  the  resistance  of  a  composite  medium  to  the  case  of  a  particle  in  a  uniform  field 
and  obtain  an  expression  for  the  relative  change  in  current  obtained  for  the  presence 
of  a  particle 


1.5S 
1  -6 


whereS  is  the  ratio  of  particle  volume  to  that  of  the  orifice.  Provided 
10 


t  \  ] 

VR  i 


l  where  li  is  the  radius  of  the  concentric  spheie  in  which  the  field 
3 


is  uniform  in  the  absence  of  the  particle.  This  expression  is  examined  experimentally 
in  a  subsequent  paper  by  Grover  et  al  (11). 
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5.4.3 
(a)  1 


Coulter  Analogue  Results 

The  Response  to  Single  Particles 

Because  of  the  physical  dimensions  of  the  orifice  it  is  difficult  to 
obtain  precisely  the  resistance  change  with  position  in  a  real  situation. 
However,  an  electrically  similar  analogue  can  be  built  in  which  the  phys¬ 
ical  dimensions  are  scaled  up  and  an  electric  field  is  maintained  across 
the  enlarged  orifice.  Large  diameter  particles  can  now  be  physically 
drawn  through  the  orifice  and  the  resistance  change  with  position, monitored. 
The  physical  dimensions  are  illustrated  in  fig.  38  which  also  demonstrates 
the  pulley  system  used  to  pull  the  particle  through  the  orifice.  As  in 
operation  of  the  Coulter  orifice  a  constant  current  was  maintained  across 
the  orifice  and  the  change  in  voltage  required  to  maintain  the  current 
constant  was  recorded  graphically  on  an  x-y  plotter.  The  position  of  the 
particle  was  monitored  by  movement  of  a  spiral  wound  potentiometer  attached 
to  the  pulley  system. 

A  typical  response  of  voltage  against  position  is  illustrated  in  fig.39 
for  spherical  particles.  From  these  curves  the  peak  response  can  be 
obtained  as  a  function  of  particle  volume.  This  variation  is  presented 
in  table  1  and  illustrated  by  fig.40  which  shows  that  the  proportion¬ 
ality  between  peak  response  and  particle  volume  is  good  for  spheres  up  to 
50%  of  the  orifice  diameter,  a  figure  in  excess  of  that  previously  accepted. 

The  results  also  show  that  the  resistance  of  the  orifice  is  changed  by  a 
particle  well  outside  the  physical  limits  of  the  orifice  and  that  the  sens¬ 
itive  region  is  a  function  of  particle  size.  To  investigate  the  reason  for 
this  further, the  orifice  and  the  electrode  system  have  been  represented  on 
graphited  paper.  A  potential  placed  between  the  electrodes  enable  lines  of 
equipotentittl  to  be  obtained.  A  particle  approaching  the  orifice  is  now 
represented  by  the  removal  of  a  disc  of  the  conducting  paper  and  the  lines 
of  equipotential  obtained  once  again.  The  results  are  sketched  in  fig. 

41  and  show  how  the  distribution  of  potential  is  changed  by  the  presence 
of  the  particle.  Whilst  it  is  not  suggested  that  the  electric  field 
obtained  in  this  way  is  the  field  in  the  cylindrical  geometry  of  the  Coulter 
Counter,  it  is  interesting  to  note  the  change  in  fit  id  occurs  and  to  note 
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FIGURE  38  Analog  "oultor  Counter  -  general  diagram 
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FIGURE  40  Typical  response  to  a  spherical  particle 


the  influence  of  the  particle  outside  the  physical  limitations  of  the 
orifice.  These  results  suggest  that  the  sensitive  region  can  be  divided 
into  two  regions  shown  in  fig.  42.  Region  I  is  the  region  of  uniform 
field  in  which  a  particle  will  give  the  peak  resistance  change.  Region  II 
is  a  region  in  which  there  is  a  response  due  to  the  presence  of  the  particle 
but  that  response  will  depend  on  the  position  of  the  particle.  Outside 
region  II,  is  the  electrolyte  in  which  the  presence  of  particles  will  not 
cause  a  resistance  change.  In  fig.  42  the  two  regions  labelled  £1  are 
shown  as  being  symmetrical  but  this  will  only  be  so  if  the  electrodes  are 
placed  symmetrical  to  the  orifice.  These  conclusions  were  shown  to  be  valid 
by  further  experiments  in  the  Coulter  Analogue.  Figure  43  demonstrates 
the  effect  of  doubling  the  orifice  length.  The  pulse  becomes  more  flat 
topped  but  is  otherwise  symmetrical,  thus  indicating  that  the  region  of 
uniform  field  is  proportional  to  the  orifice  length. 

(b)  The  response  to  two  or  more  particles 

When  two  equal  sized  particles  separated  by  a  distance  d  are  drawn  through 
the  orifice,  the  response  is  complex  and  is  illustrated  in  fig.  44.  The 
peak  resistance  change  is  plotted  as  a  function  of  the  separation  in  fig.  45 
These  resistance  changes  can  be  reconstructed  by  adding  together  the  spacial 
response  of  single  particles  with  the  peak  responses  separated  by  the 
distance  between  the  particles  fig.  46.  Further  complex  responses  will 

be  obtained  if  three  or  more  particles  are  found  in  the  sensitive  regions 
at  the  same  time.  Predicted  pulse  shapes  are  shown  in  fig.  47  for  three 

particles  and  figs  48-53  show  oscillographs  of  pulses  obtained  using  pollen 
in  a  Model  A  Coulter  Counter. 

5.4.4  The  Interactions  between  monosized  particles 

The  following  interactions  can  now  be  seen  to  occur. 

1.  A  single  particle  separated  from  all  other  particles  by  a  distance  will 

produce  a  peak  response, A  R,  proportional  to  its  volume  within  the  limit¬ 
ations  discussed  earlier.  Where  the  distance^  is  equal  to  the  distance 
between  the  outer  limits  of  the  two  regions  II  see  fig.  54. 
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Regions  of  field  about  orifice 


FIGURE  45  Psalt  voltngo  response  as  a  function  of  particle  separation 


FIGURE  50 


Two  Particles  Showing  Intermediate  Response 
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;  FIGURE  51 

! 

I 


Two  Particles  Showing  Almost  True  Double t 


FIGURE  52 


2.  Two  or  more  particles  found  within  a  distance, ^defining  region  I  and  with 
all  other  particles  at  a  distance  greater  than  t  away  will  give  a  'vertical 
interaction’  response.  The  response  will  be  proportional  to  the  total 
volume  of  the  particles  within  the  region. 

3.  With  two  particles  at  a  distance  apart  between  Jl  and  &  ,  all  other  part- 

X  £ 

icles  being  at  a  distance  grf  .ter  than  i.  away  ,  the  response  will  be  inter- 
mediate  between  that  for  a  particle  of  unit  volume  and  that  of  double  volume, 
and  will  be  dependant  on  the  precise  separation. 

4.  With  three  or  more  particles  within  a  distance  SL  ,  all  other  particles 
being  greater  than  %  apart,  the  response  is  complex  and  will  depend  on 
the  geometry. 

It  is  interesting  to  note  that  the  horizontal  interaction  of  Wales  and 
Wilson  does  not  occur  but  on  more  complex  response  giving  rise  to  ’apparent'  part¬ 
icles  of  volumes  intermediate  between  that  of  integral  number  of  particles  is 
obtained.  This  is  further  demonstrated  by  the  apparent  size  distribution  obtained 
when  counting  concentrated  dispersions  of  monosizad  particles.  Figure  55  demons¬ 
trates  the  apparent  change  in  size  distribution  of  27  v*pollen  in  a  140  Worifioe, 
where  the  only  difference  is  the  concentration  of  the  particles. 

5.4.5  The  Theory  of  Coincidence 

If  the  suspension  from  which  particles  are  drawn  has  sufficiently  large 
volume  so  that  the  removal  of  some  particles  does  not  materially  change  the  con¬ 
centration  of  particles  and  if  there  are  no  other  physical  conditions  affecting 
the  separation  between  particles,  the  distance  between  particles  will  follow  the 
Poissonian  distribution. 

The  probability  of  finding  a  separation  between  x  and  x  +  d^  is 

<$P(x)  =  p  exp  (  -  vi  x)  dx 

where  p  is  the  reciprocal  mean  separation 

ird2  No 
V 
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Interaction  1  True  single  pulses 


Since  true  pulses  occur  if  there  is  a  separation  greater  than  £.  on 
each  side  of  the  particle  the  probability  of  a  true  singlet  will  be 

-00  -00 

exp  <-  px)  dx 


P(l) 


■JF 


/" 


P(l)  =  exp(  -  2p*  2) 


pexp  (-P  x)  dx 
*2 

.  .  .(1) 


and  the  probability  of  particles  not  being  counted  as  a  single  particle 

=  1-  P(l)  =  1  -  exp  (-2  pi  }  an  approximate  form  of  ecuation  (1) 

2  2 

is  P(l)  --  1  -  2  ltd  N  o  _ 

0  2. 


so 


N=  N  (1-CN)  the  equation  of  Princen  and  Kwolek 
o  o 


Interaction  II  True  multiplet  pulses 
(a)  Doublet 

If  two  particles  are  found  at  a  separation  less  than  ^  and  on  either  side 
of  the  doublet  there  are  separations  greater  than  i  ,  then  the  response  will  be 
that  of  a  doublet. 


P(2) 


.  jr. 


exp  (-  |*x)  dx 


ii 
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exp  (-J#x)  dx  j  exp  (- |/x)  dx 

'2 


P(2)  =  exp  (-  2pi2)  (  1  -  exp  (-  p£  ) 


Triplet 


If  three  particles  are  found  within  the  distance  £  then  the  response  will 


be  a  true  triplet. 


<0 

r 
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p  exp  (-  x)  dx 


J 
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p  exp  (-p  x)  exp  (-p(y-x))  dx  dy/  p  exp(-px)dx 
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P(3) 


=  exp  (-  2p  #2)  (1  -  exp  (-  p  2^)  (1  +  p^)  ) 


Quadruplet 


Similarly 


.  r,.„,  <“■  ’tr. 


P(4)  =.  1  p.  exp  (-  ux)  cix 
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dy  fly.'  exp  (-  px)  exp  (-  pz) 


o  /o  J  o 


exp  (p(y-x-z)  x  dy  dz 


p  exp  (-  px)  dx 


2  m2 

.  P{4)  =  exp  (  -  2p  £2)  Cl  -  exp  (-  p  j^)  (1  +  p  x  —  +  .  .  .) 


Multiplet 

The  probability  of  finding  n  particles  in  length  with  separations 
greater  than  on  either  side  is 


n-2 

P(n)  =  exp  (-  2  p  ?2>  (1  -  exp  (-pJt^)  (1  + 

1 

Probability  of  Inteimedidate  Response 
Interaction  III 

-  i-*c. 


Doublets  with  separation  between  and 
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PX(2)  =  exp  (-  2P^)  (exp  (-PS.  ^  -  exp  (-P& 


and  the  probability  of  a  separation  between  d  and  d  +  dd 


dp  (2)  =  exp  (-2p£2)  p  exp  (-P  d)  dd 


Interaction  IV 

Three  further  types  of  response  are  possible  and  are  illustrated  in  Fig. 
56.  All  other  possibilities  are  allowed  for  in  interactions  previously 

considered.  The  pulse  recorded  will  depend  on  the  mechanism  of  pulse  height 
selection  used  but  it  appears  (12)  that  pulse  (a)  will  be  recorded  as  an  inter¬ 
action  III,  a  doublet  with  intermediate  response,  pulse  (b)  will  also  be  recorded 
as  a  doublet  with  intermediate  response  whilst  (C)  will  be  recorded  as  a  triplet 
with  intermediate  response  between  two  and  three. 


Probability  of  pulse  type  (a)  is 
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Pulses  of  the  shape  3(b)  and  3(c)  will  have  the  same  probability 
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Similar  analysis  will  apply  to  interactions  of  four  or  mors  particles. 

5.4.6  The  Magnitude  of  the  Coincidence  Correction 

This  theory  can  be  used  to  investigate  the  <ragni  tade  of  the  coincidence 
effeci  if  values  can  ne  assignee  co  the  critical  lengths  1  and  1  .  The  former 
lj,  is  *he  lengtl  of  uniform  field  within  the  orifice  ar.d  to  a  first  approximation 
may  be  assumed  to  equal  ..he  length  of  the  orifice.  This  approximation  to  the 
critical  length  will  be  too  large  and  will  sliehtly  over  emphasise  the  occiurence  of 
true  mu Itiplets  but  tho  true  effective  distance  of  u-iform  field  will  depend  on  the 
shape  of  the  orifice  and  the  position  of  the  electrodes  and  is  in  practise  difficult 
to  determine. 

The  total  effectively  sensitive  distance  waj  be  obtained  by  observing  the 
shape  of  the  pulses  oulained  when  monosized  par  tides  are  drawn  through  the 
orifice  as  ir.  the  usual  calibration  procedure.  By  overlaying  the  pulses  obtained 
on  the  display  of  an  oscilloscope,  an  effective  time  of  passing  of  the  particle 
is  obtained.  If  one  now  assumes  tin  t  the  particle  is  travelling  at  the  average 
velocity  of  the  fluid  through  the  orifice  an  estimate  of  the  sensitive  distance 
12  is  obtained.  Table  XIII  represents  the  results  obtained  whan  calibrating 
pollens  of  different  sizes  are  passed  through  a  range  of  orifices.  Particles  of 
diameter  small  compared  to  the  orifice  diameter  give  rise  to  ?  fairly  consistant 
estimate  of  the  sensitive  length.  Large  particles  however  give  much  larger 
estimates  of  the  length.  This  is  due  to  the  complex  paths  of  large  particles 
approaching  the  orifice  and  to  the  incorrect  assumption  that  these  particles  will 
be  travelling  at  the  fluid  velocity.  Values  of  the  sensitive  distance  obtained 
with  the  smaller  particles  have  been  taken  as  more  representative  of  the  above 
thecr  y . 

These  values  for  the  sensitive  distances  have  been  substituted  in  the 
equations  presented  above  and  estimates  of  the  occurrence  of  the  singlets,  inter¬ 
mediate  doublets,  doublets  etc.  made  for  a  series  of  initial  particle  concentrations. 
The  results  are  prevented  in  a  series  of  tables  XIII,  XIV,  XV,  XVI. 

These  tables  illustrate  that  the  probability  of  coincidence  is  not  only  a  function 
of  concentration  of  the  particles  but  is  also  a  function  of  the  orifice  tube  diameter. 
Moreover,  it  illustrates  that  particles  not  counted  as  single  particles  will  ba 
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counted  as  multiple  particles  and  so  that  ihe  distribution  obtained  from  che 
counting  of  concentrated  suspensions  my  troll  t a  in  serious  error.  The  tables 
estimate  the  probable  extent  of  that  error. 

5.4.7  Use  of  a  Multichannel  Pulse  Height  Analyser 

Ifh.Mst  there  exists  a  Coulter  Counter  which  will  divide  a  distribution  into 
15  channels  (Model  T),  this  is  not  yet  available  in  the  tl.fi,  and  so  could  not  be 
evaluated.  Ideally  however,  the  voltage  pulse  obtained  when  a  particle  n asses 
through  the  orifice  across  which  a  constant  current  is  maintained  should  be 
divided  into  a  larger  number  of  channels.  A  convenient  instrument  is  the  L;  tern 
Pulse  Height  Analyser  which  will  analyse  the  pulses  into  100  channels. 

However  ,  the  analyser  5  a_  a  unique  analogue  to  digital  converter  sequence 
which  makes  it  superior  to  other  types  oi  analyser  See  supplement . 

In  collaboration  with  ths  manufacturers  a  suitable  interface  has  been 
designed  which  will  permit  the  instrument  to  accept  the  long  pulses  and  also  to 
permit  th«  voltage  level  across  the  ICO  channels  to  be  adjusted  to  a  convenient 
level . 

The  experimental  arrangement  <3. a  So  far  proved  satisfactory.  The  time  for 
a  complete  analysis  is  short  of  J0-30  sees,  tnd  shows  that  this  instrumental  set-up 
is  ideal  for  rapid  particle  size  analysis.  However,  there  is  a  need  for  the 
comparison  between  particle  size  analyses  obtained  or,  this  instrument  and  those 
obtained  on  the  conventional  Coulter  Counter. 

Because  the  Coulter  principle  measure*  the  pavtiC’e  volume  a  100  channel 
analyser  is  restricted  in  the  range  of  particle  diameters  ,iat  can  be  resolved 
accurately  because  of  the  scale  compression  at  the  lower  size  end.  In  erder  to 
overcome  this,  a  logarithmic  amplifier  preceeded  by  an  impedance  converter  has 
been  placed  between  the  Coulter  orifice  and  the  multichannel  analyser.  This 
arrangement  increases  the  input  impedance  of  the  amplification  chain  to  ''-20fiM 
which  facilitates  the  use  of  non-aqueous  electrolyte  systems  as  are  necessary 
when  analysing  mineral  oil  base  fluijt. 

A  further  logarithmic  stage  would  enable  a  cube  root  of  the  pulse  height 
to  be  analysed  thus  giving  a  linear  scale  in  terms  of  equivalent  spherical  diameter. 


5,4,8  Summary  apd  Conclusions 


The  response  of  a  Coulter  Counter  to  spherical  particles  has  beer,  examined 
and  it  has  been  found  that  the  peak  resistance  change  of  the  orifice  is  a  linear 
function,  oi  particle  volume  up  to  50%  of  the  orifice  diameter.  However,  the 
resistance  begins  tc  change  as  the  particle  approaches  the  orifice.  The  effect  of 
this  resistance  change  as  two  or  more  particles  approach  with  different  separations 
has  been  measured  and  it  has.  been  illustrated  how  complex  pulse  shapes  can  be 
obtained. 

The  change  in  resistance  in  response  to  the  approach  of  particles  has  been 
classified  into  four  main  areas  and  the  probability  of  occurrence  of  each  event  has 
Seen  calculated.  Estimates  of  the  critical  lengths  have  been  obtained  for  a  series 
of  orifice  diameters.  These  resuits  are  presented  in  tabular  form  and  illustrate 
that  the  froblem  of  cofncidence  is  not  only  the  problem  of  lost  counts  as  previously 
considered  but  is  reallj  the  problem  of  the  distortion  of  the  particle  size  dist¬ 
ribution  obtained  due  to  the  additive  responses  that  can  be  obtained.  The  analysis 
has  only  been  applied  here  to  monosizad  particles.  The  analysis  will  be  more  complex 
with  many  sized  particles  as  there  will  be  more  possibilities  of  additive  response. 
However,  the  general  solution  will  still  apply.  The  easiest  way  of  reducing  the 
problems  associated  with  coincidence  is  to  count  with  a  low  concentration  suspension 
and  this  will  be  true  for  multi-sized  dispersions  as  well  as  :nono-sized  dispersions 
and  the  tables  presented  can  oe  used  to  select  acceptable  concentration  levels. 

A  suitable  interface  has  been  designed  to  enable  the  orifice  ’transducer' 
to  be  connected  to  a  pulse  height  analysc-r  and  tests  indicate  that  the  arrangement 
is  suitable.  Further  work  is  proceeding. 
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TABLE  XII 

Analogue  Response  to  Single  Spheres 


Particle 

diameter 

cms 

Particle  diameter 

Orifice  diameter 

Particle 

volume 

ml 

< 

Peak 

response 

mv  1 

Length  of 
response 

1 

1.2 

.105 

0.9 

6.9 

- 

2.0 

.175 

5.1 

5.2 

- 

3.0 

.264 

13.5 

12.4 

- 

3.36 

.295 

20 

17.0 

34 

3.86 

.339 

30 

29.0 

36 

4.24 

.372 

40 

37.0 

39 

4.50 

.395 

50 

44.0 

43 

5.76 

.505 

100 

99.0 

46 

6.60 

.579 

150 

150.0 

46 

7.26 

.637 

200 

210.0 

48 

8.30 

.728 

300 

324.0 

49 

9.10 

.808 

400 

450.0 

50 

9.84 

.863 

500 

620.0 

50 

10.46 

.918 

600 

750.0 

50 

1.  Analog  :e  current  maintained  at  0.5  amp. 
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Length  oi  Sensitive  Regions 


TABLE  XV 

Coincidence  Effects  for  200.0pm  Tube 


No/cc 

1000 

3000 

5000 

10,000  1 

* 

30,000 

50,000 

Singlets 

987.0 

2,884.7 

4,684.0 

8,775.8 

29,275.7 

26,025.3 

Intermediate 

Doublets 

9.2 

82.0 

224.1 

857.2 

6,371.1 

14,276.8 

Doublets 

2.7 

24.3 

65.7 

244.5 

1,647.7 

3,428.6 

Intermediate 

Triplets 

.7 

3.0 

22.5 

438.9 

1,475.2 

Triplets 

- 

.1 

.5 

3.4 

68.8 

236.5 

- 

- 

etc. 

etc. 

etc. 
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TABLE  XVI 


Coincidence  Effects  for  280um  Tube 


f  ! 

No/cc  1000 

3000 

10,000 

30,000 

s 

Singlets 

974.7 

2,778.3 

4,399.4 

7,741.7 

13,920.0 

Intermediate 

Doublets 

19.5 

170.3 

457.8 

■ 

10,194.4 

Doublets 

6.9 

58.4 

528.9 

2,662.3 

Intermediate 

Triplets 

.1 

3.4 

15.0 

102.3 

1,478.4 

Triplets 

- 

.6 

2.7 

18.5 

272.6 

- 

- 

■ 

etc. 

etc. 
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Automated  Particle  Size  Analysis  Using  a  Digital  Computer 


5.5.1  Summary 

A  Research  programme  has  recently  been  started  at  Loughborough  with  the 
object  of  developing  improved  methods  of  automatic  particle  size  analysis. 

This  paper  discusses  preliminary  work  on  an  image  analysing  system  that 
is  based  upon  digital  computation  and  shows  how  this  system  can  be  used  for 
particle  sizing. 

The  images  from  optical,  and  electron,  microscopes  are  suitable  starting 
points  for  particle  size  analysis.  These  images,  however,  contain  ambiguities 
which  may  be  easily  dealt  with  by  human  operators  but  not  by  existing  automatic 
analysers.  The  data  processing  capability  of  the  digital  computer  and  its  memory 
are  important  assets  for  an  image  analysing  systan . 

An  image  scanner  is  described  whose  motion  can  be  controlled  by  digital 
means  and  which  supplier  data  that  completely  defines  the  projected  areas  of 
particles  in  the  field  of  view. 

5.5.2  Introduction 

The  derivation  of  various  particle  parameters  from  measurements  obtained 
by  visual  examination  of  a  sample  through  a  microscope  is  a  well  established  tech¬ 
nique.  The  main  drawback  is  the  tedium  involved  in  maxing  these,  measurements. 

A  great  deal  of  effort  has  therefore  been  expended  by  many  workers  with  the  object 
of  automating  the  measurement  process.  A  number  of  commercial  instruments  have 
been  produced  but  these  have  achieved  varying  degrees  of  success. 

The  need  for  automatic  particle  size  analysis,  is  based  upon  visual  exam¬ 
ination,  is  as  important  as  ever  and  the  present  project  is  an  attempt  to  improve 
upon  existing  techniques. 

A  major  problem  that  has  to  be  overcome  in  an  automatic  system  is  the  need 
to  handle  large  quantities  of  data.  The  basic  weakness  of  existing  equipment  is 
that  the  data  is  derived  and  processed  using  analogue  computer  techniques.  Digital 
data  handling,  if  it  is  used  at  all  is  only  utilised  at  the  end  of  the  system  after 
the  primary  measurements  have  been  made  and  the  format  of  the  data  frozen. 
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Although  analogue  circuits  have  a  very  high  data  processing  rate  they  suffer 
from  the  following  limitations. 

(a)  the  equipment  increases  in  size  as  the  number  of  operations,  in  a  calculation, 
increases 

(b)  only  very  limited  data  storage  is  possible  if  cost  is  to  be  minimised 

(c)  they  are  not  as  reliable  or  as  precise  as  digital  equipment 

(d)  if  the  calculations  are  modified  then  extensive  re-construction  of  the 
equipment  will  be  necessary. 

The  value  of  digital  circuitry,  and  computation,  for  image  processing  has 
already  been  demonstrated  in  systems  handling  cloud  cover  pictures  and  identifying 
man-made  objects  on  aerial  photographs.  Additional  techniques,  relevant  to 
automated  particle  sizing,  have  been  developed  by  researchers  concerned  with  machines 
for  computer  recognition  of  tests  composed  of  handwritten,  or  printed  characters. 
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The  methods  being  developed  at  Loughborough  for  size  analysis  are  based 
upon  knowing  the  X,  Y  coordinates  and  image  density  of  any  point  of  interest 
in  the  image  plane.  This  is  subject  to  a  limit  of  resolution  fixed  by  dividing 
the  image  ,>lane  into  a  matrix  of  square  size  400  x  400.  The  availability  of  such 
data  means  that  the  perimeter  of  every  particle,  in  the  sample  being  analysed  can 
be  precisely  specified.  Using  this  data  any  desired  particle,  or  population  para¬ 
meter  can  be  computed  as  long  as  it  has  been  defined  in  terms  of  projected  area 
lengths . 

Automated  size  analysis  can  be  considered  under  the  following  headings: 

1.  Obtaining  an  image  for  scanning 

2.  The  scanning  equipment 

3.  Data  processing 

5.5.2  (i)  Obtaining  an  image 

The  success,  cr  failure,  of  the  automatic  process  depends  to  a  large  extent 
upon  the  quality  of  image  that  can  be  produced  from  the  original  particles.  Ideally 
the  image  should  consist  of  absolutely  sharp  black  objects  on  a  white  background  or 
white  objects  on  a  black  ground.  This  ideal  ic,  of  course,  rarely  attained. 

In  the  present  work  an  image  is  formed  by  projecting  a  photographic  negative 
of  the  sample  onto  the  measuring  plane  of  the  image  scanner.  The  photographic  film 
is  processed  to  obtain  maximum  contrast  and  the  projected  image  is  20  inches  square. 

5.5.2  (ii)  The  scanning  equipment 

The  image  processing  system  can  be  sub-divided  into  the  blocks  shown  in 

fig. 1. 

5.5.2  (ti)  a  The  scanner 

At  present  the  outward  form  of  the  scanning  equipment  being  used  has  been 
dictateu  by  financial  limitations  but  it  does,  however,  contain  all  the  required 
facilities.  A  diagram  of  the  equipment  is  given  in  fig.  58. 

The  image  scanner  is  basically  a  modified  X-Y  plotting  table,  fitted  with 
a  photocell  in  place  of  a  pen,  linked  to  a  digital  computer.  At  present  the 
computer  is  being  used  as  a  development  tool  since  some  of  its  tasks  would  be 
more  suitably  carried  out  by  special  hardware  in  a  final  system. 
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The  computer  is  programmed  to: 


(a)  position  the  photocell  by  sending  pulses  to  the  motors 

(b)  read  the  analogue  output  of  the  photocell  and  convert  this  into 
a  digital  number. 

(c)  pre-process  the  photocell  readings. 

A  scanning  motion  made  up  of  a  series  of  parallel  straight  lines  that 
cover  the  field  of  view  is  extremely  efficient  for  size  analysis.  The  angular 
orientation  of  the  scan  to  the  particles  can  be  varied  either  by  rotating  the 
slide  or  by  changing  the  ratio  between  the  horizontal  and  vertical  increments  for 
the  stepping  motors. 

A  scan  starts  in  the  bottom  left  hand  corner  of  the  image.  After  each 
reading  of  the  image  density,  the  position  of  the  photocell  is  advanced  by  one 
increment  in  the  scan  direction.  When  the  end  of  a  scan  line  is  reached,  the  Y- 
axis  motor  is  advanced  and  the  X-axis  scan  re-started.  The  advantage  of  using 
incremental  positioning  of  the  photocell  is  that  the  X,  Y  coordinates  are  readily 
available  in  digital  form  for  each  measured  value  of  image  density. 

Interleaved  with  the  program  for  positioning  the  photocell  and  reading  the 
phoxocell  output  is  a  £  ogram  which  carries  out  some  preliminary  processing. 

This  pre-processing  routine  detects  particle  edges  and  measures  the  particle 
interceptions  with  each  scanned  line  of  the  image  field. 

The  simplest  form  of  boundary  detection  is  to  use  a  fixed  reference  density 
to  determine  wj.ether  the  photocell  is  currently  positioned  on  a  particle  or  the 
background.  This  procedure  has  been  used  in  many  automatic  instruments  and  is 
used  in  the  present  project  when  the  image  contrast  is  adequate  and  no  anomalies 
are  present. 

The  reference  level  is  chosen  by  making  a  number  of  exploratory  traverses 
across  the  image,  recording  the  photocell  output  with  a  chart  recorder,  and  deciding 
on  a  suitable  value  after  visual  examination  of  the  results. 

Once  the  photocell  has  traversed  a  particle,  the  X,  Y  coordinates  of  its 
initial  intersection  with  the  line  of  scan  are  recorded  together  with  the  length 
of  the  intersection. 
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Since  each  particle  intersects  an  appropriate  number  of  scan  lines,  and  all 
of  these  intersections  are  recorded,  its  two  dimensional  projected  image  in  the 
plane  of  the  photographic  emulsion  is  fully  defined  in  terms  of  X,Y  coordinates 
and  chord  lengths. 

5.5.2  (ii)  b  The  processor 

Typical  output  data  from  the  preprocessor  are  given  m  table  XVI t 
pictorially  ir.  fig.  59  The  essential  difference  between  the  table  and  the 

figure,  however,  is  that  in  fig.  53  the  re-recording  onto  X,Y  coordinates  r.as 
established  for  a  human  observer  which  chord  belongs  to  which  particle  Tab.-e  S'-  i  x 
gives  the  data  as  it  is  held  in  the  computer  core  store  and  it  is  not  apparent 
to  the  machine  which  chords  constitute  a  particular  particle. 

The  first  routine  in  the  processor  is  therefore  a  sorting  routine  to  find 
out  how  many  particles  there  are  and  which  chords  belong  to  each  particle.  Ihe 
criterion  for  determining  the  association  of  chords  with  particles  is  shown  in 
fig.  60. 

5.5.3  Computation  of  size  parameters  for  individual  particles 

Table  XVII  lists  the  chords  for  the  particle  shown  pictorially  in  Fig.  61 
The  boundary  of  this  particle  is  defined,  in  terms  of  X,Y  coordinates  only,  in 
Table  XIX  Table  XIX  is  computed  by  using  the  length  of  a  chord  and  its  X,Y 
co-ordinates  at  the  left  hand  end  to  obtain  the  right  hand  coordinates. 
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TABLE  XVII  Output  from  preprocessor  for  particles  in  fig.  59 


1 

X 

Y 

L 

X  Y 

_L 

X 

Y 

L 

X _ 

_Y _ 

_ L 

153 

2 

7 

127  25 

34 

‘i 

35 

20 

160 

41 

13 

f  149 

3 

14 

61  26 

13 

61 

35 

i3 

4 

42 

30 

145 

4 

20 

128  26 

34 

120 

35 

36 

122 

42 

29 

f 

j  142 

5 

23 

58  27 

17  = 

161 

35 

1 

162 

42 

U 

;  K2 

6 

24 

125  27 

37 

4 

36 

22 

4 

43 

30 

1  140 

n 

» 

27 

59  28 

17 

63 

36 

11 

122 

43 

28 

1  141 

8 

27 

126  28 

36 

121 

36 

36 

161 

43 

12 

i  \3l> 

9 

25 

58  29 

19 

161 

36 

10 

5 

44 

32 

140 

10 

30 

124  29 

36 

3 

37 

24 

123 

44 

27 

:  138 

11 

32 

12  30 

2 

65 

37 

9 

163 

44 

10 

J 

138 

12 

32 

60  30 

18 

120 

37 

34 

5 

45 

31 

135 

13 

34 

124  30 

37 

160 

37 

13 

122 

45 

26 

134 

14 

35 

7  31 

12 

5 

38 

24 

163 

45 

10 

133 

15 

35 

60  31 

17 

65 

38 

8 

7 

46 

31 

133 

16 

35 

122  31 

3^ 

122 

38 

32 

123 

46 

25 

131 

17 

35 

7  32 

14 

160 

38 

IS 

165 

46 

8 

132 

18 

34 

61  32 

17 

4 

39 

25 

7 

47 

31 

129 

19 

35 

123  32 

36 

65 

39 

5 

123 

47 

24 

131 

20 

35 

5  33 

17 

121 

33 

31 

165 

47 

8 

129 

21 

34 

60  33 

16 

159 

39 

14 

8 

48 

31 

120 

22 

33 

121  33 

37 

5 

40 

26 

124 

48 

24 

62 

23 

3 

l63  33 

4 

68 

40 

2 

186 

48 

7 

12S 

23 

33 

5  34 

19 

122 

40 

31 

9 

49 

29 

63 

24 

7 

62  34 

15 

161 

40 

12 

124 

49 

3 

130 

24 

32 

12.1  34 

37 

4 

41 

28 

131 

49 

15 

61 

25 

10 

163  34 

6 

122 

41 

28 

168 

49 

3 

163. 
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J 
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These  chords  belong  to  separate  particles  since  x  does  not  equal  x  and 


in  fact  is  less  thaw  x„.  Also  £.  is  not  long  enough 

2  i 
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2' 


(x2y2) 
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FIGURE  60  Associating  chords  to  fora  particles 
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TABLE  XVIII 


Chords  for  particles  shown  in  Fig.  61. 


(increments  of  0.1  "  for  X  and  Y) 


X 

Y 

_ L 

X 

Y 

L 

77 

1 

7 

63 

17 

18 

74 

2 

11 

63 

18 

17 

73 

3 

13 

63 

19 

16 

73 

4 

14 

64 

20 

15 

72 

5 

16 

63 

21 

14 

71 

6 

16 

64 

22 

13 

69 

7 

18 

63 

23 

13 

68 

8 

18 

65 

24 

11 

67 

9 

17 

69 

25 

6 

67 

10 

17 

73 

26 

1 

67 

11 

16 

66 

12 

17 

65 

13 

18 

65 

14 

18 

63 

15 

19 

63 

16 

19 
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.CURE  61  Chords 


TABLE  XIX 


The  boundary  of  particles  (Fig.  61)  in  terras  of 
x,  y  co-ordinates 


X 

_ Y 

X 

Y 

X _ 

_ Y 

X 

_ Y 

X 

77 

1 

65 

14 

74 

26 

83 

13 

83 

74 

2 

63 

15 

75 

25 

83 

12 

82 

73 

3 

63 

16 

76 

24 

83 

11 

81 

73 

4 

63 

17 

76 

23 

84 

10 

80 

72 

5 

63 

18 

77 

22 

84 

9 

79 

71 

6 

63 

19 

77 

21 

86 

8 

78 

69 

7 

64 

20 

79 

20 

87 

7 

68 

8 

63 

21 

79 

19 

87 

6 

67 

9 

64 

22 

80 

18 

88 

5 

67 

10 

63 

23 

81 

17 

87 

4 

67 

11 

65 

24 

82 

16 

86 

3 

66 

12 

69 

25 

82 

15 

85 

2 

65 

13 

73 

26 

83 

14 

84 

1 
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Once  a  particle  boundary  has  been  defined  by  a  set  of  X,  Y  coordinates  it  is 
possible  to  compute  any  defined  parameter.  Such  parameters  might  be 

3.1.1  The  diameter  of  the  circle  of  equal  projected  area 

3.1.2  The  diameter  of  a  circle  having  the  same  projected  perimeter 

3.1.3  The  Feret  diameter  for  any  viewing  angle  in  the  2-D  plane 

3.1.4  Martin’s  diameter 

(a)  Circle  of  equal  projected  area 

The  projected  area  of  the  particle  is  given  by 

N 


A  z 


where  £  =  length  of  the  n  th  chord  in  terms  of  x-axis  increments 

n 

N  =  the  total  number  of  chords  for  the  particle  when  scanned  using  x 
and  y  increments  of  size  I 

The  circle  of  same  projected  area  is  given  by 


(b)  Circle  of  equal  projected  perimeter 

The  perimeter  of  a  particle  P  (fig.  62)  can  be  found  using  the  set  of 
X,  Y  coordinates  defining  the  particle  edge.  The  circle  of  equal  projected 
perimeter  is  then  given 


(c)  Feret  diameters 

Feret’s  diameter  can  be  readily  calculated  for  a  viewing  axis  which  is 
parallel  to  either  the  x  or  y  axis.  Fig.  63  illustrates  these  situations.  If  the 
viewing  axis  is  parallel  to  the  y-axis  then  Feret’s  diameter  is  given  by 


(x, 


xl) 
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point  x,  y,  on  the  particle  perimeter  is  the  point  with  the  smallest  x-co- 


g 

i 

I 

m 


ordinate  and  the  point  x y 2  is  the  point  with  the  maximum  x-coordinate. 

When  the  viewing  axis  is  parallel  to  the  x-axis  the  points  are  chosen  on 
the  basis  of  minimum  and  maximum  y-coordinates . 

(d)  Martin’s  diameter 

The  area  of  a  particle  in  terms  of  increments  is  given  by 
N 

A  =  21  Z 

,  n 
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Therefore  to  find  the  Martin’s  diameter  parallel  to  the  x-axis  it  is  only  necessary 
to  find  the  r  th  chord  where 

r 


The  computation  of  Martin’s  diameter  when  the  viewing  axis  is  parallel  to  the 
x-axis  requires  slightly  more  work.  A  set  of  chords  for  a  scan  direction  parallel 
to  the  y-axis  must  be  computed  using  the  initial  set  of  chords.  After  this  the 
calculation  is  the  same  as  before. 

5.5.4  Population  parameters 

The  commonly  accepted  frequency  distributions  for  the  scanned  sample  can  be 
computed  from  the  results  for  individual  particles.  Since  the  total  number  of 
particles  on  the  slide  has  already  been  computed  it  is  straightforward  calculation 
to  obtain  t’  a  mean  value  and  standard  deviation. 

5.5.5  Conclusion 

At  present  the  procedure  for  size  analysis  using  the  automated  procedure  is:- 

(a)  photograph  the  sample  using  either  a  microscope  or  direct  photography 
with  a  camera. 

(b)  project  the  negative  and,  using  the  data  from  the  scanner  and  pre¬ 
processor  ^nerate  che  type  of  drawing  shown  in  fig.  59.  This  is 

done  automatically  with  an  incremental  graph  plotter. 

(c)  use  the  main  processing  programme  to  compute  the  relevant  parameters. 

Our  entire  system  obviously  requires  much  effort  before  we  can  feel  satis¬ 
fied  with  it.  Apart  from  changes  in  the  outward  form  of  the  hardware  the  soft¬ 
ware  must  be  improved.  Some  particular  areas  of  interest  are:- 

(a)  touching  particles 

C b)  filtering  out  the  effect  of  dust  and  scratches 

(c)  particles  surrounded  Ly  a  halo  or  containing  areas  of  the  same  contrast 
as  the  background. 

(d)  better  methods  for  edge  detection 

(e)  on-stream  three  dimensional  analysis 
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6.1 


Introduction 


Filtration  is  a  means  of  separating  solids  from  liquids  and  is  mad<=-  use  of 
in  hydraulic  systems  to  remove  contaminant  particles  from  the  hydraulic  oil.  Many 
parameters  are  involved  in  this  operation  and  the  design  of  filtration  equipment  is 
an  attempt  to  optimise  these  for  a  particular  situation.  This  chapter  consists  of 
two  sections,  the  first  using  the  technique  described  in  chapter  4  to  describe  the 
flow  of  liquid  through  a  porous  bed  which  could  be  either  a  filter  cake,  filter 
medium  or  combination  of  the  two.  The  second  section  proposes  a  retention  efficiency 
curve  for  characterising  the  efficiency  of  a  filter  and  describes  how  this  curve  may 
be  obtained. 

The  pressure  distribution  in  a  bed  of  particles  is  complex  and  is  amenable 
to  control  for  two  main  reasons.  On  a  particulate  scale,  the  force  is  handed  on 
from  particle  to  particle  through  the  point  contacts.  An  individual  particle  may 
deform  elastically,  plastically  or  it  may  fracture.  Its  behaviour  will  depend  upon 
its  own  si zo  and  geometry,  its  own  physical  properties  and  those  of  its  surrounding 
neighbours.  Thus,  the  pressure  distribution  in  a  bed  can  be  greatly  controlled  by 
changing  the  form  of  the  original  particles  or  alternatively  mixing  in  other 
particles.  This  is  the  process  which  is  commonly  known  as  pretreatment. 

On  a  moire  macroscopic  scale,  the  pressure  distribution  is  completely  con¬ 
trolled  by  the  geometry  of  the  equipment.  Whenever  there  is  a  solid  boundary,  the 
pressure  distribution  is  redistributed  due  to  friction  against  the  boundary  and  there 
is  here  great  scope  for  more  sophisticated  design  of  filtration  equipment  where 
traditionally  almost  all  of  the  mechanical  pressure  is  allowed  to  be  transmitted  to 
the  septum. 

The  movement  of  particle  beds  would  seem  to  be  important  in  the  design  of 
filtration  equipment  also  for  two  reasons.  Firstly,  the  pressure  distribution  in  a 
static  bed  of  particles  is  immediately  changed  when  the  bed  starts  to  flow.  More 
important,  the  residual  moisture  in  a  dried  cake  is  held  in  the  form  of  pendular 
moisture  at  the  point  contacts.  It  is  clear  that  if  this  moisture  is  to  be  further 
reduced  substantially,  then  the  particles  must  move  relative  to  each  other  and  thus 
release  this  moisture.  Hie  relative  movement  can  be  caused  by  a  variety  of  forces, 
for  example  body  forces,  mechanical  shear  forces  or  vibration. 
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The  third  mechanism  which  would  seem  to  be  of  interest  to  filtration 
engineers  is  that  of  explaining  the  complex  fluid  flow  patterns  in  a  separation 
system.  For  example,  the  fluid  flow  pattern  in  a  deep  bed  filter  must  clearly  be 
contrived  to  bring  each  particle  into  contact  with  a  solid  surface.  Thus,  a  tortuous 
bed  is  required  rather  than  straight  channels  and  a  precise  description  of  the  fluid,  flow 
pattern  in  terms  of  the  pore  distribution  is  required.  This  description  is  also 
necessary  for  those  types  of  separator,  such  as  centrifugals,  in  which  the  particle 
trajectories  are  calculated  on  the  assumption  that  the  pressure  of  the  particles 
does  not  disturb  the  fluid  flow  pattern.  It  is  known  that  at  very  modest  concentra¬ 
tions  that  the  pattern  is  disturbed  and  that  return  flow  occurs  but  this  cannot  yet 
be  predicted  quantitatively. 

In  the  following  section  an  example  is  given  of  the  type  of  statistical 
solution  which  can  be  made  of  fluid  flow  patterns.  In  this  type  of  solution,  it  is 
not  necessary  to  be  able  to  calculate  the  velocity  and  pressure  at  each  point  in  a 
system.  Quite  simple  geometries  pose  major  mathematical  problems  in  order  to  obtain 
a  classic  solution  of  this  kind.  However,  a  statistical  solution  deals  only  with  a 
linow ledge  of  the  distribution  curve  of  the  variables  and  in  many  cases  only  the  mean 
values.  This  approach  has  been  accepted  for  many  years  in  describing  the  properties 
of  a  gas  by  the  kinetic  theory.  However,  in  that  case  the  geometry  of  the  container 
is  not  important.  Starting  with  a  statistical  characterisation  of  the  geometry  such 
as  proposed  in  4.1  this  can  be  effected. 


6.2 


Viscous  Flow  of  Fluid  in  an  Irregular  Pipe 


6.2.1.  This  technique  of  characterising  a  void  space  will  first  be  applied  to  a 
pipe  which  is  irregular  in  cross-section  but  which  is  straight.  In  this  case,  the 
fluid  velocity  and  the  pressure  gradient  lie  in  the  same  direction  at  all  points  and 
are  both  parallel  to  the  axis  of  the  tube.  It  is  only  necessary  to  consider  the 
effect  of  the  shape  of  the  cross-section  on  the  flew  pattern  and  the  filaments  vill 
only  be  drawn  within  the  plane  of  this  cross-section.  The  filaments  will  be  drawn 
in  every  possible  direction  within  this  plane  but  the  distribution  will  be  normal¬ 
ised  so  that  first  moment  of  the  distribution  is  proportional  to  the  cross- 
sectional  area.  A,  of  the  pipe.  Thus,  for  all  the  filaments  lying  in  any  one 

J  d.J  =  A  .  6.2.1 

thickness  of  the  filaments 
length  of  filament 
number  distribution  of  filaments 


direction 


where 


At 

At 

J 


»(J) 


When  the  distribution  of  filaments  is  plotted  for  all  possible  directions  it 
will  be  normalised  back  to  conform  with  equation  6.2.1. 


The  pressure  gradient  is  equal  at  all  points  in  the  cross-section  and  is 
given  by  the  equation: 


1 

V 


dp 

dz 


=  d  u  +  d  u 

2  2 
dx  dy 


6.2.2 


Consider  a  filament  AB,  lying  in  the  cross-section  of  the  channel  as  shown 
in  Fig,  64.  The  velocity,  u,  of  the  fluid  varies  with  position  along  the  length 
of  this  filament  and  may  be  plotted  as  a  velocity  profile  stretching  between  points 
A  and  B.  It  is  tacitly  understood  that  the  velocity  is  perpendicular  to  the  cross- 
section  and,  therefore,  may  be  considered  to  be  a  scalar  quantity  within  the  plane 
of  the  cross-secxion.  The  first  and  second  gradients  of  the  velocity  are,  however, 
parameters  which  vary  along  the  line  AB  and  which  determine  the  shape  of  the  velocity 
profile  between  AB.  Considering  now  the  fluid  to  be  flowing  through  all  the  possible 
filaments,  any  element  of  fluid  is  passing  through  an  area  in  which  there  is  a 
statistical  distribution  of  three  parameters  which  are  velocity  and  the  first  and 
second  gradients  of  that  velocity.  For  convenience  we  shall  denote -2H-  by  u  and — Jiby  u.  It 
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is  emphasised  that  this  in  no  way  implies  a  variation  with  time  but,  in  this 
context,  is  used  to  denote  the  gradients  of  a  velocity  profile  which  is  itself  a 
section  of  the  total  velocity  profile.  Each  of  the  three  parameters  has  a  distri¬ 
bution  function  which  expresses  the  probability  of  finding  a  particular  value  of  the 
parameter  at  some  point  on  one  of  the  filaments.  These  distributions  are  illustrated 
in  Fig.  65. 

The  velocity  u  will  at  all  points  be  positive,  the  distribution  of  the  first 
gradient  will  be  symmetrical  about  the  origin  and  the  second  gradient  will  be  at  all 
points  negative  if  they  are  measured  in  a  consistent  direction. 

The  total  area  under  each  distribution  function  must  be  the  same  since  It  is 
equal  to  the  total  length  of  line  contained  in  this  area. 


Thus:- 

f  P 

(u)  du 

-r  * 

(u)  du  = 

P  r  (il)  d 

=  A 

J  all 

A 

Jaii 

A 

J  All  A 

At 

Where  At  -  the  thickness  of  one  filament. 

We  shall  show  that  since  the  t!  ree  fields  must  fit  together  into  the  same 
space,  it  is  possible  to  relate  the  mean  velocity  to  the  mean  value  of  the  second 
differential  by  relating  each  to  the  first  differential.  Thus,  the  relationship 
between  the  flowrate  and  pressure  drop  in  the  channel  will  be  deduced,  since  the 
second  differential  is  related  to  the  pressure  gradient  by  (6.2.2.).  Since  the 
filaments  traverse  the  area  in  every  possible  direction  then  the  x  -  y  co-ordinates* 
are  no  longer  relevent,  and  we  shall  consider  the  distance  along  a  filament  from 
some  reference  point  to  be  the  only  co-ordinate.  This  distance  will  be  denoted  by 
x. 

When  all  the  possible  filaments  are  drawn  in  the  cross-section  each  point 
is  traversed  once  and  only  once  in  each  direction.  Since  the  pressure  gradient, 
dp 

is  also  a  constant  at  every  point,  then  equation  3  reduces  to*- 
dz 

2  ~2~ 

d  u  d  u  =  -  1_  JL  dp  =  (t  ........  6.2.3. 

dx^  dy^  2  u  dz 

where  the  bar  denotes  the  average  value  of  the  second  gradient  taken  over  all 
directions.  This  relationship  applies  to  any  smaller  area  drawn  within  the  cross- 
section  which  is  t ran versed  by  filaments  in  every  direction. 
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In  order  to  choose  the  reference  point,  consider  again  the  velocity  profile 
shown  in  Fig.  64  stretching  between  two  points  AB.  We  then  know  two  facts. 
Firstly  that  the  velocity  is  zero  at  points  A  and  B.  Secondly,  that  there  must  be 
some  intermediate  point  at  which  the  first  gradient  is  zero  and  hence  the  velocity 
is  a  maximum.  This  point  will  be  denoted  D  and  will  be  taken  as  the  reference  point 
for  this  filament.  The  reason  for  this  is  tha-d  the  second  gradient  can  be  integra¬ 
ted  along  the  filament  from  this  point  to  any  other  point  to  give  the  first  gradient. 

The  velocity  is  a  maximum  at  this  point  only  for  this  particular  filament. 
For  the  other  filaments  which  traverse  the  same  point  xn  a  different  direction  it 
will  not  be  the  maximum  point  except  for  the  special  point  which  is  the  point  of 
maximum  velocity  for  the  whole  profile.  In  this  derivation,  only  channels  which 
have  a  single  point  of  maximum  velocity,  C,  are  considered.  Channels  which  are  very 
complex  in  snaps  may  have  several  such  points  and  the  implication  of  this  is  that 
some  of  the  filaments  will  have  several  maxima  in  their  velocity  profile. 

The  mean  value  of  the  second  gradient  which  is  proportional  to  the  pressure 
gradient  by  equation  6.2.3.,  will  now  be  related  to  the  mean  velocity.  Consider  the 
point  E  chi  the  line  AB  which  is  distance  j  from  the  reference  point  D. 


Then:-  -uj 


6.2.4 


where  ftj  is  the  first  gradient  of  the  profile  at  point  E.  Considering  still 
Fig.  64  draw  a  line  of  constant  velocity,  u,  which  encloses  a  portion  of  the 
cross-section  and  allow  this  area  to  be  traversed  by  all  the  possible  filaments  in 
every  direction.  Then  the  relationship  given  in  (6.2.4)  applies  to  each  of  these 
filaments.  However,  each  filament  is  divided  into  two  parts  by  a  point  corresponding 
to  point  C  and  the  sectioned  filaments  will  have  a  number  distribution  function  which 
can  be  denoted  by  n,  (j),  j  .  Then  summing  equation  over  all  the  section 
fi laments :- 


—  In  (j)  i  ^  dj 

J  all  j 


J, 

I, 


n  (j)  ft  .  dj 
3 


n  (j)  j  dj 
all  j 


Jn  ( j)d  j  ^  f*  u  dx 
a>1  j  Jo 

Jn  (j) d  j  .  j 

fill?  ^ 
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n  (j)j  dj 
all  j 
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The  right  hand  side  of  this  equation  is  simply  the  second  gradient  integrated 
along  all  the  filaments  divided  by  the  total  length  of  filament. 


J 

J. 

J 


n  (j)  u .  d j 
0 


n  (j)  j  dj 
all  j 


n  (j)  dj 


-  all  j 


I 


n  (j)  dj 
all  j 


u  =  constant 


S  f  n  (j)  j  dj 

J  all  j _  U  "j 

f  n  (j)  dj 

J  all  j 


The  term  j  is  the  average  length  of  the  sectional  filaments. 


The  left  hand  side  of  this  equation  is  the  average  value  of  the  first 
gradient  at  the  ends  of  tne  filaments. 


u .  =  u  j  . .  6.2.5 

J 

This  value  of  the  first  gradient  is  the  average  with  respect  to  number  at  the 
end  of  the  filaments.  However,  it  is  not  the  average  with  respect  to  the  perimeter 
of  the  line  of  constant  u,  since  the  various  filaments  cut  the  perimeter  in  different 
directions.  Consider  an  element  of  the  perimeter,  ds,  as  shown  in  Fig.  66. 

Suppose  that  the  first  gradient  normal  to  this  element  is  given  by 

™  =  u  which  can  be  considered  to  be  constant  for  the  element  ds. 
dr  r 

Then,  in  direction  0  , 
u  -  cos  6 

and  the  number  of  filaments  which  cut  ds  at  this  angle  is  given  by:- 

dn  =  —  cos  6 
fit 
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Thus  the  number  mean  of  u  taken  over  all  directions  is  given  by:- 


“N 


fA 

I  2udnd© 

J  o 

rA 

I  2dnd0 


urds 

t 

ds 

t 


fA  2 

I  2 cos  6d9 


J« 


fA 

I  2cos9d0 

J° 


The  mean  of  u  with  respect  to  perimeter  is  given  by:- 

f4  r 

j  2udsd©  I  2cos9d9 

Jo  ds  ur  J  o 


if 


A 

2dsd 0 
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u  = 


8u 

2NTC 


6.2.6. 


Substituting  equation  (6.2.5) 


-u 

s 


8 

2ii  J 


6.2.7. 


This  relationship  applies  to  each  element  of  the  perimeter,  ds,  and  is  therefore 
true  for  the  whole  perimeter. 

In  order  to  deduce  the  flow  through  the  channel  it  would  be  possible  to 
integrate  this  first  gradient  of  velocity  twice  more  to  give  velocity  and  then  flow. 
However,  it  is  easier  to  integrate  the  volume  contained  under  the  velocity  profile 
with  respect  to  velocity  rather  than  with  respect  to  area.  Thus,  if  the  velocity  has 
a  value  u  when  the  area  contained  within  the  line  of  constant  us  is  a,  then:- 


Q 

or  Q 


6.2.8a 

6.2.8b 
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This  is  illustrated  in  Fig.  67  and  we  shall  take  the  (6.2.8b)  as  being  the 
easiei*  to  evaluate. 


The  variables  involved  are  now  superficially  a  and  u.  Both  these  variables 
are  functions  of  the  number,  n,  and  the  mean  length  j  of  the  filaments  cutting  the 
area  contained  by  the  line  of  constant  velocity.  From  (6.2.1)  it  is  clear  that:- 


nj 


a 

At 


6.2.9 


Now  consider  the  surface  of  constant  u  and  one  of  (u  -  du)  which  enclose  an 
area  da  as  shown  in  Fig.  67  The  mean  value  of  the  first  gradient  of  velocity  in 

this  area  da  is,  in  the  limit,  the  same  as  the  mean  with  respect  to  the  perimeter. 
This  is  illustrated  in  Fig.  68  if  it  is  realised  that  the  total  length  of  line 
.traversing  the  element  ds  dr  is  the  same  in  every  direction. 


If  all  lines  traversing  the  element  da  are  taken  then  the  velocity  difference 
between  the  ends  of  each  is  -du  except  for  the  few  lines  which  have  both  ends 
terminating  on  the  (u-du)  line.  In  the  limit  these  lines  disappear  and:- 


Ju  dx  = 
5 

all  n 


u  da 
s — 
At 


-ndu 


6.2.10 


Differentiating  equation  6.2.9 

da  =  At  (ndj  +  jdn)  6.2.11 

Combining  (6.2.10)  and  (6.2.11) 

u 

-du  =  —  (ndj  +  jdn)  6.2.12 

n 


(6.2.7)  and  (6.2.12)  may  be  combined  to  give:- 


2 

Ju  =  —  _  (jdj  +  dn)  6.2.13 

2  u  n 

71 

Combining  equations  (6.2.8b),  (6.2.9)  and  (6.2.13) 


dQ  = 


2  3 

-«At  (nj  dj  +  j  dn) 
u 


6.2.14 


Mathematically,  n  and  j  are  independent  variables  but  physically  they  will 
combine  to  produce  the  maximum  value  of  Q  which  can  satisfy  (6.2.14).  The  condition 
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for  maximum  Q  is  seen  to  be  that:- 


n  =  N  dJ 

1  3 


Substituting  (6,2.15)  and  (6.2,16)  into  (6.2.14) 
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t! 


At 


N 
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V  O 


3 

j  dj 


4  U  £t  MJ 

v.  2 


6.2.16 


Equation  (6.2.3)  will  be  substituted  in  (6.2.17),  remembering  the  total  area 
of  the  pipe  will  be  substituted  by  tie  relationship:- 

N  JAt«  A 


.*.  O  =  A  i  1 o  — 2  6.2.18 

V  dz  *  2 

There  is  one  final  substitution  to  make.  J  is  not  a  complete  filament. 

Each  complete  filament,  9,  has  been  divided  into  two  parts  to  produce  J.  Therefore 
the  mean  filament  9  has  twice  the  value  of  U 

D  -  2J 

Substituting  into  equation  (6.2.18) 

2 

u  =  — I  i.  D  6.2.19 

2u2  y  dz 

This  equation  is  a  general  relationship  for  a  pipe  of  any  shape  within  the 
limitation  that  it  has  only  a  single  point  of  maximum  velocity.  It  is  easily  shown, 
for  example,  that  for  a  circular  pipe,  the  mean  filament  is  given  by:- 

D  *  «r 

2 

where  r  is  the  radius  of  the  pipe. 

(6.2.19)  then  degenerates  to  the  Hagen- Poiseui lie  equation. 
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6.3 


Filter  Testing 


6.3.1.  Introduction 

In  the  past,  the  testing  of  filters  has  amounted  to  little  more  than  the 
determination  of  the  fractional  amount  of  the  feed  solids  which  has  been  retained  by 
the  filter  with  perhaps  some  attempt  to  estimate  the  maxium  size  of  particle  which 
would  be  found  downstream  of  the  filter.  This  type  of  testing  has  involved  close 
definition  of  the  conditions  of  test  and,  in  particular,  specifications  of  the  test 
solids  by  nature  and  size  distribution.  The  results  of  these  tests  are  useful  for 
the  comparison  of  filters  only  when  the  filters  are  used  under  similar  conditions. 

For  practical  use  the  results  have  only  indicative  value  since  the  conditons  of  use 
and  the  solids  to  be  filtered  will  in  general  differ  widely  from  those  of  the  test. 
More  informative  data  of  the  filter  performance  can  be  obtained  by  determining  the 
actual  point  efficiency,  that  is,  the  fraction  of  each  size  of  particle  which  is 
retained.  This  can  be  expressed  as  a  retention  efficiency  curve  (point  efficiency 
plotted  against  particle  size) ,  the  parameters  of  which  can  be  investigated  as 
functions  of  test  variables. 

6.3.2.  Determination  of  the  Retention  Efficiency  Curve 

Given  the  mass  fraction  of  the  solids  retained  by  the  filter  and  the  size 
distributions  (usually  cumulative  fractional  weight  oversize  or  undersize)  of  the 
feed  and  the  solids  passing  or  retained,  the  two  commonest  methods  ^f  determing  the 
retention  curve  are  either  to  construct  a  histogram  of  retention  efficiency  by 
considering  various  size  ranges  or  to  calculate  the  point  efficiencies  from  the 
mass-frequency  distributions.  Both  methods  are  straight  forward  but  rely  on  having 
accurate  size  distributions,  particularly  at  the  ends  of  the  size  distributions.  In 
the  histogram  method  rather  wide  size  ranges  have  to  be  considered  at  the  extremes  of 
the  ranges  in  order  to  obtain  accuracy  in  the  calculated  mean  efficiencies,  and  in 
the  direct  determination  method  it  is  more  difficult  to  determine  the  values  of  the 
mass-frequency  distributions  at  extreme  sizes  so  that  the  calculated  point 
efficiencies  are  in  greater  error. 

A  suggested  alternative  method  is  to  use  "mono  sized”  fractions  of  powders 
so  that  only  the  fractional  weight  retained  has  to  be  determined  and  the  point 
efficiency  is  found  directlyt  However,  objections  to  this  method  arise  from  the 
difficulty  of  obtaining  fractions  of  small  size  range  in  quantity,  that  the  system 
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Let  D 


Particle  size 


Q(D) 

R(D) 

F(D) 

f 


* 

<!> 


Fractional  weight  less  than  D  in  Feed. 

Fractional  weight  less  tnan  D  in  total  retained  solids. 

Fractional  weight  less  than  D  in  filtrate. 

Fraction  of  feed  passing  unfiltered  to  filtrate. 

Fraction  of  feed  retained  by  all  mechanisms  other  than  retention  at 
the  medium. 

Total  weight  of  retained  solids 
Total  weight  of  feed  solids 

Weight  of  solids  retained  by  the  filtering  action  of  medium 
Weight  of  feed  solids  presented  to  medium 


nCD)  =  Overall  point  retention  efficiency 


nr(D)  =  Point  retention  efficiency  due  to  filtering  action  of  medium. 

The  distribution  of  the  solids  will  then  be  shown  in  Fig.  70  the  quantities 
Q(D),  F(D) , R(D),ij>  being  experimentally  determined. 


For  any  size  D,  we  have  the  mass  balance: 
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is  highly  idealised,  and  any  effect  due  to  a  wide  size  distribution  of  the  feed 
solids  is  eliminated. 

4 

For  these  reasons,  therefore,  the  method  which  is  proposed  is  an  adaption  of 
a  technique  recently  developed  in  these  laboratories  for  the  analysis  of  the 
behaviour  of  classifiers.  Essentially  tnis  consists  of  plotting  the  size 
distributions  of  the  material  retained  or  passing  the  filter  against  the  feed  size 
distribution,  and  analysing  the  curve  obtained,  as  explained  below.  A  further 
advantage  of  the  method  is  that  it  allows  considerably  more  information  about  events 
occurring  inside  the  filter  to  be  obtained  than  by  the  previously  described  methods 
which  present  only  an  overall  picture. 

6.3.3.  Theory  of  Method 

In  the  simple  case,  any  two  of  the  three  size  distributions  which  are 
obtained  from  the  test  are  related  by  the  retention  efficiency  curve.  In  Fig.  69 
the  diagonal  OA  would  be  obtained  if  no  filtration  occurred;  if  the  filtration  were 
such  that  all  material  greater  than  some  size  given  by  (l-.$  )  =  Q(D)  were  retained, 
then  the  lines  OB  AC  would  be  obtained.  0  then  equals  the  fraction; 

n^ss  of  solids  retained  for  cumulative  undersize  distributions, 
mass  of  feed  solids 

In  practice,  curves  such  as  OFA,  OGA,  will  be  found. 

For  the  general  case,  we  must  also  consider  the  situations  when ; 

(a)  the  filter  is  faulty,  that  is,  the  medium  has  broken  down  at  some  point  or 
by-passing  occurred  so  that  a  fraction  of  solids  with  the  feed  size  distribution 
passes  with  the  filtrate; 

(b)  A  fraction  of  the  feed  solids  is  held  up,  either  due  to  the  design  of  the 
system  or  by  the  action  of  the  cake  which  has  formed.  ' ..  effect,  this  material  has 
not  been  filtered  by  the  medium. 


Sufficient  data  are  therefore  available  to  determine  f ,  r,  41  , .  However,  as  $ 
can  be  experimentally  determined,  a  check  of  f  and  r  is  avail  le. 

Fig.  71.  illustrates  a  plot  for  the  general  case. 

A  further  check  of  the  tangents  is  available  by  consider.  ;  their  intersections, 
H,  K.  on  Fig.  71  It  is  easily  shown  that 


5h  -  \ 

We  have  also,  that 

=  »!> 


=  1  ~  -  f 

1  -  f  -  r 


(1-f-r) 


6.3.7 


i  •  i 


and 


f 

I-*, 


1-f-r 

1-fr-f 

1-f-r 


Si  Sc 


6.3.8 

6.3.9 


To  determine  the  point  efficiency  due  to  the  medium,  n  ,  equation  (6.3.1)  or 

T 

JR  dF 

(6.3.2)  can  now  be  used  with  the  values  of  —  or  determined  at  the  appropriate 

points. 

If  the  overall  point  efficiencies  are  required,  these  are  determined  also  from 


(6.3.1)  and  (6.3.2),  putting  r=f=0. 
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Distribution  of  solids  in  filtration 


FIGUSF  70 


Identification  of  Origin  of  Contaminant  Particles 


^.l.l  Introduction 

Particles  presented  in  hydraulic  oils  originate  from: 
j  "Built-in"  dirt.  This  is  swarf  from  the  initial  machining  of  the  components, 
that  generated  during  assembly  of  the  system  and  dirt  of  external  origin 
present  in  components  before  assembly.  It  is  recognised  that  any  frictional 
contact  will  give  rise  to  particles  and  for  this  reason  very  critical 
components,  e.g.  propellant  valves,  are  constructed  with  no  screw  threads 
within  the  system.  Although  this  dirt  is  usually  soon  removed  from  the 
system  by  the  internal  filters  it  may  do  great  damage  during  commissioning 
if  it  occurs  between  a  filter  and  a  sensitive  component. 

(b)  Extraneous  dirt.  This  will  enter  the  system  through  breathers-,  wiper 
seals  and  be  picked  up  on  the  disconnects  of  both  the  system  rnd  its 
servicing  equipment.  It  is  characteristic  of  the  environment  in  which  the 
system  is  operating  and  is  usually  silicacaeous  nature  with  some  fragments 
of  cloth,  paint,  etc.  from  servicing  operations. 

(c)  System  generated  dirt.  This  may  be  metal  particles  arising  from  component 
wear  or  elastomer  fragments  from  O-rirgand  seals. 

7.1.2  Identification  of  Particles 

Specifying  the  origin  of  particles  is  facilitated  by  identifying  the 
materials  of  which  they  are  composed.  For  example  silica  is  likely  to  be  extra- 
veous,  bronze  from  a  pump  or  bearing  and  nitrile  rubber  from  an  elastomer  seal. 
Microscopic  examination  on  membrane  filters  under  the  reflected  light  conditions 
specified  in  ARP  598  gives  some  information,  e.g.  bronzes  can  be  differentiated 
from  white  metals,  transparent  particles  may  either  be  silica  or  shreds  of  '*Teflon", 
but  for  example  two  different  types  of  bronze  or  aluminium  alloy  could  not  be 
*  di f f erent iated . 

Vickers,  Inc.,  Troy,  Mi.  (personal  communication  -  R.  Leslie)  have  been 
successfully  using  emission  spectroscopy  to  identify  the  elemental  composition  of 
particles  of  >  40pm.  However  as  particles  of  much  smaller  dimensions  are  of  great 
interest  in  hydraulic  oil  contamination  it  was  decided  to  investigate  the 
possibilities  of  electron  probe  x-ray  analysis. 
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7.1.3  Electron  Probe  X-ray  Analysis 

In  this  technique  of  analysis  the  specimen  to  be  examined  is  placed  in  a 
vacuum  environment  and  bombarded  with  high  energy  (10-30  KV)  electrons  causing  the 
sample  to  emit  x-radiation.  The  x-rays  produced  are  characteristic  of  the  elements 
present  and  may  be  differentiated  using  an  x-ray  monochromater  and  detector.  Many 
instruments  of  this  type  require  a  flat  specimen  surface  which  constrains  the 
geometry  of  the  system.  However,  this  is  not  possible  with  particles  and  hence 
the  instrument  chosen  was  the  Stereoscan  Scanning  Electron  Miscroscope  Model  2A 
with  x-ray  accessories  (Cambridge  Scientific  Instrument  Company,  Cambridge,  England). 
Fig.  72  is  a  conventiona  -canning  electron  micrograph  demonstrating  the  topography 
of  the  particles  and  figs.  73  to  76  the  scans  for  Al,  Fe,  Si  and  Cr.  No  Cu  was 
found.  This  is  a  preliminary  report  and  improvement  in  the  specimen  preparation 
technique  should  give  better  particle  distribution  in  the  field.  With  calibration 
quantitative  analyses  would  be  possible.  However  these  results  do  show  that  it  is 
possible  to  identify  at  least  the  major  elements  present  in  ly  particles. 

7.1.4  Technique 

In  electron  probe  analysis  it  is  necessary  to  coat  non-conducting  particles 
with  a  conducting  layer  (usually  carbon  or  Au/Pd  alloy)  to  enable  the  induced  charge 
to  leak  away.  An  attempt  was  made  to  examine  particles  retained  on  a  Millipore 
filter  by  mounting  the  filter  onto  a  specimen  mounting  stub,  coating  wiln  Ai/Pd 
and  examining.  However  the  particles  still  charged  and  the  cellulose  acetate 
filter  disintegrated  in  the  high  beam  current  necessary  to  induce  x-ray  emission. 

The  next  technique  tried  was  an  attempt  to  overcome  both  of  these  problems. 

A  pj.ece  of  silver  foil  (the  specimen  mount  is  aluminium  and  this  was  being  analysed 
for)  was  coated  in  'Silver  Dag',  a  colloidal  silver  suspension,  the  membraae  pressed 
particle  side  down  onto  it  and  when  dry  the  membrane  dissolved  away  with  acetone. 
Experience  however  showed  it  to  be  very  difficult  to  satisfactorily  remove  the 
membrane  without  detaching  the  colloidal  silver  from  the  foil. 

An  alternative  technique  tried  was  to  remove  the  particles  from  the  oil  by 
centrifuging,  wash  in  hexane  using  a  centrifuge,  disperse  ultrasonically,  pour  a 
drop  of  the  suspension  onto  the  silver  foil,  mount,  coat  and  examine.  Figs. 72 
76  were  prepared  by  this  technique.  It  has  the  slight  disadvantage  over  the 
Millipore  technique  in  that  examination  by  optical  microscopy  is  not  so  easy. 
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7.2 


Techniques  for  measuring  particle  and  fluid  behaviour  in  a  flowing  suspension 


7.2.1  Summary 

A  technique  has  been  developed  whereby  all  the  important  parameters  in 
suspension  flow  can  be  measured  for  each  phase  separately  without  any  disturbance 
of  the  flow.  The  fluid  velocity  profile  is  measured  with  a  technique  using  a 
photochromic  substance  dissolved  in  the  fluid.  The  particle  velocity  is  measured 
by  filming  the  motion  of  coloured  particles  introduced  into  the  bulk  of  glass 
particles,  which  have  the  same  refractive  index  as  the  liquid  and  thus  cannot  be 
seen.  The  concentration  profile  is  measured  by  a  standard  x-ray  attenuation 
technique. 

7.2.2  introduction 

Existing  methods  of  measuring  the  motion  of  the  two  phases  in  a  flowing 
suspension  have  two  major  disadvantages.  Firstly  they  involve  disturbing  the  flow 
in  one  way  or  another  and  secondly,  and  more  important,  it  is  often  not  possible 
to  measure  the  behaviour  of  the  two  phases  separately.  Attempts  to  measure  the 
flui-i  velocity  with  a  pitot  tube  for  example  2iave  tended  to  give  too  high  values 
caused  by  the  deceleration  of  the  particles  in  the  stagnation  region  just  in  front 
of  the  probe.  In  the  case  of  the  pitot  tube  there  is  also  the  problem  of 
particles  clogging  the  probe  which  will  inhibit  its  use,  particularly  for  measure¬ 
ments  of  fine  suspensions.  The  hot  wire  anemometer  is  another  example  where  the 
presence  of  the  solids  largely  affects  the  measurements.  A  different  class  of 
techniques  has  involved  colouring  the  liquid  by  injection  of  dye,  but  this  again 
disturbs  the  flow  and  the  particles  may  obstruct  the  coloured  liquid. 

To  overcome  these  difficulties  we  have  developed  a  new  method  whereby  all 
the  important  parameters  in  suspension  flow  can  be  measured  for  each  phase  separately 
without  any  disturbance  of  the  flow. 

A  number  of  aromatic  compounds  are  known  to  be  photochromic,  i.e.  they 
change  colour  when  irradiated  with  light  of  certain  wave  lengths.  The  change  of 
colour  being  instantaneous  enables  one  to  dye  certain  well  defined  parts  of  a 
liquid  in  motion  without  disturbing  it.  The  coloured  liquid  can  subsequently  be 
traced  photographically  or  by  inspection. 

Using  a  suspension  of  transparent  particles  in  a  liquid  of  the  same 
refractive  index,  the  irradiating  beam  will  go  through  the  suspension  without 
deviating,  leaving  all  fluid  elements  coloured  along  the  path  and  leaving  the 
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particles  unchanged.  The  matching  of  refractive  indices  and  using  some  coloured 
particles  enables  one  to  follow  the  translation  and  rotation  of  particles 
irrespective  of  their  position  in  the  cross  section.  The  concentration  of 
particles  can  be  measured  by  x-ray  attenuation.  The  attenuation  is  a  function  of 
the  total  mass  transversed  by  the  x-ray. 

7.2.3  Photochromic  Technique 

The  object  was  to  find  a  photochromic  substance  giving  rise  to  an  intense 
colour  persisting  long  enough  for  our  purposes  in  a  liquid  which  not  only  was 
suitable  from  this  point  of  view  but  also  had  the  same  refractive  index  as  some 
suitable  transparent  solid.  Literature  studies  indicated  that  DKBP*  would  be  a 
suitable  compound.  It  is  colourless  in  the  stable  lorm  (except  in  solutions 
with  high  pH  values),  and  the  unstable  isomer  is  intensely  blue.  There  is  some 
uncertainty  as  to  the  reaction  that  actually  takes  place.  The  following  two 
proposals  have  been  made: 


(2,4  dinitrobenzyl)  2  -  (2,4  dinitrobenzyl)  pyridine 


Whichever  reaction  that  actually  takes  place  the  coloured. form  is  going  to  he 
easily  oxidized,  and  therefore  the  solutions  used  may  be  expected  to  deteriorate 
when  exposed  to  both  light  and  oxygen.  This  has  in  fact  proven  true  in  our 
experiments. 

The  rate  of  the  fading  reaction  has  been  investigated  previously  by  others. 
It  is  a  first  oroer  reaction  with  the  reaction  constant  k.  Since  the  reaction 
involves  the  transfer  of  a  proton  it  follows  that  the  pH  value  should  be  of 
importance  in  a  water  solution.  Higher  pH  values  slow  down  the  reaction  (2).  Too 
high  pH  values  will  transform  the  compound  into  its  anionic  darn,  gieen  form  however, 
in  which  case  no  change  in  colour  will  take  place.  The  importance  of  the  nature  of 
the  solvent  and  the  temperature  is  shown  in  the  following  table  XX. 
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TABLE  XX  Solvent  conditions  for  photochromic  dyes 


The  absorption  spectrum  of  the  stable  form  shows  a  peak  in  the  region  from 
3000  X  slightly  into  the  visible  region.  There  are  other  peaks  at  smaller  wave¬ 
lengths  but  the  energy  involved  is  then  so  high  that  it  may  decompose  the  compound. 
Thus  the  light  must  be  filtered  so  as  not  to  transmit  any  large  amounts  below 
3000  X. 

The  visible  absorption  spectrum  of  the  unstable  blue  form  is  quoted  in 
literature.  There  is  a  peak  around  5500  X,  hence  the  blue  colour.  There  is  also 
an  indication  of  a  peak  in  the  same  region  as  that  of  the  stable  form,  i.e.  3000- 
4000  8.  This  indicates  that  the  penetration  of  a  UV-  neam  into  the  liquid  may  be 
rather  limited  and  that  there  will  be  an  optimum  concentration  of  DNBP  depending 
on  the  required  depth  of  penetration  and  the  intensity  of  the  light  source. 

As  described  below  we  decided  to  use  pyrex  and  a  toluene-ethanol  mixture. 
The  concentration  was  chosen  with  guidance  from  literature  to  0,1%  by  weight. 


Penetration  and  persistance  experiments  were  made  in  room  temperature  with  slurry 
and  pure  liquid  in  a  pyrex  tube  of  dimensions  to  be  used  in  later  experiments.  The 
internal  diameter  was  40  mm  and  the  wall  thickness  was  2.5  mm.  The  light  sources 
were  different  flash  tubes  emitting  200-400  joules  in  ca  1  millisecond.  The  wave¬ 
length  of  the  light  varied  over  a  wide  range  from  far  down  in  the  UV-region  up  in 
the  infra-red.  The  percentage  in  the  desired  region  (3000-400oX)is  not  known. 

The  blue  colour  penetrated  into  the  middle  of  the  tube  being  less  intense 
the  further  av/ay  from  the  light  source.  It  penetrated  approximately  the  same 
length  with  or  without  solids  present.  A  large  piece  of  pyrex  glass  inserted 
between  the  light  source  and  the  tube  had  no  measureable  effect.  The  results 
seem  somewhat  contradicting  but  the  accuracy  of  measurement  was  probably  not  good 
enough.  Photographic  recording  of  the  results  will  be  made  which  will  improve 
the  accuracy  and  also  enable  measurements  to  be  maae  very  shortly  after 
irradiation.  The  colour  persisted  from  one  to  several  seconds  depending  on  the 
initial  colour  intensity. 

Attempts  are  being  made  to  collect  the  UV-light  by  the  use  of  reflectors. 
Since  UV  is  absorbed  in  ordinary  glass  lenses  one  is  limited  to  the  use  of  MgFg- 
coated  reflectors  or  quartz  lenses. 

Preferably  the  light  should  be  collected  in  the  form  of  one  intense  beam 
so  as  to  enable  one  to  draw  a  pencil  line  through  the  fluid  under  investigation. 
Alternatively  one  can  use  a  cylindrical  reflector  in  conjunction  with  a  line 
source  whicn  will  produce  line  image  in  the  fluid. 

A  pulse  laser  would  solve  the  problem  of  collecting  the  light.  However, 
the  required  energy  per  pulse  is  larger  than  that  which  the  common  lasers  provide. 
The  laser  in  the  department  gives  pulses  of  about  1  joule.  The  frequency  would 
have  to  be  doubled  for  the  purposes  which  would  decrease  the  energy  to  about 
0,05  joule  per  pulse. 

7.2.4  Matching  of  Refractive  Indices 

The  solid  must  be  chosen  so  as  to  have  the  same  refractive  index  as  a 
suitable  liquid  and  at  the  same  tia*  transmit  light  of  3000  X  -  4000  %.  A 
transmission  spectrum  was  made  for  pyrex.  The  absorption  is  approximately  constant 
through  the  visible  region  down  to  3200  X  where  it  increases  rapidly.  Practically 
nothing  below  3000  %  is  transmitted.  Using  pyrex  glassware  will  thus  serve  as  a 
filter  for  the  non  wanted  lower  wavelengths. 
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Glass  spheres  are  available  as  load  glass  (ballatini)  or  as  soda  glass. 

The  ballatini  spheres  could  not  be  made  to  vanish  due  to  the  refractive  index  not 
being  constant  throughout  the  beads.  Etching  with  hydrofluoric  acid  did  not 
increase  the  transparency.  Sodium  glass  beads  are  only  available  in  a  very  poor 
quality  and  would  be  even  harder  to  make  vanish.  Both  the  lead  and  the  soda  glass 
would  in  addition  absorb  large  amounts  of  UV  light  and  would  be  unsuitable  in  high 
concentrations  for  that  reason.  Pyrex  glass  was  broken  into  particles  and  vanished 
sufficiently  as  particles  between  20Cp  and  30Qi  (by  sieving)  in  a  40  mm  tune  at  a 
concentration  ot  about  50%  by  volume. 

Since  the  photochromic  technique  works  well  with  ethanol  (n  =  1,3624)  the 
aim  was  to  find  a  liquid  with  the  refractive  index  higher  than  that  of  the  glass 
(n  =.  1,4700)  so  that  a  fair  amount  of  ethanol  could  be  added  to  adjust  the  refractive 
index.  Benzene  and  its  derivatives  are  the  only  solvents  with  high  enough  refractive 
index  with  the  exception  of  some  rather  toxic  and  unpleasant  bromocompound .  They 
have  refractive  indices  around  1.5.  Mixtures  containing  either  benzene,  toluene, 
xylene  or  benzyle  alcohol  were  tested  for  the  photochromic  reaction.  They  were 
equally  good  and  as  good  as  pure  ethanol  with  the  exception  of  the  last  one  which 
did  not  work  at  ell.  From  the  point  of  view  of  safety  we  chose  to  use  toluene. 


The  proportions  of  the  slurry: 

15,5%  by  volume  of  ethanol) 
84,5% . toluene) 

pyrex 


n  =  1,470 


n  =  1,470 


7.2.5  Concentration  measurement  by  x-ray  attenuation 

Concentration  measurement  by  x-ray  attenuation  is  not  a  new  technique  and 
therefore  it  will  only  be  mentioned  briefly  here. 

The  x-ray  source  is  chosen  so  that  as  large  a  change  in  counts  per  minute 
as  possible  is  obtained  for  a  certain  percenta3e  change  in  concentration  of  solid. 
The  attenuation  follows  the  Lambert-  Beer  law  and  therefore  this  occurs  for  the 
value  of  Jo/J  =  e.  The  fact  that  the  glass  walls  also  cause  attenuation  complicates 
matters  somewhat.  Since  both  pathlength  and  concentration  will  vary  in  these 
experiments  anyway,  one  can  disregard  the  effect  of  the  wall  if  it  is  comparatively 
thin  and  one  does  not  transverse  the  tube  toe  far  from  the  centre. 
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A  =  sensitivity  =  (£i)  (  ) 

Jo  ZB 

Jo=  number  of  counts  per  unit  time  with 
no  attenuation 

J  =  number  of  counts  per  k  unit  time 
m  =  mass  in  the  path  of  the  beam 


The  figure  illustrates  how  the  sensitivity  changes  with  the  mass.  If  the  attenuation 

follows  the  equation 
-km 

J  =  Jo  e 

one  can  deduce  the  following  equation  for  the  sensitivity 

.  ~  -tan 

A  =  K  m  e 

In  the  present  case  we  have  chosen  Americium  241  which  is  a  fairly  low  energy 
emitter.  The  path  through  the  centre  of  the  pipe  will  contain  a  mass  varying  with 
the  concentration  in  a  way  whereby  we  will  operate  on  both  sides  of  and  on  the 
peak.  The  recommended  equivalent  length  of  Aluminium  is  20-30  mm  and  the  actual 
equivalent  length  is  20-35  mm  if  the  concentration  varies  between  0  and  50%  by 
volume.  If  shorter  pathlengths  are  chosen  the  sensibility  will  be  high  only  if 
the  concentration  is  high.  In  the  actual  experiment  with  horizontal  pipes  one 
can  expect  concentration  to  be  low  in  the  upper  parts  of  the  pipe  where  the  path 
length  is  short  which  is  a  less  favourable  situation.  In  the  middle  and  lower 
parts  of  the  pipe  the  sensitivity  will  be  around  its  peak  however. 

7.2.6  Particle  translation  ar-*.  rotation 

Individual  particles  may  be  traced  by  introduction  of  coloured  particles 
and  photographing  them  with  a  higher  speed  film  camera.  Segregation  in  size 
fractions  etc.  may  also  be  detected. 
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Summary  and  Conclusions 


This  report  begins  with  an  analysis  of  the  present  state  of  the  art  with 
respect  to  particulate  contamination  of  hydraulic  flight  control  sub-systems.  This 
analysis  is  then  used  to  propose  an  overall  strategem  within  which  work  should 
proceed.  Because  of  the  complexity  of  both  particle  properties  and  hydraulic 
systems,  progress  in  this  direction  will  inevitably  be  slow  but  this  is  the  only 
long  term  method  likely  to  lead  to  a  priori  predictions  of  behaviour.  This  approach 
does  not  invalidate  the  continued  use  of  more  empirical  forms  of  study  which  enable 
specific  problems  to  be  solved,  but  the  latter  technique  is  in  the  long  term  waste¬ 
ful  of  effort. 

Until  the  necessary  techniques  for  describing  of  particles  with  respect  to 
their  bulk  behaviour  can  be  established,  the  methods  of  analysing  system  behaviour 
or  studying  component  behaviour  must  be  capable  of  application  to  any  future 
description  of  particle  properties. 

This  point  is  exemplified  in  Chapter  3,  Mathematical  Modelling  of  Hydraulic 
Systems,  where  a  stepwise  model  is  established  and  in  which  both  the  contaminant 
concentrations  within  the  system  and  the  characteristics  of  the  components  are 
expressed  in  a  general  matrix  form  which  is  applicable  to  any  method  of  particle 
size  analysis.  This  modelling  technique,  given  sufficient  information  about  the 
components  in  the  system)  enables  the  point  at  which  any  position  in  the  system 
reaches  a  critical  contamination  level  to  be  calculated. 

In  order  to  study  experimentally  both  operational  and  laboratory  systems  it 

is  necessary  to  be  able  to  remove  a  representative  sample  of  fluid  and  analyse  that 

& 

sample.  Chapter  2  states  the  requirements  of  a  sampling  device  and  describes 
techniques  for  both  the  theoretical  and  experimental  study  of  sampling  devices.  As 
far  as  can  be  established  none  of  the  sampling  devices  currently  used  meet  all  of 
these  requirements,  but  the  magnitude  of  the  errors  involved  are  not  known. 

Chapters  4  and  5  discuss  the  meaning  of  particle  size,  propose  a  new 
description  of  particle  size  and  give  accounts  of  the  study  of  three  well  known 
(e.g. ,  the  Coulter  Counter,  the  HIAC  Counter,  and  ARP  598)  and  one  more  specialised 
technique  of  particle  size  analysis.  Specifically,  the  assumption  of  equal  dist¬ 
ribution  of  particles  on  the  membrane  in  an  ARP  598  analysis  is  confirmed  to  be 
valid  for  one  size  of  particle  (and  would  not  be  expected  to  differ  for  other  sizes 
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